RAYLEIGH-TYPE SURFACE QUASIMODES 
IN GENERAL LINEAR ELASTICITY 



SONKE HANSEN 

Abstract. Rayleigh-type surface waves correspond to the charac- 
teristic variety, in the elliptic boundary region, of the displacement- 
to-traction map. In this paper, surface quasimodes are constructed 
for the reduced elastic wave equation, anisotropic in general, with 
traction-free boundary. Assuming a global variant of a condition 
of Barnett and Lothe, the construction is reduced to an eigenvalue 
problem for a selfadjoint scalar first order pseudo-differential oper- 
ator on the boundary The principal and the subprincipal symbol 
of this operator are computed. The formula for the subprincipal 
symbol seems to be new even in the isotropic case. 



1. Introduction 

Rayleigh [R] discovered the existence of surface waves which propa- 
gate along a traction-free flat boundary of an isotropic elastic body and 
which decay exponentially into the interior. The propagation speed of 
the surface wave is strictly less than that of body waves. Barnett and 
Lothe [LB 76] showed that Rayleigh-type surface waves can also exist 
at flat boundaries of anisotropic elastic media. 

The goal of this paper is to construct, for elastic media which are 
not necessarily isotropic, Rayleigh-type surface quasimodes which are 
asymptotic to eigenvalues or resonances. We use a geometric version 
of semiclassical microlocal analysis. 

The Rayleigh wave phemomenon of isotropic elastodynamics was 
explained by Taylor |Tay79| as propagation of singularities, over the 
elliptic boundary region, for the Neumann (displacement-to-traction) 
operator. Nakamura [Nak91j generalized this to anisotropic media, 
using the theory of Barnett and Lothe. Assuming isotropy of the elas- 
tic medium, Cardoso-Popov |CP92j and Stefanov [SteOOj constructed 
Rayleigh quasimodes. 
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Let (M, g) be an oriented Riemannian manifold with non-empty com- 
pact smooth boundary X. The (infinitesimal) displacement of an elas- 
tic medium occupying M is a vector field u on M. The Lie derivative of 
the metric tensor is a symmetric tensor field, Def u = C u g/2, called the 
deformation (strain) tensor caused by the displacement u. The elastic 
properties are defined by the elasticity (stiffness) tensor. This is a real 
fourth order tensor field C G C°°(M; End(T°' 2 M)), e \-> Ce, which 
maps into symmetric tensors and vanishes on antisymmetric tensors. 
We assume positive definiteness of C, i.e., (e\f)c = (Ce\f) defines an 
inner product on the space of symmetric tensors e and /. Here (-|-) 
denotes the inner product on tensors induced from g. This assumption 
is often called the strong convexity condition. If coordinates are given, 
then the components of C satisfy symmetries, C 1 ^ = & = C i , 
and eijCkt > if is a nonzero symmetric tensor. (We use the 
summation convention.) Denote the Riemannian volume elements on 
M and X by dV^ and dVx, respectively. The elasticity operator L 
and the traction T are defined, on compactly supported vector fields, 
by 

(1) / (Def u\ Defy) c dV M = [ (Lu\v)dV M + [ (Tu\v)dV x . 
Jm Jm Jx 

A positive mass density p G C°°(M) and the elasticity tensor C de- 
fine the material properties of the elastic medium. If the surface X is 
traction-free, then vibrations of the medium are solutions of the follow- 
ing eigenvalue problem: Lu = \ 2 pu in M, Tu = at X. See [MH83J 
for linear elasticity in the language of Riemannian geometry. 

The principal symbol of L, and of the /i-differential operator h 2 L, 
equals the acoustic tensor, c(£) = c(£, £) G End(CT x M), £ G T*M; 
see fl3Tl) . Here the associated acoustical tensor c(£, rj) G End(CT x M), 
£, rj G T*M, is defined as follows: 

(2) (c(£,t))v\w) — (y (8) r}\w <g) 

(Using g, we identify vectors with covectors.) The ik-th covariant com- 
ponent of c(£, 77) equals C^^rji. 

The existence of Rayleigh waves depends on the characteristic va- 
riety, S, of the surface impedance tensor, z. To define z, we first 
recall the definition of the elliptic boundary region, £ C T*X. Let v 
denote the unit exterior conormal field of the boundary X. Identify 
T*X = u ± C T* X M. By definition, £ G S iff c(f + sp) - p Id is positive 
definite for real s. From the factorization theory of selfadjoint matrix 
polynomials one gets g(£) G End(CT x M), £ G £ n T*X, such that 

(3) c(t + sv)-pld= (sld-q*(0)c(u)(sld-q(0), 



SURFACE QUASIMODES 



3 



s G C. Moreover, the spectrum of q(£) lies in the lower half-plane, 
specg(£) C C_, and these properties determine q(£) uniquely. The 
surface impedance tensor z is defined as follows: 

(4) z{Z) = ic(v)q(Z)+ic(v,Z), £e£. 

The significance of z results from the fact, proved in Lemma [181 that z 
is the principal symbol of a parametrix of the displacement-to-traction 
operator. In physics, the meaning of z is that it relates the ampli- 
tudes of displacements to the amplitudes of tractions (forces) needed 
to sustain these. 

The surface impedance tensor is Hermitian, and positive definite for 
large |£|, [BL85| Theorem 6]. If dimM = 3, then 

(U) z(£), £ e £■> nas a t m ost one non-positive eigenvalue. 

This property expresses the uniqueness of Rayleigh-type surface waves, 
[BL85j Theorem 8]. In case dimX ^ 3, we shall assume (jUJ) as a 
hypothesis. The characteristic variety of z, 

E = {£g£ ;detz(£) = 0}, 

is a smooth hypersurface, transversal to the radial directions of the 
fibers of T*X. Compare |BL85(, Theorem 7]. Rayleigh waves exist 
only if E is not empty. We shall make the stronger assumption that E 
intersects every radial line: 

(El) EnR + ^0 if £ G T*X \ 0. 

Compare [BL851 Theorem 12], [Nak91l Theorem 2.2], jKNOOl (ERW)]. 
Assuming flU} and flEll . there exists a unique p G C°°(T*X\0), p > 0, 
homogeneous of degree 1, such that 

(5) S=p- 1 (l). 

See Proposition [6j Furthermore, the kernel of z defines a line bundle, 
kerz -> S, over the compact base E. We shall require that its first 
Chern class vanishes: 

(E2) ker z ->■ E is a trivial line bundle. 

In particular, the bundle is assumed to possess a unit section. Property 
(IE2p is stable with respect to homotopies in the material properties; 
see Corollary [7J In the case of isotropic elasticity with positive Lame 
parameters, ([U|) . flElj) . and flE2j) hold. Moreover, 

E = {c r |£| = 1} c £ = {c s \£\ > 1}, p(£) = c r \C\. 

Here c r is the propagation speed of the Rayleigh surface wave which is 
strictly less than the speeds of the body waves, < c r < c s < c p . See 
Example [HJ 
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Next we state the central result of this paper: The traction-free 
surface eigenvalue problem can be intertwined with a selfadjoint eigen- 
value problem on the boundary. We employ a semiclassical pseudo- 
differential calculus, with distributions and operators depending on a 
small parameter, < h < 1. We write Ah = Bh iff the Schwartz kernel 
of Ah — Bh belongs to C°° with seminorms satisfying Oc°°{h°°). 

Theorem 1. Assume dimM = 3 ; or ([U]). Assume (1E1I) . (1E2I) . Given 
a unit section v o/ker z -> £, there exists a selfadjoint, elliptic operator 
P E \E fl (X; fi 1 / 2 ), independent of h, and operators, 

B h : L 2 (X; CT X M) -» L 2 (M; CTM), \\B h \\ = 0(h 1/2 ), 

J h ,J h E fi 1/2 , CT X M), J£j h e/Mzc at £, 

(/i 2 L -p)B h = 0, TS ft J ft = J fe (P - /T 1 ), 

and -B^lx = H in o neighbourhood of E. T/ie principal symbol of P 
equals p of (JSj). Furthermore, there is a formula, (1501) . /or t/ie subprin- 
cipal symbol p sn h of P. If v is changed to another unit section, e llfi v, 
then the subprincipal symbol changes to p su b + {p, <^} ; where {p, (p} de- 
notes the Poisson bracket. 

This result is known in the isotropic case, |CP92j . [SteOOj . except for 
the assertions about the subprincipal symbol. 

The operator Bh is a parametrix of the Dirichlet problem near S; 
see Proposition [T71 Its range consists of functions which are smooth 
in the interior of M, supported in a preassigned neighbourhood of the 
boundary, and which decay like e - Sdlst ( x )/ h j n to the interior. 

Ignoring finitely many eigenvalues the spectrum of P consists of a 
sequence of positive eigenvalues pj t oo. Applying Theorem [1] to an 
associated orthonormal system of eigenvectors we obtain, in Proposi- 
tion [211 a sequence of quasimode states: Luj — p 2 puj = Oc°°{hf) with 
boundary tractions equal to zero. Moreover, the quasimode states are 
well-separated. The construction also works when starting with a se- 
quence of almost orthogonal quasimode states of P. 

The unbounded operator D L 2 (M; CTM; p dV M ), u H> p^Lu, 
with domain D = {u 6 ; Tu = 0}, is symmetric and nonnegative. 
The associated quadratic form is given by the left-hand side of (jTJ). 
Denote Lt the Friedrichs extension of this operator. For a selfadjoint 
operator A with spectrum consisting of a sequence of eigenvalues ac- 
cumulating at +oo, denote iV^A) the usual counting function for the 
eigenvalues of A. The following lower bound on Nl t (X) is an example 
application of our results. 
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Corollary 2. Assume M compact, dimM = 3, and flEll) . ( 1E2I) . Lei 

P be the selfadjoint operator given in Theorem^ For every m > 1, 
Nl t {^) — Np(X — A~ m ) is bounded from below. 

Rayleigh waves have been studied in several papers with the empha- 
sis of getting information about resonances in scattering theory, [S V95j . 
[S V96j . |SjV97| , |Ste00j . and, for anisotropic media, |KN00] . Stefanov 
[SteOOj uses Rayleigh quasimodes to derive lower bounds on the num- 
ber of resonances. See the remark at the end of section [9] about going 
from quasimodes to resonances. 

The subprincipal symbol p su b affects the eigenvalue asymptotics of 
P, |DG75j . and it enters quasimode constructions, |CP92j . The sub- 
principal symbol occurs in the final formulas via integrals, such as 
f s , x Psub and J 7 Psub) where 7 is a closed bicharacteristic. We point 
out that these integrals do not depend on the choice of the unit section 
v in Theorem [H although p su t, itself does. An important aim of the 
present work is to give explicit formulas for the subprincipal symbol 
of P. These seem to be new even in the isotropic case which is dealt 
with in more detail in Proposition [25J The main difficulty comes from 
the fact that an invariant notion of subprincipal symbol has only been 
available for scalar operators. To overcome this obstacle we adapt and 
systematically use the geometric pseudo-differential calculus of Shara- 
futdinov [Sh a05al ISha05bj which assumes given a differential geometric 
structure. The principal and subprincipal symbol levels are contained 
in the leading symbol of a (pseudo-) differential operator. 

The paper is organized as follows. In section [2] the surface impedance 
tensor is studied; in particular, a selfcontained treatment of Barnett- 
Lothe theory is given. The leading geometric symbols of some differen- 
tial operators are computed in section [31 In section H] we geometrically 
decompose the elasticity operator near the boundary into normal and 
tangential operators, keeping track of leading geometric symbols. Sec- 
tion [5] gives, microlocally at the elliptic region S, a factorization of 
h 2 L — p into a product of first order operators. Using the factorization, 
we construct in section O a parametrix for the Dirichlet problem mi- 
crolocally at £ . The displacement-to-traction operator Z is defined in 
section and its leading geometric symbol is determined. In section [S] 
we derive a diagonalization of Z, and we prove Theorem [TJ In section M 
we construct localized traction-free surface quasimodes, and we prove 
Corollary [2j In section [10] we calculate, for an isotropic elastic medium, 
the subprincipal symbol of P. The appendix [A] contains a detailed ex- 
position of Sharafutdinov's geometric pseudo-differential calculus in a 
semiclassical setting. 
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2. The Surface Impedance Tensor 

First we collect some well-known facts about spectral factorizations 
of selfadjoint matrix polynomials. Refer to [ GLR8 2, Chapter 11]. Let 
V be a finite-dimensional complex Hilbert space, and f(s) = as 2 + bs + 
c G End(V) a quadratic polynomial in the complex variable s. The 
spectrum of / is the set of s G C such that ker/(s) ^ 0. Assume 
that the leading coefficient of /, a, is nonsingular. Then the spectrum 
is finite. Assume that / is selfadjoint, f(s)* = f(s), and that, in 
addition f(s) is positive definite for real s. The spectrum of / is a 
disjoint union <j + U cr_, where a + and cr_ are contained in the upper 
and lower half-planes, respectively. There is a unique q G End(V) such 
that f(s) — (s — q*)a(s — q), and the spectrum of q equals <r_. If 7 is 
a closed Jordan curve which contains cr_ in its interior and a + in its 
exterior, then 

(6) ql f(s)- 1 ds= (fsfisy'ds. 

The integral on the left is nonsingular. Jordan-Keldysh chains are a 
means to compute q. In particular, one has qv = sv if f(s)v = and 
Ims < 0. Moreover, the solvency equation f(q) = holds. 

The following representation of the factor q by integrals is important. 
We shall also apply it later to establish symbol properties. Denote 
i = a/— 1 the imaginary unit. 

Lemma 3. Let f and q be as above. Then 

(7) aqfo = -7rzld+/i, 

where fo = /(s) -1 ds is selfadjoint and positive definite, and 

h= f saf(s)- 1 ds+ [ s-^a-ZOO)/^)" 1 ^. 

J\s\<l J\s\>l 

The integrals converge absolutely in End(V). 

Proof. Let 7^ denote the negatively oriented closed contour composed 
of the semicircle {\s\ = R,lms < 0} and the interval [— R, R]. The 
integral representation ([6]) holds with 7 = 7^ if i? is sufficiently large. 
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We have /(s)" 1 = s~ 2 a~ l + 0(|s|~ 3 ) as \s\ — > oo. It follows that 



Urn j> f(s)' 1 ds = J™ f{s)- x ds, 



and 



lim d> saf(s) 1 ds = — 7ri Id + lim / saf(s) 1 ds. 
R ~*°° Jir R ^°° J-R 

Using s 2 af(s)~ 1 — Id = (s 2 a — f(s))f(s)~ 1 we obtain 

/ saf(s)- 1 ds= [ s -i( s 2 a _ /(s) ) /(s) -i ds . 

Jl<\s\<R Jl<\s\<R 

This proves the formulas. The remaining assertions follow from these 
and the positive defmiteness of f(s). □ 

Let £ G T*X, and denote v G T*M the unit exterior normal. Set 
a = c(v), cti(0 = c(z/, £), and a 2 (£) = c(£). Note that ai(£)* = c^;^)- 
The polynomial 

(8) /(s) = c(£ + sz/) - p = as 2 + (ai + a*)s + a 2 - p, 

/(s) = /(s,£), has values in End(CT.j.M). It is selfadjoint with real 
coefficients. By definition, £ G £ iff /(s) is positive definite for sel. 

Lemma 4. TTie elliptic region £ is an open subset ofT*X with compact 
complement. Moreover, £ is symmetric and star shaped with respect to 
infinity, i.e., G £ whenever £ G £ and t real, \t\ > 1. 

Proof. By positive defmiteness of C, there exists 5 > such that 
g(v,c(rj)v) > 5\v®7] + ri®v\ 2 for (co-)vectors ^,77. The symmetrization 
of a non-zero real elementary tensor is non-zero. Therefore, with a new 
5 > 0, in the sense of selfadjoint maps, 

(9) c{ V ) > 5\r]\ 2 Id. 

Since |£ + sz/| 2 = |£| 2 + s 2 the first assertions follow. The symmetry 
and the star-shapedness follow from c(trj) = t 2 c(r]). □ 

If £ G £, then ([7j) holds with q = g(£), = /_? (£)• The spectral 
factor g solves Q; using current notation: 

(10) as 2 + (ai + a{)s + a 2 — p = (s — q*)a(s — q). 

The spectrum of q lies in the lower halfplane, and q is uniquely deter- 
mined by these properties. Notice that q is a smooth section of the 
bundle vr* End(CT x M) -> £, where ?r : £ C T*X _> X denotes the 
canonical projection. 
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The surface impedance tensor, denned in (BJ, equals z = i(aq + ai). 
Lemma [3] implies 

(11) ^o = 7rId+z(A + aa/ ). 

Since the fj's are real, this gives the decomposition of z into real and 
imaginary parts. Following [MF04J, we shall use the Ricatti-type equa- 
tion 

(12) (z + ial)a~ 1 (z — ia\) = a 2 — p 

to deduce properties of z. Equation (|12p follows upon insertion of 
q = —a~ l {iz + ai) into the solvency equation associated with ffTU]) . 

(13) aq 2 + (ai + a\)q + a 2 - p = 0. 
A consequence of f|T2|) is 

(14) (iq)*z' + z'(iq) = <g + q*a[ + (a 2 - p)' + q*a'q, 

where the prime denotes the derivative with respect to some cho- 
sen parameter. The spectra of q and q* are disjoint. Therefore, the 
Sylvester equation {iq)*x + x{iq) = i{xq — q*x) = y has a unique so- 
lution x for given y. The solution is, in fact, given by an integral, 
x = J° exp(irq)*y exp(irq) dr. It follows that x is positive definite if 
V is- 

Proposition 5. The impedance tensor z(£), £ G £ , has the following 
properties, 
(i) z(£) is selfadjoint. 

(ii) z(£) is positive definite if |£| is sufficiently large. 
(Hi) Rez(^) is positive definite. 

(iv) z(£) has at least two positive eigenvalues if dim M > 3. 

(v) (d I ' dt)\ t= \t~ x z(t£) = z — z is positive definite. 

(vi) The complex conjugate z(£) = z(— £). 

We call £(£) = (d/dt)\ t= iz(t^) the radial derivative of z at £. It 
follows from (jvj) that i is positive definite on the kernel of z, ker z. 

Proof. To prove @ we follow the arguments in |MF04[ Theorem 2.2]. 
First note that ( |T2|) remains true if z is replaced by z* . Subtracting the 
two equations we get the Sylvester equation (iq)*(z— z*) + (z — z*){iq) = 
0, implying z — z* = 0. 

It follows from (fTTj) that Re z = 7t/q . This proves ( Imj) . 

Suppose dimM > 3. Aiming at an indirect proof of (jiv)l . assume that 
z(£), £ e h as a ^ m ost one positive eigenvalue. Then there exists 

w G CT X M such that z(£) is negative semidefinite on the orthogonal 
complement w 1 . Choose a real vector n ^ which is orthogonal to 
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both Rew and Imw. Then v G w 1 , and (Re z(£)v\v) = (z(£)v\v) < 0, 
contradicting the positive definiteness of Re z. 

Next we prove (jvj) following the method of [MF04, Theorem 2.3]. 
Since Oj(£) is homogeneous of degree j in £, equation (fr2|) implies, 

(r 1 ^) + iol(0)a- 1 (r 1 z(tO - iai(O) = a 2 (0 - t~ 2 P . 

Taking the derivative with respect to t at t = 1, we get 

— 2) + (i — = 2p. 

By the remarks following f JT4|) . we see that z — z is positive definite. 

We now prove Note /(s, -£) = f(s,£), /,•(-£) = ( -l) J 7j(fl , 
and Ol (-0 = -%(£). Using ([II]) we derive z(-£)/o(£) = *(0/o(0- 
Since /o is real and nonsingular the formula follows. 

It remains to prove (Jn|). Let 77 G T^X, |r/| = 1. It suffices to show that 
Zoo = li m t|oo t^zitrf) exists and is positive definite. Set q t = t^qitrf), 
t > 1 large. From ffTO]) deduce 

as 2 + (01(77) + a\(q))s + a 2 {rf) - t~ 2 p = (s — q* t )a(s - q t ), s G R. 

Using ([7]) and dominated convergence in the integrals giving fj we 
infer that q^ = lim^ q t exists. In particular, t~ l z(tr]) converges to 
Zoo = i(aqoo + a i(v)) as t t 00 • Let y G T X M such that (zoo2/||/) < 0. 
We must show y = 0. Set w(r) = exp(irg 00 )?/, r < 0. The solvency 
equation ( ITBl holds with g replaced by goo, p = 0. Therefore, aD 2 w + 
{ai+a\)D r w+a2W = holds, where we use the abbreviation a 3 - = 07(77). 
Take the inner product in CT X M with to and integrate. A partial 
integration gives 




(aD r w + aiw\D r w) + (D r w\aiw) + (a,2w\w) dr 



= i(aD r w + aiitj|iw)|_ = (zoayly) < 0. 

Set W(r) = w(r) <g> 77 + D r w(r) <S> v G End(CT x M). Recall a = c(u), 
ax = c(v,r)), a 2 = c(rj), and (j2J). We have shown: 

/ (W(r)\W(r)) c dr < 0. 
J —00 

Recall that C is real, and that ( | )q is an inner product on symmetric 
2-tensors. It follows that the symmetrization of W(r) vanishes for all 
r < 0. In particular, 

(15) (w(r) (g) i] + D r w(r) <g> v \ ( <g> v + v <g> C) = 

for C G CT^M, r < 0. Recall (77^) = 0. Setting ( = v, we derive 
-D r (w(r)|z/) = (D r w(r)|z/) = 0. Since w(r) — > as r — >■ —00, we obtain 
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(w(r)\u) = 0. Now, ( 1T5]) simplifies to (D r w(r)\() = 0. Since ( is 
arbitrary, this implies, successively, D r w = 0, w = 0, y = 0. □ 

If dimM = 3 then (JU|) holds. This follows from flrv]) . 

Proposition 6. Assume (jH])- T/ien t/ie characteristic variety of z, 
£ = {detz(£) = 0}, is a smooth hypersurface in £. Each radial line 
M + £ C T*X intersects £ in at most one point, and the intersection is 
transversal. The kernel of z\x defines a line bundlekerz -» S. Assume, 
in addition, (lEip . There is a unique p G C°°(T*X\0) , homogeneous of 
degree one, such that E = Moreover, p > 0, and p(— £) = 

Proof. From the assumption and ([v]) of Proposition [5] it follows that 
(d/dt) det z > if t£ G S, t > 0. In particular, zero is a regular 
value of det z. Hence £ is a codimension one submanifold transversal 
to the radial field. Since M + £ D £ is connected, a given radial line 
intersects S in at most one point. Because of (|U]) and the selfadjoint- 
ness of z, zero is simple eigenvalue of z. It follows that kerz -> E is a 
line bundle. Now assume also (lEip . Then each radial line intersects £ 
in a unique point. Define p as follows. For ^ £ G T*X set = 1/t 
if t£ G E, t > 0. Smoothness of p follows from the implicit function 
theorem. The evenness of p is a consequence of fjvTj) . The other proper- 
ties of p are obvious. Clearly, the homogeneity and p|g = 1 determine 
p uniquely. □ 

Corollary 7. Let p t and C t be homotopies of the mass densities and 
the elasticity tensors, < t < 1. Assume that the associated surface 
impedance tensors z t and their characteristic varieties £t satisfy ((HI) 
and (lElj) /or ever?/ 1. 27ie line bundles ker zo So one? ker zi -> Si are 
zsomorpazc. 

Proof. The factorization ()3]) and the definition of the impedance ten- 
sor imply that z t depends continuously on the homotopy parameter 
t. It follows from Proposition E] that the characteristic varieties are 
canonically diffeomorphic to the sphere bundle SX. We deduce that 
the Chern classes of the bundles ker z t -> SX do not depend on t. The 
assertion follows from this. □ 

Example 8. We consider, as special case, an isotropic elastic medium. 
We shall verify (JH]), f lEip . and (1E2j) . The elasticity tensor reads, in 
component notation, 

(16) C*" M = Aa V + H(9 ik 9 ji + 9 li 9 jk ), 
where A, fi denote the Lame parameters. Equivalently, 

(17) c(£,n) = Ae®n + /in®e + /io(e,")W. 
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Positive definiteness of C is equivalent to p > 0, AdimM + 2p > 0. 
We make the stronger assumption A,p > 0. Let £ G r*X. We list 
the eigenvalues s G C and the eigenvectors v G CT X M of the quadratic 
polynomial c(£ + si/) — p: 

(a) (A + 2p)(|£| 2 + s 2 ) - p = and u = £ + si/, 

(b) p(|£| 2 + s 2 ) - p = and v = s£ - 

(c) p(|£| 2 + s 2 ) — p = and v is orthogonal to £ and za 

Introduce c p = \/(X + 2p)/p and c s = \fpjp, the speeds of pressure 
and of shear waves, respectively. Assume that £ G £. This is equivalent 
to c s |£| > 1. The above eigenvalues and eigenvectors diagonalize q, 
q(£,)v = sv if Ims < 0. Denote V the subbundle of 7r*(CTxM) -> S 
spanned by v and £, and l^ -1 its orthogonal bundle. Fix the orthonormal 
frame u, £ = £/|£| of V, and choose an orthonormal frame of l/ -1 . In 
block decompositions of matrices we let the indices 1 and 2 correspond 
to V and l^- 1 , respectively. We denote (e)ij the block ij of the matrix 
which represents the endomorphism e. Observe that q leaves V and 
V 1 - invariant, (q)n — = (9)21- A simple computation gives 

lei 



(*?) 
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Uty/l — t —%{b — ut) 
i(b-t) ty/l - ut 

Here t = (c s |£|)~ 2 , u = (c s /c p ) 2 = p/(A + 2/i), 6 = 1- VI - ut^/T^t. 
Moreover, {iq)22 equals — £ times the unit matrix. The maps 

a = c(z/) and ai = c(i/, £) also leave V and l^ -1 invariant. We compute 



(19) (z)n 



ty/T^t -i (2b - 1) 
i(2b - t) ty/1 - ut 



and (2)22 = p (2(7)22- The determinant of z equals (p|£|\/l — t) dimV 
times 

(20) det(z) n = p 2 )^ 1 (V(l - t)(l - ut) - (2 - t) 2 ). 

Given k g]0, l/2[, the unique zero t g]0, 1[ is found as the solution of 
Rayleigh's cubic equation, (0 (24)], = ((t-2) 4 - 16(1 -t)(l - ut))/t. 
Define the Rayleigh wave speed c r = c s Vt G l7°°(X). Set p(f) = c r |£|. 
The characteristic variety £ equals {p(£) = 1}. Thus (jUJ) and (lElj) 
hold. Obviously, i(26- + - G kerz(£), £ G E. Observe that 

(21) 2(26 - t) = t(2 - t) on S. 
Thus 

(22) z(2-t)z/ + 2 v / l-t£Gker^, i = (c r /c s ) 2 , 

is a nowhere vanishing section of the kernel bundle. Hence also (1E2[) 
holds. This example is of course well-known. 
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Remark. The identity ( ITTT) goes back to Barnett and Lothe; compare 
[LB76[ (3.18)]. It is key to proving, in dimension three, the unique- 
ness of subsonic traction-free surface waves |BL85l Theorem 8]. The 
second assumption in Proposition [6] is needed to prove the existence 
of Rayleigh surface waves. Compare with [BL85, Theorem 12], where 
existence criteria are given in terms of the so-called limiting velocity 
which corresponds to the boundary of the elliptic region. See |Nak91t 
Theorem 2.2] for the Barnett-Lothe condition in a microlocal setting, 
and the real principal type property of the Lopatinski matrix it entails. 
See [Tan07] for a recent exposition of Barnett-Lothe theory, and for a 
treatment of isotropic and transversely isotropic media. 



3. Connections and Geometric Symbols 

The elasticity operator is defined in terms of the Levi-Civita con- 
nection and of the elasticity tensor. We use the geometric pseudo- 
differential calculus of Appendix [A] to define and compute the leading 
symbol of the elasticity operator. The leading symbol includes the prin- 
cipal and the subprincipal level. The calculus depends on the choice of 
connections. 

Equip M with the Levi-Civita connection of g. Let exp denote its 
exponential map. If x, y G M, then denote by [y <- x] the shortest 
geodesic segment from x to y, assuming its interior does not intersect 
the boundary, and that it is unique. 

Let E -y M be a (complex) vector bundle with connection V E . 
Denote G End(£ , a; , E y ) the parallel transport map along a given 
curve 7 in M from x to y, e.g., T ^ x y The connection can be recovered 
from its parallel transport maps: 

(23) Vfs(x) = -^\ t=Q T^ exPxtv] s(exp x tv). 

Denote n*E -> T*M the pullback of E to the cotangent bundle 
7r : T*M — > M. Let a be a smooth section of tt*E -> T*M. Following 
|Sha05a| [Sha05b] . we introduce the vertical and the horizontal covariant 



derivative of a. The vertical derivative v \/a(x, £) G E x ®T X M, at 
f G T*M, is the derivative of the map T*M ->■ E x , £ ^ a(x,£). The 
definition of the vertical derivative depends only on the linear structure 
of the fibers of T*M. The horizontal derivative h Va(x,£) G E X ®T*M 
is the derivative at v = of a map T X M — > E x , 

d 

(24) h Va(x, = — I _t£ ex „ v] a(exp x v, tLJ?^). 



dv 
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The horizontal derivative depends on the Riemannian structure and on 
the connection V B . In the scalar case, E — C, in local coordinates, 

h Va(x,0 = {d Xj a(x,0 +T%(x)Z k ds t a(x,Z)) dx* , 

where T^- denote the Christoffel symbols of the Levi-Civita connection. 
Writing a local section of n*E as a sum of products ai(x, £)a,2(x) where 
a\ is scalar and a 2 a section of E one readily derives local formulas 
for the horizontal derivative in terms of connection coefficients. The 
vertical and the horizontal derivative extend to first order differential 
operators, "V and ^V, which map sections of 7r*(E(g)T r ' s M) to sections 
of n*(E®T r+1 ' s M) and of 7i*(E®T r ' s+1 M), respectively. The operators 
"V and commute. It suffices to prove this when E is the trivial line 
bundle, E — C In this case the assertion is easily checked in normal 
coordinates. 

Let F -> M be another vector bundle. Let A : C°°(M; E) -> 
C°°(M; F) be a differential operator of order m. We introduce a small 
parameter, < h < 1, and we replace A by the /i-differential operator 
h m A. Then A G \I> m ' (M; E, F) as a semiclassical (pseudo-) differential 
operator. Refer to Appendix |A] for an exposition of Sharafutdinov's 
geometric pseudo-differential calculus in a semiclassical setting. The 
formula (167]) for the geometric symbol, <Jh{A) G S m '°, simplifies to 

(25) a h (A)(x,Os = M^' CXP ' ly)/ SU s )\ y = x > 

where £ e T*M, s G E x , and z = y/—l. The geometric symbol extends 
by continuity to the boundary of M. In symbol computations we track 
the leading geometric symbol, defined before Proposition [571 In the 
following, the symbol of an operator is always its geometric symbol. 

For the Laplace-Beltrami operator one has cr/ l (— h 2 A)(x, £) = |£| 2 . 
This is readily checked using normal coordinates. 

From ([25]) and flUD deduce 

(26) a h {-ihV E ){i)e = e®ieE x ®T*M. 

As before, to ease notation, we usually do not write the base point x 
into the arguments of tensors and symbols. 

If E -> M is a Hermitian vector bundle then we define, using the 
volume element <1Vm, the Hilbert space L 2 (E). Assume E and F are 
Hermitian vector bundles having metric connections. The leading sym- 
bol of the formal adjoint A* of A is given by 

(27) a h (A*) =a h (A)* -ihtT( v V h V(T h (Ay). 
See Proposition 1271 
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Equip the bundle E £g> T*M with the induced Hermitian structure 
and the induced connection. The connection is metric. Observe that 
the horizontal derivative of ah{—ihV E )* vanishes. Therefore, (1271) and 
flU imply 

(28) a h ((-ihV E y)(0(e®v)=9(Z,v)e, t,V G T^M, e e E x . 

By Proposition [28] the leading symbol of a composition is given as 
follows: 

(29) a h (AB)=a h (A)a h (B)-ihtr( v Va h (A). h Va h (B)). 

The trace is the contraction of the TM®T*M factor which is produced 
by a pair of vertical and horizontal derivatives. The dot terminates a 
differentiated expression, serving as a closing bracket. 

Let C e C°°{M;End(E ® T*M)). View C as an operator which 
acts by multiplication on sections of the bundle E ®T*M ->• M. Let 
V denote the connection on the bundle End(E £g> T*M) -> M induced 
from the Levi-Civita connection and from V s . Define sections c, dive 
of 7T* End(E) -> T*M as follows: 

c(Oe=(e,7r*C(e<8»0). 

(divc)(Oe = ^(^', (7r*V v .C)(e ® 0>, 
j 

where the angular brackets denote contractions on covectors, using g. 
Furthermore, (vj) and (rf) are any dual frames of TM and T*M. 

Lemma 9. o h (-h 2 V E * o C oV E ) = c - ihdiv c + 0(h 2 ). 

Proof. Observe that a h (C) = n*C, and h \7n*C = tt*VC. The symbol 
( 1281) is linear in £. Its vertical derivative is obvious. Using ( J29|) . the 
symbol a of —ihSJ E o C is found to be 

a(0(e®?7) = (e,7r*C(e®77))-i/iJ](77 J ',7r*V, ) .C(e®7 ? )). 

i 

Here (uy) and (r]- 7 ) are as in the definition of dive. The horizontal 
derivative of the symbol of —ihV E vanishes. Therefore, 

a h (-h 2 V E * o C o V E )(Oe = a(0a h (-zhV E )(0e = a(£)(e <g) 0, 

where we used ( 1291) . □ 

Now assume = CTM and C the elasticity tensor. Identify 

End(CT°' 2 M) = End(CTM ® CT*M). 
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Let L the elasticity operator defined in (TI]). Recall from Riemannian 
geometry the following relation between the Levi-Civita connection and 
the Lie derivative: 

(30) (£ u g)(v, w) = g(V v u, w) + g(v, V w u), 

for (real) vector fields u,v,w. Using the symmetries of the elasticity 
tensor we get 

L = Def oC o Def = (-iV)* o C o (-iV). 
We obtain the following corollary to Lemma [9j 

(31) a h (h 2 L~ p) = c - pld-ih&w c + 0(h 2 ). 

If respresent C with respect to some given local coordinates, then 
QUID reads 

a h (h 2 L - p)(0 lk = Cv%& - p5 lk - V^lh&^tt + 0{h 2 ). 

The vertical bar followed by j means covariant differentiation with 
respect to the j-th coordinate. If the elastic medium is isotropic then 
the leading symbol becomes 

*h(h 2 L - p)(0 = p{c 2 p \i\ 2 - 1)P(0 + p{c 2 s \i\ 2 - l)(Id ~P(0) 
-z/i(VA®£+(V/i®£)* + (£,V//)Id), 

where P(£) = | <8> | denotes the orthogonal projection to the propaga- 
tion direction £ = £/|£|. 

4. The Elasticity Operator in a Boundary Collar 

In a boundary collar, ] — e, 0] x X C M, we write the elasticity op- 
erator L in terms differential operators on X having coefficients which 
depend on r e /, the negative distance to X. 

Let N(x) G T X M denote the unit exterior normal at x G X. There 
exists e > such that, if we set J =] — e, 0], the exponential map of the 
Levi-Civita connection defines a diffeomorphism onto a neighbourhood 
of X in M: 

J x X — > M, (r, x) t-> y = exp(riV(x)). 

Essentially without losing generality, we assume that this map is onto 
M. The inverse map is y ^ (r,x), where — r = d(y,X) is the distance 
from y to X, and x = p(y) is the unique point in X closest to y. 
The distance function r satisfies the (eikonal) equation |Vr| = 1 in M. 
Extend X to M by X = Vr. Also introduce the unit conormal field 
v = dr. The level hypersurfaces 

M r = {yeM;r + d(y,X) = 0} 
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are diffeomorphic to X = Mo. The shape operator S = ViV is a 
field of symmetric endomorphisms of TM, g(Su,v) = g(u,Sv). The 
second fundamental forms of the level hypersurfaces M r assign (u, v) h-> 
—g(Su, v) (Weingarten equation). The dependency of the metric tensor 
on r is given by the formula {C^g){v,w) = 2g(Sv,w). This formula 
follows from (|30|) . Introduce J e C°°(I r x X), the solution of d r log J = 
trS, J\ r =o = 1. Then we have the following formula for the volume 
form of M: 

(33) / f{y)dV M (y)= [ [ f(exp(rN(x)))J(r,x)dV x (x)dr, 

J M J I J X 

f E Cp°(M). See |Pet98t Ch. 2] for the geometry of hypersurfaces using 
distance functions. 

Let E -> M be a vector bundle with connection V s . Denote E r -> M r 
the bundles induced by the inclusions M r C M, reJ. Set i?x = -So- 
Let u G C°°(M] E) be a section of i£. Using parallel transport in E 
along the geodesies which intersect the boundary orthogonally, define 
u:I^C°°(X;E x ), 

u( r )(x) = u(r,x) = T { ^ y] u(y), if y = exp(riV(x)). 

The map 

(34) C°°(M; E) — )■ C oc (J, C°°(X] E x )), u^u, 

is an isomorphism of Frechet spaces. The isomorphism commutes with 
bundle operations such as tensor products and contractions. 

The covariant derivative in normal direction is transformed into d r 
under the above isomorphism: 

(35) V|w(r) = d r u{r), r el. 

To see this, consider the geodesic / —> M, r \-t y(r) = exp(rN(x)). 
The tangent vectors are y(r) = N(y(r)). Using f l23|) . it follows that 

(V£ (v(r)) u)(y(r)) = ^\ s=r r§ ir ^y( s)] u{y{s)) 
j£ d I , ., > 

This implies (ESJ). We have V N N = SN = 0. It follows that d r N = 0, 
and d r v = 0. Abusing notation, we write d r to denote Vfr. 

Define e u e = e g) v and i v {e ® r\) = (rj, v)e. Notice that e v and i v 
commute with d r . 

Let F -> M be a another vector bundle with a connection. Let 
5 : C°°(M;E) -> C°°(M;F) be a differential operator. Assume that 
i? is tangential. This means, by definition, that i? commutes with the 
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distance function r, [B,r] = 0. Then, for every r £ I, B restricts 
to an operator B r : C°°(M r ;E r ) -»■ L7°°(M r ;F r ), B r U = (Bu)\ Mr , 
where u is a section of E -> M which extends a given section [/ of 
£7 r -> M r . The assumption [-B,r] = implies that B r is well-defined. 
Parallel transport along the geodesies orthogonal to X defines bundle 
isomorphisms E r = Ex and F r = Fx- Via these isomorphisms the B r 's 
induce differential operators B(r) : C°°(X;E X ) -> C°°(X;Fx-), called 
associated with B, such that Bu(r) = B(r)u(r), r G /. Each B(r) is 
a differential operator having coefficients which are C°° with respect 
to r. Conversely, an operator B is tangential if it is given in this way 
by a family of differential operators {B(r) ; r G /} with coefficients 
depending smoothly on r. 

Lemma 10. Let E M be a real vector bundle with connection V E . 
Then 

(36) V E = e u d r + B, 
where B is tangential. Moreover, B(0) = V Ex . 

Here E cg> T*M carries the induced connection. The lemma extends, 
by decomposition into real and imaginary parts, to complexifications of 
real bundles with connections. In particular, it holds for complexified 
tensor bundles with the Levi-Civita connection. 

Proof. Let P ± ,pW E C°°(M; End(TM)) denote the orthogonal pro- 
jectors onto the span of iV and onto its orthogonal complement, iV -1 , 
respectively. Identify E®T*M with Hom(TM, E). Let u e C°°(M; E). 
We have the following decomposition in C°°(M; Hom(TM, E))\ 

\/ E u = (V E u)P L + (V E u)P 11 = (V%u) ®v + Bu. 

This defines B, and implies fl36l) . Note that B is tangential. We have 

S(0)u(0) = Bu\ x = {{V E u)pW)\ x = (V Ex (u\ x )){P ll \x). 

This proves the asserted formula for -B(O). □ 

Assume E ->■ M a Hermitian bundle with a metric connection. Us- 
ing fl33l) . and the fact that parallel transport preserves inner products, 
we have 

(37) / ( u \v) E dV M = I I {u\v)E x JdV x dr, 

J M J I JX 

if u, v G C^°(M; E). Formal adjoints of differential operators on M 
are taken with respect to these inner products. The inner product 
of sections u and v of E x -> X is J x (u\v)e x dVx- Formal adjoints 
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of operators A(r) associated with a tangential operator A are defined 
with respect to this inner product. 

Next we prove a formula which expresses the elasticity operator L as 
a quadratic polynomial in D r = —id r with tangential coefficients. Now 
assume E = CTM, and let B as in (I36I) . Define tangential operators, 

A = L u Ce U} A 1 = -ii u CB, A 2 = B*CB. 

The order of Aj is j. Moreover, A\ = iB*Ce u . 

Proposition 11. The elasticity and traction operators defined in (CQ) 
are as follows: 

L — (D r — i tr S) (A D r + A x ) + A{D r + A 2 , 
-iT = A (0)D r + Ai(0). 
Furthermore, A*(0) = A^O)* . 

Proof. Let u,v £ C 6 °°(M; CTM). It follows from ([30]) and the symme- 
try properties of the elasticity tensor that 

f ( Def u | Def v) dV M = [ (CVu \ Vv) dV M - 

J M J M 

Inserting ( 1361) and using the definition of Aj, the right-hand side equals 

/ {i v CVu\d r v) dV M + [ {B*CVu\v)dV M 
Jm Jm 

(A d r u + iAiu\d r v) dVu + / (— iA\d r u + A 2 u\v ) dVu 
M Jm 



Partial integration with respect to r is done, using (1371) . as follows: 

/ (w\d r v)dV M = [ [ {w\d T v)J dV x dr 
Jm J i Jx 

= [ (w{0)\v(0)) dV x - [ [ ((d r \ogJ)w + d r w\v)JdV x dr. 
Jx J i Jx 

Summing up we have 

/ (Def« | Defv) c dV M 
Jm 

= / ((D r - itiS)(A D r + A x )u + A\D r u + A 2 u \ v) dV M 
Jm 

+ / (A o (0)(d r u)(0)+iA 1 (0)u(0)\v(0))dV x . 



' x 

Comparing with the formulas for L and T follow. The last assertion 
follows because J = 1 at X. □ 
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Next we compute the leading symbols of the operators (associated 
with) Ay The symbols are r-dependent sections of 7r*End(CTxM) -> 
T*X. Dropping tildes, the symbol of A equals 

a h (A ) = a = c{u) G C°°{I, C°°{T*X; tt* End(CT x M)). 

Introduce the divergence of the acoustic tensor restricted to X as fol- 
lows: 

(div x c)(C)v = (V a , {^ Va C){v ® C)>, 

if C ^ T*M, v G T X M, x G X. Here (v a ) and (j] a ) are any dual frames 
of TX and T*X. If local coordinates are chosen such that r is one 
coordinate and the other coordinates are constant along the geodesies 
orthogonal to X, then (divx c)(C)* fc = C l \ a Ct- Here the summation 
convention is used with Latin indices refering to all coordinates, and 
Greek refering to all coordinates except r. We also need the contraction 
(C,S) G C°°(M; End (TM)) of the elasticity tensor with the shape 
operator, in coordinates, 

//"> G\*^ riiikl o c 

\y,o) = u bj£, bj£ = Uj\£. 

(Because of VatS 1 = one can also write Greek coordinates instead of 
j and £.) 

Lemma 12. Let a% and ci2 denote the principal symbols of the h- 
differential operators hA\ and h 2 ' A2, respectively. At r = 0: ai(£) = 
c(z/, £), and 02(C) = c (0- On the leading symbol level, o~h{hAi) = ax, 
o~h{hA\) = a* — ihax-, and ah(h 2 A 2 ) = a 2 — iha 2 - + 0{h 2 ), where, at 
r = 0, 

ai - = (dw x c)(u) + 7r*(C,S), a 2 _(£) = (div* c)(£). 

Proof. By Lemma [TUl we have 

hAi(0) = l u C o (-iTiV), h 2 A 2 (0) = (-ihV)* 0C0 (-ihV), 

where V = V TX is the Levi-Civita connection of the boundary. We 
compute the leading symbol of hAi(0) using the composition formula 
Recall The vertical derivative of the symbol of i v C vanishes, 
Hence 

a h (hAx)(0)(0 = o- h (Mi(0))(0 = c(v,£), £ G T* X M. 

The formula for ah{h 2 A2(0)) follows from Lemma [91 In view of ff2Tj) . 
a\- = tr "V^VaJ. Since a\{£) = c(£,v) = (£, n*{Ce v )) is linear in £, 
its vertical derivative is immediate. Hence 
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Now, V v (Ce v ) equals {S7 v C)e v plus a contraction of C with V^z/ = Sv, 
proving the formula for ai_. □ 

If the elastic medium is isotropic, ( Tl6|) . then a straightforward com- 
putation shows that, at r = 0, 

(div x c)(C) = (VA <g> C) + (V/i ® 0* + (C, V/i) Id, 
(C, 5) = (A + //)£' + (/itrS) Id. 

Here VA, V/i G TX C TxM are the gradients of the Lame parameters 
restricted to X. 



5. Microlocal Factorization 

We factorize, microlocally in the elliptic region, the /i-differential 
operator h 2 L — p into a product with right factor hD r — Q, where Q is 
a tangential /i-pseudo-differential operator such that the spectrum of 
its principal symbol is contained in the lower halfplane, C_. 

As in the previous section we identify M with a boundary col- 
lar / x X, and sections of CTM -> M with r-dependent sections of 
CTxM -¥ X. Operators are polynomials in D r /h = hD r with tan- 
gential /i-(pseudo-)differential operators as coefficients. The latter are 
quantizations Bh = Op h (bh) G ^tang, of tangential symbols, 

h e SZn g = C°°{I, S m ' k {T*X; vr* End(CT x M))). 

By Proposition [Til the principal symbol f(s,£) = c(£ + sv) — p of 
h 2 L — p at £ + is a second order polynomial in s. View s as the 
symbol of D r /h- The coefficients are /i-independent tangential symbols. 
By (jnj), there exists a constant 5 > such that 

(38) /( S ,0>^(l + |s| 2 + |£| 2 )ld, s6K, 

holds if £ is sufficiently large. If F C £ is closed and R > 0, then 
F\ {|£| > -R} is compact. Hence there exist < e',5 such that (1551) 
holds uniformly for (r, £) G [— e', 0] x F. We say that a property holds 
at the elliptic region £ if it is true in every open subset of I x £ where 
fl38|) holds uniformly. 

Recall from section [2] that we have a unique spectral factorization 
([TO]) at £. 

Lemma 13. Lei q = £ G £ ; i/ie unique solution of the spectral 
factorization /(s, £) = (s — <7(£)*)a(s — <z(0)> s P eC( K0 c C_. T/ien 

9 £ ^tang Ot£. 
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Proof. By Lemma [3] we have aq = —nif^ 1 + fif$ l with integrals fj = 
/j(0 defined there. Using (!38l . we can estimate /o(£) = / /(s,£) _1 ds 
as follows: 

/oo 
5- 1 (l + \s\ 2 + \er 1 ds = 7 r/5(0, 
-oo 

(£) = (1 + l^l 2 ) 1 / 2 . The integrand /(s,0 _1 remains integrable after ap- 
plying d r , h \7, and v V finitely many times. Therefore these derivatives 
can be interchanged with the integral. In view of the symbol properties 
of /, we deduce, using estimates as above, fo G S^ng a ^ Using an 
upper bound /(s,£) < (Hfls) 2 + (£> 2 ) Id , w e derive /„(£) > 5(f) -1 Id, 
again in the sense of selfadjoint maps. Therefore /o is an elliptic sym- 
bol, and /q -1 G 5^ at £. 

Write f\ = f w + f n , where / 10 (£) = J] s |< 1 sa/(s, £) _1 ds, 

J|s|>l 

Recall s 2 a — f(s) = —s(ai + a*) — (a 2 — p). Reasoning as in the proof 
of f G S^g, we see that the integrand of /n and its derivatives are 
integrable. Moreover, we deduce fn G S^ang- It is easy to see that 
fio G S^g. Therefore, at £, /i G S^g. The lemma follows. □ 

For a /i-tempered family (w^) G h~°°C~ 0O (X) the semiclassical wave- 
front set WF h (u h ) C T*X U S*X is defined, [Ger88j . [SZ02] . Below we 
deal with operators associated to symbols which are not defined on all 
of T*X but only at £. These operators are defined microlocally in £ 
by letting them operate on the subspace of distributions (uh) which 
satisfy WF h (uh) C S, modulo the space h^C 00 . 

Lemma 14. Let q be as in Lemma [XI Microlocally at £, 

(39) h 2 L - p = (D r/h - Q")4)(A. /fc - Q), 

where Q,Q^ G ^tan g ? suc h that Q — Op h (q),Q^ — Op h (q*) G ^tang- #ere 
A is as in Proposition [ill 

Proof. Initially we set Q = Op h (q) and = Op h (q*). At £, 

(40) h 2 L-p= (D r/h - Q*)A (D r/h -Q)+R 1 + R D r/h , 

where Rj G ^tang- Here we used the formula for L given in Proposi- 
tion [TTJ Observe that, if A G ^tan g , then the commutator [D r / h , A] 
belongs to ^tang -1 - Aiming at an inductive construction, we assume 
that (|4"01) holds for some positive integer k with Rj G ^(ang ~ k - The 
spectra of q and q* are disjoint. It follows that the equation sq — q*s = r 
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has, at £, for every symbol r £ S' m a unique solution s £ S*" 1 Ap- 
plying this construction to the principal symbols of the Rj's, we find 
operators Sj £ ^tang"* sucn that SjQ — Q^Sj — i?j £ ^tang Set 

Qx = Q - Ao^SoQ + 50, Q} = Q s + (Q»S + Si)V- 

Then 

{D r/h -Q{)A Q {D r/h -Q 1 ) 

= (D r/h - Q s )A (D r/h - Q) 
+ (S Q - Q i S )D r/h + (S 1 Q - Q i S 1 ) 
+ [Dr/h, SoQ + Si] - (Q*S + S 1 )Ao 1 (S Q Q + Si). 

Replace Q and Q$ by Q\ and Q\, respectively. Then, by the symbol 
calculus, (flOl) holds with smaller errors, Rj £ ^tara' - *" 1 - The proof is 



V J c * tang 

completed using asymptotic summation. □ 

It follows from the foregoing construction that the symbol of Q is 
classical. 

6. A Dirichlet Parametrix 

Microlocally at £, we solve, constructing a parametrix, Bf = u, the 
Dirichlet problem h 2 Lu — pu = 0, u\x = f ■ We adapt the method of 
|Tay96[ 7.12] to our setting. 



Denote S™ is C C°°([— 1, 0], C°°(T*X; tt* End(CTxM))) the space of 
symbols b(s,T]), — 1 < s < 0, r\ £ T*X, which satisfy the estimates 



for all nonnegative integers r, j, and I. Let S^ is denote the correspond 



pois 



ing space of /i-dependent symbols b h . Observe that g(s(r])) £ S® ojs if 



g(t) = |tpe e *, £ > 0, j a nonnegative integer. 

We continue to work in a collar I x X G M. Choose a cutoff Xo 



as in (164T) . Let 6 > 0. Given 6^ £ S 1 ™', introduce the operator B 



pois 



h 



Op 5h (b h (r/h)) as follows: 

B h f(r,y) = (2nh)- n / / r —/-<<-,,, ; :,/< ^ r , 

(41) 




■ 6^(r//i, y, rj)T^l x J£ vV] f(ex.p y v) dv drj, 

r & I, n — dimX. We call Bh a Poisson operator with symbol bh and 
(exponential) decay 5. The arguments in |Tay96 Ch. 7 Prop. 12.4] ap- 



ply to give B h : L 2 (X) -> H~ m+1/2 (IxX) with norm 0(h- k+l / 2 ). (The 
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Sobolev spaces are defined using hD instead of D.) If < 8' < S, 
j E N, then r j B h E Op a> S^ ,k ~ s . Moreover, B h f E C°° in r < 0, 
and Bhf(r) decays together with its derivatives as e s ' r ^ h , uniformly if 
/ ranges in a bounded subset of L 2 (X). We call /i-dependent oper- 
ators negligible if they have Schwartz kernels which are smooth and 
O c °°(MxX)(h°°). We write A = B iff A - B is negligible. Note that 
Bh in f l4Tj) is negligible if there exists e > such that bh(s,rj) = if 
-e < s < 0. 

We need to handle the composition of a Poisson operator with a 
tangential operator. The following lemma deals with this when the 
symbols are classical, i.e., they possess asymptotic expansions in powers 
ofh. 

Lemma 15. Let < 5' < 5. Let A = Op h a(r) and B = Op Sh b(r/h), 
where a = a(r, r]) E and b(s, rf) E S™ ois are h-independent symbols. 
Then AB = Op s , h c(r/h), where c = Ch E S™£ s '° has an asymptotic 
expansion c ~ Ylj>o & c j> c j £ S™oL ~^ ■ The principal term equals 

c (s lV ) =a(0, 77)6(5, r,)e^- 5 '^. 

Proof. Using Taylor expansions, a(r,rj) = J2j<n r ^ a j(v) + r7Va Ar( r ; v)i 
and the properties of r^B noted above, we may assume without loss 
of generality that a does not depend on r. Arguing as in the proof 
of Proposition [281 we can write, at least formally, AB = Op Qh c(r/h), 
where 

c(s, x, = (2nh)-^ J e^ h a(x, ^r^^h^ y, C)e^> 

T x xT*xT x xT* 

■ M(x, w + v, v) d(v, r), w, 

ip as in (1721) . We use the standard arguments in handling compositions 
of symbols: dyadic decompositions and the method of (non-)stationary 
phase. We infer that there exist e > and dj E S™^ ~ J , d (s,r]) = 
a(r])b(s,r)), such that for every N, 

eM = ( h j dj{s,rj))e SB U> + d Nh (s, V )e es ^ , 

j<3N 

where d^h £ ^^a~ ~ ■ Observe that (0 / (v) is uniformly bounded 
from below if £ and rj range in the same dyadic shell. Above we have 
chosen e less than 5 times this bound. Define Ch(s, rf) as the product of 
an asymptotic sum 

T,j>o hjd A s ,v) with the symbol e**"*'^ E S° ois . 
It follows that AB — Op s , h c{r/h) belongs to Op e h S™J~^~ N ~ f° r every 
N. Thus AB = Op s , )h c(r/h). " □ 
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Let q and Q as in Lemma [TU If r\ G £ ranges in a set having a 
positive distance to the complement of the elliptic region, then there 
exist positive constants 8$ and M such that 

(42) | e »*9(o,»7)| < Me sSo{r >\ s<0. 

This follows from the fact that the spectrum of q(0, rj)/ (rj) is contained 
in a compact subset of the lower halfplane then. We shall solve (D r m — 
Q)B = 0, B\ r=0 = Id, microlocally at £. On the symbol level we have 
to solve linear ordinary differential equations with constant coefficient 
matrices. The following assertions are true microlocally in £ where 
(g2D holds. 

Lemma 16. Let < 5 < 5 . Let r G and v G S m . Let b(s,rj) be 

the solution of the initial value problem 

(43) d s b( Sl ri) = (iq(0,r]) - 5(r]))b(s,T]) +r(s,T]), -1< s < 0, 
and 6(0, rf) = viji). Then b(s, rf) G S 



m 

pois ' 



Proof. Note that the coefficient matrix of ff43l) does not depend on s. 
Representing b by Duhamel's formula and using fT42l) we derive the 
estimate 

r° 

\b(s,r))\ < M\v(r))\ + M / e (<5 °-' 5)s<r ' ) |r(s,?7)|ds 



< M\v(r])\ + (M/(5 - 5)) sup \r{s,r])\/(r)). 

s<0 

Moreover, we can estimate d s b(s, rj) by estimating the right-hand side 
of (j4"3l) . Differentiating fj4*3l we derive linear ordinary differential equa- 
tions for 9j( 1 'V)' 7 ( ?i V) £ &('5, r]). These equations are of the same struc- 
ture as ( 14"3"|) with the same coefficient matrix. The asserted symbol 
estimates are obtained recursively. □ 

Proposition 17. Let < 5 < 5q, and e > 0. There exists B G 
Op<5,h Spois with Schwartz kernel supported in — e < s < 0, such that, 
microlocally at £, (D r /h — Q)B = and B\ r= o = Id. Moreover, (h 2 L — 
p)B = 0. 

Proof. It follows from Lemma [T5l that, for a classical symbol b G 5*™'^, 
< 5' < 5, modulo negligible operators, the composition (D r / h — 
Q)Op Sh b(r/h) equals Op s , h c(r /h), c G S^ s 1,k . Moreover, c is classi- 
cal, and, modulo S™^ 1 , 

c(s, 77) = ( - id s b(s, 77) - i5(r})b(s, rf) - q(0, r])b(s, r]))e {S - s ' )s{r,) . 
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Fix a sequence (Sj), 5 < Sj + i < 5j. Using Lemmas [15] and [16] we 
recursively find /i-independent symbols bj G S &i| r =o = Id, bj\ r=0 = 
if j > 1, such that -Bj = ft/ -1 Op Sjh bj(r/h) satisfy 

(D r/h -Q)( Bl + ... + Bj) G Op,, h Sl;t J , <5'< 5 r 

Now B is contructed using asymptotic summation. The last assertion 
follows from the factorization ([391) . □ 

7. The Displacement-to-Traction Operator 

In this section we deal with operators on the boundary X. Therefore, 
in the following, operators and symbols £1X6 ; cLS cL rule, evaluated at 
r = 0. 

Let B denote the Dirichlet parametrix given in Proposition [T7] and 
T the traction defined in ([T]). The operator Z = hTB is called the 
semiclassical displacement-to-traction operator, or Neumann operator, 
at £. By Propositions [TT1 and [TTl we have, if WF/ l (/) C 8, 

Zf = (iA D r/h Bf + ihA x Bf)\ x = iA (0)Q(0)f + ihA^f- 

Therefore, Z = iA Q + ihAi, and Z is, microlocally in £, a pseudo- 
differential operator of class ty 1,0 . The symbol of Z is classical since 
the symbols of Aj and Q are. 

Lemma 18. The displacement-to-traction operator Z is, in S, up to 
a negligible operator, formally selfadjoint. The principal symbol of Z 
equals the surface impedance tensor 

(44) z = i(aq + a x ) G S\£\ vr* End(CT x M)). 
The leading symbol of Z is z + hz~, where z_ G S* , 

(45) z_q - q*z_ = itv(S)z + id r z - a 2 _ - a x _q + tr ( v \7 q* -a h V q) ■ 

Proof. Let f h f 2 G L 2 (X;CT X M), WF h (fj) C £, and set Uj = Bf). 
By®, 

/ {Zf\\f 2 )&V x - [ (h\Zf 2 )dV x 
J x J x 

= h^ 1 / (ui\h 2 Lu 2 — pu 2 ) — (h 2 Lui — pui\u 2 ) dV^-- 

J M 

It follows from Proposition [T7] that the right-hand side is 0(h°°), uni- 
formly if the /j's range in a bounded set and have /i-wavefronts con- 
tained in a common closed subset of £. Thus Z* = Z in £. Recalling 
Z = iAoQ + ihAi, we infer from the symbol calculus that z = i{aq + a{) 
is the principal symbol. 
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It remains to prove the formula for z_. Write the leading symbols 
of Q and as q + hq^ and q* + hq'L, respectively. It is easy to see 
that z_ = iaq-. Recall the formula for L in Proposition [TTJ The 
factorization ( 1391) is equivalent to 

(D r/h - ihti(S))hA 1 + {hA\ - ihtr(S)A )D r/h + h 2 A 2 - p 

= -D r/h A Q - Q*A D r/h + Q*A Q. 

This in turn is equivalent to the following two equations of tangential 
operators: 

hA x + hA\ - ihtr{S)A + A Q + Q*A = 0, 

[D r/h , A Q + hA^ - ihtT(S)hA! + h 2 A 2 -p- Q^AqQ = 0. 

On the principal symbol level these equations become a± + a\ + aq + 
q*a = and a 2 — p — q*aq = 0. These equations agree with (JTDj) . On 
the leading symbol level the equations become, after division by h, 

-idx- -itx(S)a + aq^ + qta - ztr ( v \7q*- h Va) = 0, 
— d r z — itr(S)ai — ia 2 - — q*aq^ — qi_aq + ztr [ v \/ q* . h \/ aq) = 0. 
Elimination of gl from these equations gives 

{aqJ)q — q*(aqj) = ia^q + tr(S)z + d r z + ia 2 _ —iti ( v \7 q* .a hj \7 q). 
Formula (H5j) for Z- = iaq_ follows. □ 

In principle z_ is found as the unique solution of the linear equation 
(|45|) . The right-hand side of the equation consists of known quantities 
and their first order derivatives. Refer to section [10] for an algorithm 
computing Z- if the elastic medium is isotropic. 

8. DlAGONALIZATION OF Z 

Assume ((HI) and (1E1I) . By Proposition [6] the kernel ker z defines a 
line bundle over the characteristic variety E = of the surface 

impedance tensor z. Since zero is a simple eigenvalue of z at E, there 
exist e > and an open neighbourhood K C £ of E such that 2(£), £ G 
has exactly one eigenvalue A (£) of modulus < e. (In the following, 
K is to be replaced by a smaller neighbourhood when necessary.) The 
line bundle E = ker(z — A ) -> if is a subbundle of n*CT x M = 
Hom(C, tt*CTxM). The orthoprojector onto this bundle is given by a 
contour integral, u = (27ri)~" 1 ^j A | =e (A — z)^ 1 dA. Denote u\ = Id — u 
the orthoprojector onto the orthogonal bundle, E\. 

Assume also (1E2j) . Choose a unit section v of kerz E, |u| = 1. 
Using uq, extend v to a unit section of E -> K. Call this section also 
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v. Clearly, Uq = v <E> v*. If R G \1/ ' denotes the inverse of a square 
root of the scalar operator Op h (v)* Op h (v), then V = Op h (v)R satisfies 
V*V = Id, i.e., V is an isometry. 

Lemma 19. Choose V G ^f°'°(K; C, CTxM), with principal symbol v, 
such that V*V = Id. Set U = VV* , U x = ld-U . There exist B G 
^" 1 '- 2 (X; CT X M), B* = B, and R G tf" 1 >~\X\ CT X M), R* + R = 0, 
such that, microlocally in K, 

(46) (Id - R*)Z (Id -R) = U (Z + B)U + U X (Z + B)U X . 
In particular, 

(47) (Id -R*)Z (Id -R)V = VV*(Z + B)V. 

The leading symbol of the scalar operator V*(Z + B)V G iff 1 ' equals 

(48) A + h(z-v\v) - thtr (v* v Vz. h Vv + v Vv* . h V X -v) . 

Here, as in Lemma lT^ z + hz- denotes the leading symbol of Z . 

Proof. To prove (j4"6l) we adopt ideas of jSteOOj . The operators Uq and 
XJ\ are orthogonal projectors, U* = Uj = Uj, and U\Uq = 0. Write Z = 
UqZUq + U\ZUi + B, where B = U ZUi + U\ZU§. Since UjZ = zuj and 
u x Uo = we have B G \I/ '~ 1 . Let hb, b = b* G S°, denote the principal 
symbol of B. Define the section Zj = z\e 4 of End(Ej). The spectra of z 
and Z\ are disjoint. Therefore the Sylvester equation sz — z x s = Uibu 
has a unique solution s which is a section of Hom(£'o, E x ). We extend 
s to a section of End(7r*CTx-M) by s = u\SUq. Then sz — zs = u x buo, 
and s G S' 1 . Define S = Op h (hs) and R = UqS*^ - U^Uq. Then, 
R* = -R and B = U Q BU X + U x BUq = R*Z + ZR modulo if- 1 '- 2 . 
Therefore, with a different B G <ff~ N+1 ~ N , N = 2, and Z = Zi = Z, 
we have 

(49) (Id -R*)Z(Id -R) = U Z U + U 1 Z 1 U 1 + B. 

If iV > 2 then, using the same construction as before, we find -Ri G 
q-N,-N^ R * = _ Ri guch that UqBUi + c/jSC/q = + modulo 
1 j r -A',-Jv-i_ H ence we g e ^ (|49p with R and Zj replaced by R + -Ri and 
^ + B, respectively. The new error B belongs to \fr- N - N - 1 , Iterating 
this construction and using asymptotic summation (|46|) follows. Since 
U V = V, gSD implies (H7jl . 

Observe that the leading symbols of V*(Z + .B)V and V*ZV are 
equal. The principal symbol equals (v \zv) = Ao because \v\ = 1. We 
write the leading symbol of V as (l+Zj/y^+Ztw;, where v*w = (w\v) = 0. 
Note (v\zw) = 0. A straightforward symbol computation, using (1681 
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and ( I7TT) . gives 

a h (V*ZV) = X + h(z^v\v) + /i( 7 + 7)A 

-ihtx (v* v Vz. h Vv+ v V h Vv*.zv+ v Vv*. h Vzv) 

modulo 0(h 2 ). From V*V = 1 it follows that the leading symbol of 
V*V equals unity. Since \v | 2 = 1 is the principal symbol, this implies 

Kl + 7) - ihtT ("V ft Vw*.w + w V«*. h Vu) = 0. 
Therefore the expression for the symbol of V*ZV simplifies to 
a h (V*ZV) = \ + h(z_v\v) 

+ ihtr (\ v Vv*. h Vv - v* v Vz. h Vv - v Vv*. h Vzv) 
modulo 0(h 2 ). Using h \J zv = \ h Vv+ ft VA .v we deduce ( 148ft . □ 

Denote ^phg the class of /i-independent pseudo-differential operators 
A with polyhomogeneous symbols, a ~ Sj<m a i> a j homogeneous of 
degree j. When regarded as an /i-dependent operator, A G \l/ m ' m has 
the classical symbol Ylj<m^ l ~^ a j- m ^ ne nex t lemma, following |PV99j 
and [SteOO] . we use this relation to conjugate the scalar operator con- 
structed in Lemma [T9l into hP — 1, where P is h- independent. 

Recall that f2 1//2 -> X denotes the bundle of half-densities. 

Lemma 20. There is a selfadjoint operator P G \l/ ph (X; fi 1 / 2 ) wift 
principal symbol p, and an operator A G \l/ ' /rom half-density sections 
to scalar functions, elliptic nearH, such that A*V*(Z+B)VA = hP—1 
in a neighbourhood o/E. T7ie subprincipal symbol of P equals, on E ; 

Psub = (if|i>) _1 ( Re(,2_t;|t;) + Imtrff^V^-^Vv)) 
(50) \ 

+ Imtr h Vp-"Vf*-f- 

Here z denotes the radial derivative of z. If instead of v another unit 
section v = e tLp v of ker z -> E is used to define V , and thus P, then 
the principal symbol of P remains unchanged, whereas the subprincipal 
changes to p su b = p su b + {p, V 9 } onT,. Here {p, <p} denotes the Poisson 
bracket. 

Obviously, P is elliptic and bounded from below. 

Proof. The radial derivatives of p and of A = (zv\v) are, at E, equal 
to 1 and (zv\v) > 0, respectively. Therefore, near E, Oq\ = p — 
1 for some a G C°°, a > 0. Set Z = A*V*{Z + B)VA , A = 
Op h (ao). Choose Pi G \E f1,0 (formally) selfadjoint with leading symbol 
p — z/itr^V ft Vp)/2. The selfadjoint operators Z and Pi — 1 have the 
same principal symbol, p — 1. Therefore, the imaginary parts of their 
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leading symbols are equal. It follows that the principal symbol go of 
Q = Z — (Pi — 1) G equals, on E, = (zv\v)~ l times the real 
part of the coefficient of h in (HB1) . 

Define p G C°°(T*X \ 0), homogeneous of degree 0, and r_i G S^ 1 
such that go — Po + 2(p — l)r_i holds in a neighbourhood of E. Then 

Op ft (r_i)*)£(l - /i P/l (r_!)) = A + hP - 1 + 

where Po is selfadjoint with principal symbol po- Proceeding induc- 
tively, we obtain selfadjoint operators Pj G \I/ JI ' with classical symbols 
such that, for N < 1, 

(1 - hR* N )Z(l - hR N ) = h h ~ j Pj ~ 1 + h ~ N QN, 

N<j<l 

where Qn G ty N '°, R N g \l/ _1,0 . Therefore, there is an /^-independent 
operator P G ^ hg such that (1 - hR*)Z(l - hR) = hP - 1 near E. 
Moreover, P = Pi + hP modulo ^/~ 1, ~ 2 . The symbol of P equals 
p — itr^V h Vp)/2 +po modulo It follows from Corollary l29l or 
rather its analogue for /i-independent operators, that p is the principal 
symbol of P and p su b = Po its subprincipal symbol. By construction 
Po — go on E. Formula ( )50l) follows from the formula for go mentioned 
earlier. 

Note {p, if} = tr ( v Vp- h V<f- h Vp- v Vf)- The last assertion of the 
lemma follows from fl50|) . using v* v \7 z.v = ^V^o- □ 

Proof of TheoremUi The following assertions hold microlocally in a 
neighbourhood of E. It follows from lemmas [19] and [20] that, if A~* 
denotes a parametrix of A*, (Id -R*)Z(ld -R)VA = VA~*(hP - 1). 
Define J h = (ld-R)VA and J h = {Id-R^VA^*. We have J h , J h G 
* 0,0 , J h - Jh e Moreover, JJJJ h is elliptic. By definition of Z, 

TBhJh — Jh{P — h^ 1 ), where R>h is the Dirichlet parametrix given in 
Proposition \T7\ Combining the results in section [6] with lemmas [19] and 
1201 the theorem follows. □ 

9. Construction of Quasimodes 

Given P of Theorem [1] we associate to the sequence of positive eigen- 
values of P a sequence of quasimodes of Lt- We follow |Ste00l sect. 4], 
differing in some details, however. 

Let P, Bh, and Jh as in Theorem [1] Assume given a sequence of 
quasimodes, (fij), with almost orthogonal quasimodes states: 

(51) Pfj - fijfj = O c ~{hf), (fj\f k ) - S jk = 0((hj + h k )°°), 
fj G L7°°(X; n 1 / 2 ), 0<Hj< ->• oo, h, = fij 1 . 
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We define quasimode states for the traction-free boundary problem. 
By Theorem [T] the traction tj = TB h .J hj fj = Oc°°{h°°)- Choose u'j = 
Oc^iti**) satisfying A (0)d r u'j\x + tj = and u'j\x = 0. Define Uj G 
C?(M;CTM), 

(52) u j = hJ 1/2 (B h .J h .f j + u' j ). 
By Theorem (TJ 

(53) Luj - fjtjpuj = O c ™(h°°), T Uj = 0, 

and ||wj||l 2 — We can assume that the Uj are supported in 

a given neighbourhood of X. Using the ellipticity of L, we deduce 



-2\ 



J 



\\uj\\ff2 = 0(h 

To go from quasimodes to eigenvalues or, in scattering theory, to 
resonances, it is desirable to be able to decompose the quasimodes 
into well-separated clusters. In addition, the quasimode states of each 
cluster should be linearly independent, and remain so after applying 
small perturbations. 

Proposition 21. Let the assumptions of Theorem [1\ hold. Assume 
given quasimodes fij = hj 1 > of P as in (I5TI) . and define Uj as in 
fl52|) . Then fl53|) holds. Let m > dimX. There exist 5 > and a 
covering of {fij} by a sequence of intervals [ak, 6fc] C M +; such that 

b k + 25b- k m '^ x < a k+1 , h-a k < b~ k m . 

Let wj G H 2 (M; CTM) such that, for some N > 0, 

\\ Wj \\ H2 = 0(hJ 2 ~ N ), Wj - Uj = L ,{h)^ x+N ). 

Then, for large k, {wj} ak <fij<b k is linearly independent. 

Proof. Property (|53|) is clear by the arguments already given. 

It is well-known that a quasimode sequence (151]) is asymptotic to a 
subsequence of the sequence of eigenvalues of P. The latter satisfies 
the Weyl asymptotics. Hence we have a Weyl estimate j < Cfij imX . 
It follows that every interval [a, b], 1 < b, of length > L has a subin- 
terval of length > Lb~ AimX jC which contains no quasimode fij. The 
existence of intervals [a^, bk] having the stated properties follows from 
this observation. Compare [St e99l Proof of Theorem 2] . Define the set 
of indices of the k-th cluster: Ik = {j ; \ij G [ak, bk]}. 

Choose a left inverse K h G ^°'°(X; CT X M, ft 1 / 2 ) of J h , K h J h = Id at 
E. Since J^Jh is elliptic at E, Kh is readily found. By the remark after 
Lemma [26] we may assume that h i— > Kh is Lipschitz into the space of 
bounded operators on L 2 . 



SURFACE QUASIMODES 



31 



Denote 7 : v i-> u |x the trace map. By (1521) . ft,- 7^ = J/^/j + ju'j. 
From f l5T]) it follows that WF^ fj C E. Therefore, 

h) ,2 K hj ju,j = fi + Oc^hf). 
If C denotes a Lipschitz constant of ft 1— >• i^, then, for j, I G 7^, 

||i^ - i^/ l J| i 2_ i . i 2 < C|ft^ - ftj| < Ca^ 2 ]/^ - /ijl < Ca k 2 b k m . 
Therefore, 

ftflK^ - K h .) 1Uj \\ L 2 = 0(b~ m ~ 2 ) } j,ie h, 

if k is sufficiently large. The assumptions on Wj imply \\wj — Uj\\jp. = 
0(h 1 j +dimX ). Here we use the estimate \\v\\ 2 H1 < C\\v \\ L 2 \\v\\ H a. Apply- 
ing the trace theorem, ||7iUj — jUjW^ = 0(h^ +dimX ). Summarizing the 
estimates, we have shown that, for some e > 0, 

\\h] /2 K helWj - fjWv = 0(h £ e +dimX ), jje I k . 

Because of almost orthogonality of the fj and the Weyl estimate, we 
can apply |Ste99[ Lemma 4]. We obtain, for every £ G I k , the linear 
independence of {Kh e jWj}j e i h when k is sufficiently large. Since K^pj 
is linear, also {wj}j e i k is linearly independent. □ 

Proof of Corollary^ We apply Proposition [211 with \ij f 00 the se- 
quence of positive eigenvalues of P, counted with multiplicities, and 
{fj} a corresponding orthonormal system of eigenvectors. Fix m > 
dimX. Let [a*,, bk] be the intervals, clustering {fij}, given in the propo- 
sition. The quasimode states defined in fl52|) belong to the domain of 
the selfadjoint operator Lt- Let iik denote the spectral projector for Lt 
of the interval [a' k , b' k ], where a! k = a k -5b k m - dimX , b' k = b k +5b k m ~ diraX . 
The intervals [a' k , b' k ] are pairwise disjoint. Set Wj = n k Uj if \ij G [a k , b k ]. 
A well-known argument, using the spectral theorem, gives 

6 2 b k 2m - 2d ™ x \\wj - Uj\\ 2 L , < \\(L T - fifaWh = <D(b k °°) 

if Hj G [a k ,b k ]. Since L T is elliptic, we have ||iUj||#2 = 0(fi 2 ). Now 
Proposition I2TI with N = 0, implies that, for k sufficiently large, the 
rank of n k equals Jj{j ; fa G [afc,6fc]}. Hence an increase by n of Np 
over [a k ,b k ] leads to an increase > n of Nl t over [a^,^]. Taking into 
account the widths of the intervals, the corollary follows. □ 

Remark. The foregoing arguments also apply to give lower bounds for 
the counting function of resonances. In this case, ir k is the projector 
onto the space of resonant states which correspond to resonances in 
rectangles [ak, b k }+i[0, s k }. To satisfy the assumptions in Proposition!!?!! 
for Wj = ir k Uj, one establishes resolvent estimates. See [SV96j . [TZ98J, 
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[Ste99] . and |SteOO] . for ways from quasimodes to resonances. The 
clustering method was developed in this context, |Ste99j . to handle 
multiplicities appropriately. Resolvent estimates for anisotropic elastic 
systems are given in jKNOOj. 



10. The Isotropic Subprincipal Symbol 

In this section we assume that the elastic medium is isotropic. We 
evaluate the subprincipal symbol of P, p su h, starting from the general 
formula fl50|) . 

We continue with Example [H referring to the notation introduced 
there. The kernel bundle ker z is a line subbundle of V, the subbundle 
of CT X M spanned by u, | = |. Abbreviate (p3D and (USD as follows: 

Cl -<2 
<2 C3 



(iq) 



11 



IK 21 



-IK 12 

K-22 



It will be convenient to use the slownesses relative to the Rayleigh 
wave speed, a s = c r /c s and a p = c r /c p . Then t = a 2 , "•' 



a 2 ) 1 / 2 t 

u s) 1 'p 



a 



2X1/2 



S = { c r|CI = !}■ Moreover, we set t s = (1 

We first show how to evaluate (z_v\v), v G kerz, z_ as in (H5l) . 

Lemma 22. Set K = (iq) n . Define Yj by ([55]), (J55J, and fl5S]). Let 

X = (xjk) the selfadjoint 2x2 matrix which is the unique solution of 

(54) XK + K*X = -2Yi -Y 2 - Y* + Y 3 

Let v = V\v + v 2 £ G kerz. Then 



v* 

I 3 ■ 



2 Re(z_f \v ) = xn|fi| 2 + x 22 \v 2 \ 2 + 2 Rex^Wi^a- 



Proof. Set x = Then 2Ke(z^v\v) = (xv\v). By ( H5i) . x satisfies 

the uniquely solvable Sylvester equation x{iq) + {iq)*x = iy + (iy)*, 
where y equals the right-hand side of (jl5"|) . Since q leaves V and V 1 - 
invariant, X = (x)n = (x^) is the unique solution of (1541) provided 
the right-hand side of the equation equals (iy + (iy)*) n - The latter 
holds if 

Y 1 = (tr(S)« + ^) U) F 2 = (a x _ig) F 3 = (ztr "V(i?)*.a fe Vzg) u . 



Observe that the a 2 _ term of (1451) drops out because of the skewness 
of (ia2-)n- In the following we derive formulas for Yj. 

The basis vectors v and ^ do not depend on r. Therefore, (d r z)u = 
d r (z)n. We obtain 



(55) 



Fx = tr(S) 



Ci 

<2 



<2 



+ 



drCl 

id r ( 2 



-id r ( 2 
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Using Lemma [12] and the remark following it we obtain a formula 
for (ai-)n. Clearly, (a 1 _iq) u = (ai_)n(ig)n. We derive 



(56) 



Y, 



(t VA) n tr S + ( A + n) (1 SO \ [ i«2i ^22 
It remains to determine I3. Fix an orthonormal frame (77^-) of T%M, 
r/i — u, r]2 — i- To compute the contraction we use the frame (T)j)j>2 
of T*X, and the dual frame. We compute derivatives of 

iq = \£\Vl - t(Id -1/ ® 1/ - 1 <g> 

Set Sjk = (Sr}j,r}k). A calculation using h \7 v = S and ^V^ = gives 
( h V j iq) u = h Vj(iq)u + s^b^M, j > 2, where 

(ut-b)y/T=t 
i{ut — t) 



M 



21 



{ut-b)y/T=t 

Regard the coefficients Kj k as functions of c s ,c p , |£|. Then h Vj(iq) 11 = 
(rjj, Vc s )K s + (r/j, Vcp)K p , where K s and denote the partial deriva- 
tives of (ig)n with respect to c s and c p , respectively. In particular, 

Define »i = [(^ — 6)^/1 - £, — i(b— ut)). The row fc > 2 in ( h Vj 
equals Sj^b^l^Wi. 

The vertical derivative of a function k which, when restricted to a 
fiber depends only on |£|, is given by its radial derivative: 

(57) v V n K=\£\-%ri)k. 

A calculation using v \7 u = and "VC = — £ <8> gives 

("Viigjn = "Vi(ig)n = \^%jK, j > 2, 
where we have set 

Ku —iK-13 

Define w 2 = [i(b - *), y/l-ut - VI - t\. The row k > 2 in ("Viig) 21 
equals Sji t b~ 1 W2- 

Denote A = (a) n = diag[A + 2jj,,fi\. Note that (0)22 equals \x times 
the unit matrix. Summing over j > 2 we derive 

F 3 = i K*A(\i\- 1 (i,Vc s )K s + Ifr 1 ^, Vc p )if p + 8226"^) 

+ i>6 -2 |f Ktr(S') - S22M <8> iui, 

evaluated at S. □ 
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Denote v the unique unit section of ker(^ — A ) satisfying (£\v) > 0, 

w = 7" 1 «2^+(Ci-A )e), 

7 > such that \v\ = 1. We compute the w-dependent terms in the 
right-hand side of ( l50i) . 

Lemma 23. On E, Imtr h \7p. v \7v*.v = 0, and 

I67 2 Imtrty^V^V^) 
(59) = m 3 /U - 1 c r V s 6 (4 - a s 2 )(2 - a 2 )(2r s C 3 - (2 - a 2 )^)^ 

+ 2m 3 c r a s 6 (2 - o?)(5^ - 4 - a s 4 ) tr'(S), 

where tr'(S) = tr(5) — S22, S22 = (£, anc ^ m — 

Proof. Set 71 = (2/7 and 72 = (Ci — A0V7. We continue to use the 
frame (77^- ) . For j > 2 we have 

Note that ^Vj v*.v is real. Hence Imtr h \7 p. v \7 v* .v = 0. We need the 
vertical derivative of z. To compute it we proceed in the same way as we 
did when computing the derivatives of iq. Recall that z equals Id 
on V ± , where C" 1 = A*|f |V1 = *• We obtain { v V j z) ll = v Vj(z) n . 
Moreover, the column k > 2 in ("Vj2)i2 equals <5jfc|£| _1 times the 
transpose of the row vector [—1(2, (3 — C" 1 ]- We get 

Im v* v Vj z. h Vj v = 71 Re v* v Vj z.Srjj 

= 71(72 B VjCs- 7i "ViC2)s 2 j 

+ 7l|er 1 (72(C3 - C X ) - 7lC 2 )^(l - Sy). 

Summing over j > 2 we obtain 

7 2 Imtr(^V^V^) = CiCa'Vrf Cs - C^V^G 

+ c r C 2 (Ci(C3-C ± )-C2 2 ) tr'(S). 
The first term on the right equals 

m 2 \i\- l s 22 {2b - ^(tVT^Cs - (2b - *)C a ). 
Moreover, using the definition of b, we calculate 

C 3 - C X = m(v / T^ - v 7 ! 3 *)- 
Using (|2TJ), 46 = t(4 - i), we derive ([59]). □ 

The restriction to X of the radial derivative of the eigenvalue Ao = 
(zv\v) = clq 2 (p — 1) equals A = (zv\v) = Oq 2 because p = 1 on E. 
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Lemma 24. On S 7 

7 2 A = m 3 a!?(4 - (j 2 )t s (t p /t s + c s t s /c p t p + a 2 - 2), 
where m = fi\^\/b. 

Proof. The section w = i&v + Cil; equals jv on E. Therefore, to second 
order on E, 7 2 A = (zw\w) = Ci det(z)n. Inserting (|20|) . 

7 2 A = m6tv / T 3 t(4 v / (l - t)(l - ut) - (2-t) 2 ). 

Recall p = c r |£| = a s t~ 1/2 , £ = {t = o" 2 }. The rule of de l'Hospital 
gives 

, 4:J(l-t)(l-Ut) - (2-t) 2 t n, . , 9 \ 

I™, — ,_ 1/2 \ = (Tp Ts + C S T S /C P T P + a] - 2 . 

Summarizing, the formula for 7 2 Ao = 7 2 Ao/(p — 1) follows. □ 

Inserting the formulas of the lemmas of this section into the general 
formula fl50|) for the subprincipal symbol of P we obtain a formula for 
the subprincipal symbol in the isotropic case. 

Proposition 25. Denote X = (xjk) the 2x2 matrix solving (|54|) . Se^ 
N = t s (t p / t s + c s T s /c p T p + a 2 — 2) . Lei P fre £/ie operator of LemmaWJ\ 

determined by the unit section v of ker z having positive £ component. 
The subprincipal symbol of P is given as follows. 

IQNpsnb = (cr/2fj)(x u (2 - a 2 ) 2 + 4x 22 (l - a 2 ) 2 + 4Imx 12 (2 - ct 2 )t s ) 

+ / i" 1 c 2 (2 - a 2 ) (2t s ( 3 - (2 - a 2 )( 2 ) (S& |) 

+ 2c r (4 - a 2 y\2 - a 2 )(5a 2 - 4 - o*)(ix(S) - (SH)) 

Proof. OnE,ti) = 7D = (mi/2) (z(2-t)z/+2 v / l - *Q • Using lemmas[22] 
and [241 we calculate 16iVRe(z_if |u>)/ 7 2 Ao. The result is the first term 
on the right-hand side of the claimed formula. Similarly, we obtain the 
other terms combining the lemmas [23] and [2H □ 

The constituents of the above formula for p su b are curvature and 
velocities (Lame parameters), assumed known. It seems difficult to 
analyze the formula further unless it is specialized to particular cases. 
However, it should be noted that the formula allows explicit numerical 
evaluation of p su b- Therefore it can be used when solving transport 
equations for Rayleigh wave amplitudes numerically with a (seismic) 
ray tracing program, say. Formulas for the amplitudes of Rayleigh 
waves were given by Babich and Kirpichnikova in [BK04] . 
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Appendix A. Geometric Pseudo-differential Calculus 

Pseudo-differential operators on manifolds are usually introduced by 
reducing to the euclidean case via partitions of unity, |Hor85j 18.1], 
[EZ] . The principal symbol of a pseudo-differential operator is invari- 
antly defined. If the operator acts on sections of the line bundle of half- 
densities then there also is an invariantly defined subprincipal symbol, 
[Hor85l Theorem 18.1.33], [SZ021 Appendix]. 

In the body of the paper we explicitly track, down to the subprin- 
cipal level, symbols of operators acting between vector bundles. To 
achieve this we use Sharafutdinov's geometric pseudo-differential cal- 
culus |Sha05a| ISha05bj . The purpose of this appendix is to recall this 
calculus, presenting a semiclassical variant. Since we have to refer, in 
the main part of the present paper, to proofs of the calculus, we give 
a rather detailed presentation. The calculus depends on a symmetric 
connection of the manifold and on metric connections of the (Hermit- 
ian) bundles. We make the stronger assumption that the manifold is 
Riemannian and that the symmetric connection is the Levi-Civita con- 
nection. The important features of the calculus are a symbol isomor- 
phism modulo order minus infinity, and complete symbol expansions 
for products and adjoints given solely in terms of geometric data. Us- 
ing connections to develop a pseudo-differential calculus and to prove 
the existence of a complete symbol isomorphism was done earlier by 
Widom, [Wid80j. This was further developed by Pffaum who gave a 
convenient quantization map from symbols to operators, [Pfl98j . Shara- 
futdinov gave symbol expansions in terms of geometric data. 

Let X a compact Riemannian manifold without boundary, dim X = 
n. The exponential map, exp, of the Levi-Civita connection defines 
a diffeomorphism, (x,v) i— > (x,y) = (x,exp x v), between a neighbour- 
hood of the zero-section of the tangent bundle T = TX and a neigh- 
bourhood of the diagonal in X 2 . In the proofs of the propositions 
below we need the following properties of exp. In local coordinates the 
exponential map satisfies 

(60) (exp x . vY = x i + v i - r jk (x)v j v k /2 + 0(\v\ 3 ), 

where T l - k denote the Christoffel symbols. Normal coordinates centered 
at x satisfy (exp x t>) 1 = v l . There exist < r < R < inj(x), the 
injectivity radius of X, such that the equation 

(61) exp exp ^ z = exp x w 

defines, for every v G T x = T X X, \v\ < R, a diffeomorphism w i— > 
z = z(x,v,w) from an open neighbourhood of the origin, contained 
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in {|iu| < R} C T x , onto the ball {\z\ < r} C T y , y = exp x v. 
This map is used below to change variables of integration. Obviously, 
z(x,0,w) = w. A computation in normal coordinates centered at x 
shows that 

(62) (z' w )- 1 z = w-v + 0{{\v\ + \w\) i ) asv,w^0. 

Recall, from section El the notation for segments and for parallel trans- 
port maps. In local coordinates, 

(63) {T^^wy = w*- V) k {x)wh k + 0(\v\ 2 ). 

Let E -y X and F -* X be Hermitian vector bundles with metric 
connections. Recall from section [3] the definition ( T2"4"l) of horizontal 
derivatives and the definition of vertical derivatives. A C°° section 
a of the bundle n* Hom(E, F) -> T*X is called a Hom(£', F)-valued 
symbol of order m G R, a G S m = S m (T*X; vr* Hom(£, F)), iff for all 
nonnegative integers j and £, 

sn P (l + \^\y~ m \( v Vy( h vYa(x,0\<oc. 

These are the usual type 1, symbol estimates. The symbol space 
S m is a Frechet space. The space S m ' k = S m ' k (T*X; n* Hom(E, F)) 
of /i-dependent symbols of order m and degree k is the Frechet space 
of families G S m such that {h k ah ; < h < 1} is bounded in 
S m . We call a G S' ni ' fc classical if there exists an asymptotic expansion 
a ~ ^ . h?~ k dj with /^-independent symbols a 3 - G S m ~K 

In the following lemma we define, in a semiclassical setting, the quan- 
tization of symbols according to Sharafutdinov's geometric pseudo- 
differential calculus. We relate this definition of /i-pseudo-differential 
operators to the definition in the euclidean situation. For semiclassical 
analysis, in particular, for the class \l/ m ' fc = Op ft S m > of /i-pseudo- 
differential operators, including mapping properties, and for frequency 
sets (/i-wavefront sets), refer to |Ger88j . [Tv^8] . [DS99] . [SZ02] . [EZ] . 
The class of negligible operators, vl/ -00 - -00 , consists of /i-dependent op- 
erators whose Schwartz kernels are C°° with 0(h°°) seminorms. 

Fix xo £ C^°(TX), real-valued, |u| < r on the support of xo(%,v), 
such that Xo = 1 m a neighbourhood of the zero-section in TX. 

Lemma 26. Let ah G S m ' k be a Hom(E, F) -valued symbol. Then 



(64) 




A h u h {x) = (2irh)- n / / e-^'^xo^v) 



T* TV 



a h (x, r])T [x ^ exPxV] u h (exp x v) dv dr], 
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defines an h-pseudo- differential operator Ah G \& m,fc (X; E, F). Given 
a point x there exists a geodesic ball U centered at x, and a symbol 
a\ G S m,k such that, for Uh compactly supported in U , 



where y = exp^f and y' = exp x w. Moreover, at x, = modulo 
gm-2,k-2^ Every h-pseudo-differential operator is, modulo negligible 
operators, of the form f )64|) . 

The measures in (1641) are the normalized Lebesgue measures of the 
euclidean spaces T x and T*. 

Proof. We shall drop the subscript h from the notation. Fix x £ X. Let 
U denote a geodesic ball with center x and radius < R. In the following 
we assume that the support of u is a compact subset of U. In (IMj) we 
replace the variables x,v,rj by y,z,(- Next we change variables in the 
integral Au(y) such that the domain of integration does not depend on 
y. Set y = exp x v. Define z = z(x, v, w) by flBTl) . Using the symplectic 
map (w,$) i — y (z, C), C — t (. z ' w )~ 1 $i we S e t 



ip = (C,z) = (z'^^z). Since z = if and only if v = w, we have 
ip{v,w,-&) = (ip(v,w)$,w — v). Here if) = ld+0(\v\ 2 + |w| 2 ) by ( 162]) . 
Decreasing the radius of U and making the linear change of variables 

9 = i[)(v, w)$, we get 



Ji(v,w) = 1 + 0(M 2 + \w\ 2 ). It follows that A restricted to U is a 
/i-pseudo-differential operator of class fy m > h . As it stands the symbol 
depends on v, 9, w. Using the standard symbol reduction procedure we 
obtain a u (exp x v, 9). Moreover, the asymptotic expansion implies that, 
at v = 0, a u - a G S m ~ 2 ' k - 2 . 

Note that Au(y) = if the distance between y and suppw is > r. 
Using a partition of unity, we infer that the class of operators given by 




T* T a 




where the kernel K is given by 
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( |64|) equals the class of /i-pseudo-differential operators with Schwartz 
kernels supported in small neighbourhoods of the diagonal. □ 

Standard arguments show that up to a negligible operator A h = 
Op h (ah) does not depend on the choice of the cutoff Xo- The space 

V m > k {X; E, F) = Op h S m ' k + ^-°°'- 00 . 

is the space /i-pseudo-differential order m and degree k. We denote the 
geometric symbol ah(Ah) = ah- 

Remark. An operator = Op h (ah) G \l/ ' is L 2 bounded, uniformly 
in h. We also need that h y A^ is a Lipschitz map from ]0, 1] into the 
space of bounded operators on L 2 if the following assumption holds: 
h i — y (X}i is a map into S° with bounded derivative. In fact, then 
A hl - A ho = Op h (b h ) dh, where b h e S°'°, b h (x, rj) = d h b h (x,r)) + 
(?7, v \7bh(x,r])). To see this, change variables, in flMI) . from r] to £ = 
rj/h. The assumption holds if is classical and given as a Borel sum. 

In the following, we often suppress from writing the /i-dependence 
of symbols, operators and distributions. Moreover, when dealing with 
integrals like (16^1) . we move, without explicitly writing this, the x- 
dependency from the domain of integration into the integrand using 
arguments as in the proof of the lemma. 

Lebesgue measure dv on T X X and Riemannian volume are related 
by / f(y)dV x (y) = J f(exp x v)J (x,v)dv, y = exp x v, with Jacobian 
J = 1 + 0(\v\ 2 ) at v = 0. Let A = A h as in (Ep. The Schwartz kernel 
K A of A, 

Au(x) = / K A (x,y)u(y)dV x (y), K A (x,y) G Rom(E y ,F x ). 
Jx 

equals in a neighbourhood of the diagonal a partial Fourier transform 
of the symbol, 

(65) K A (x,y) = (27rh)- n J e^* a *- l M h a(x,rj) d V tP(x,y)r^ y] . 

Here ip(x, y) = Xo( x , v )l Jo(x, v), y = exp x v. The symbol a is recovered 
via the inverse Fourier transform: 

(66) a(x,t) = J e^/ h (xoJo)(x,v)K A (x,ex Vx v)r [ l Px ^ x] dv 

modulo S 1-00 ' -00 . The correspondence between an operator A = Op h (a) 
and its full symbol a, named the geometric symbol of A, defines the 
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complete symbol isomorphism 

^ m ' fc (X; E, F) Z^" 00 '- 00 S m ' k (T*X- Rom(n*E, n*F))/ S^ 00 ^ 00 . 

The geometric symbol can also be computed by applying the operator 
to suitable testing functions as follows. 

(67) a(x,Z)a = M e ^^ lv)/h ^^ l v)rS^]»)\v^- 

Here A y means that A acts on functions of the variable y. In particular, 
in case E = C, the geometric symbol is obtained at the center of normal 
coordinates x J when A is applied to ^ x3 ' h and evaluated at x- 7 =0. 

We derive symbol properties and expansions using the method of 
stationary phase: 

(det(iJ/27ri/i)) 1/2 J e l ^ x))/h a{x) dx 

= exp (2' 1 ih(H~ 1 d, d)) (e lp(x)/h a(x)) \ x=Q 

j<3N 

if if G C°°, real- valued, (p'(x) = iff x = 0, H = ip"(0) non-singular, 
95(0) = 0. The remainder p(x) = tp{x) — (Hx,x)/2 vanishes to third 
order at x = 0. The expansion has the advantage, when compared to 
that obtained using the Morse lemma, of giving an efficient algorithm 
for computing the asymptotic series. See |Hor83I Theorem 7.7.5] where 
the expansion is arranged in powers of u~ l = h. 

We are mainly interested in the leading symbols of operators. We 
call the residue of a in S m,k / S m ~ 2 ' k ~ 2 the leading symbol of an operator 
Op h (a) G ^/ m ' k . The principal symbol is, of course, the residue in 

gm,k / jgm—lfk—l 

Proposition 27. Let A = Op h (a) as in (IMl) with geometric symbol 
a G S rn ' h . The formal adjoint A* G ty m,k (X; F, E) has the geometric 
symbol 

(68) b = a*-zhti v V h Va* mod 5 m ^ 2 ' fc - 2 . 

If a is classical then so is b. 

Notice that v \7 h V a* is a section of vr*(Hom(F, E) g> T ® T*). The 
trace is taken of the T ®T* part. 

Proof. The formal adjoint of A is defined by 

/ { Ul (x)\Au 2 (x)) F dV x (x)= [ (A*uM\u2{y)) E dV x {y). 

J X J X 
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The Schwartz kernel satisfies K A *(x,y) = K A (y,x)*. Recall that par- 
allel transport preserves inner products. It follows from (j6"5|) that 

K A .(x,y) = (2nh)- n [ e^^^r^aiy, r,)* d V ^(y,x), 



T* 
y 



and K A *(x,y) = if the distance between x and y is > r. Set y = 
exp x v. Define z G T y by exp y z = x. After a linear change variables 
from r] G T* to ( = t (exp' x (v))r] G T* we have 

K A .(x,y) = (2nh)- n J e^^T^afo, 77)* d( ^(y,x)/J 1 (x, v), 

T* 

X 

with Jacobian J\(x,v) = 1 + 0(|u| 2 ). Define 

b(x,0= [ e^/ h ( Xo J )(x,v)K A *(x,y)T^ x] dv. 



T x 



Inserting K A * we have 



(69) 
where 



b{x,£) = (2nh)- n J J e ilp/h aJ d(dv, 

T x T* 



<p = (Z,v) + ( V ,z) = -(£ -£,v) + ({,*), 

a = r^ y] a(y,r])*T^ x] = r^ F ' E) a(y , 77)* , 

J = xo(x, v)J (x, v)ip(y, x)/ Ji(x, v) = 1 + 0(\v\ 2 ), 

and $ = $(x,f) = exp^.(f) _1 z +v. A computation in normal coordi- 
nates centered at x shows that $ = 0(\v | 3 ) as v — y 0. If = then 
z = 0, hence i> = 0. It follows that the critical points of (p are defined 
I'.v V 0, C, (. 

Apply the method of stationary phase to fl69|) and deduce 6 G S m ' k . 
Moreover, the following asymptotic expansion holds: 

(70) 6~E^<- fl cAy(e 4 ^>/ fc a)|^. 

i 

Differentiation of the exponential factor brings out a non-zero factor 
only if it comsumes at least three derivatives with respect to v and 
at most one derivative with respect to (. It follows that the sum is 
asymptotic. Moreover, b is determined modulo S m ~ 2,k ~ 2 by the terms 
in the asymptotic sum with j < 2, b = a* — ih{d^, d v )a. Observe that 

r [Lex P ^] exp^(w) = ld Tx +Oi\v\ 2 ) as v -> 0. 
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It follows that d v a\ = h \7a*(x,()- Hence b = a* - ihti v \7 h \7a*. 
The Schwartz kernels of Op h (b) and A* are equal in a neighbourhood 
of the diagonal. Therefore A* - B G □ 

Proposition 28. Let A G ^ niA ' kA (X; F,G) and B G V TlB ' kB (X; E, F) 
with geometric symbols a and b, respectively. Set k = ftA + ks, ui = 
rriA + ms- Then AB G ^ m,k (X; E, G) with geometric symbol 

(71) c = ab-ihtr( v \7a. h Vb) 

modulo S m ~ 2 ' k ~ 2 . If a and b are classical then so is c. 

Again the trace is taken of the T £g> T* part, and the dot terminates 
differentiated expression. 

Proof. Setting y = exp x v , C = AB is given by 



Cu{x) = (2irh)- 2n //// e- l({v ' v)+ ^' z))/h a{x,r]) 




T x xT*xT y xT* 

■ T^ y] (b(y, C)T§^ cxPy z] u(exp y z)) dz d( dv dry. 

Here and in the following we do not write the cutoff factors. Let 
z = z(x,v,w) be the solution of exp y z = exp x w. The symplectic 
change of variables (w, #) >->• (z, (), ( — t (z' w )~ 1 {}, preserves the volume 
form. We get Cu(x) = J T K c (x,exp x w)u(exp x w)J (x,w)dw, with 
Schwartz kernel 

K c (x, exp x w)J (x, w) 

c = a(x,T])r^^ E ' F h(y,C)M(x,w,v). Here M(x,w,v) G GL(E X ) 
denotes the parallel transport in E along the geodesic triangle x — > 
exp x w — > exp x v — > x. It follows that the symbol of C equals 

(72) c(x, = (2nhy 2n j e^ h co d(v, V , w, ■&), 

T x xT*xT x xT* 

(p = (£, w) — (77, v) — (C, z). We introduce w — v as a new variable, w. 
Then ([72]) holds with 



<p = -(ri-Z,v) -{•&-£,&) + (>&,$), 

>m(. 

-v) 



c = a(x, v)rZ°™] E ' F) b(y, ()M(x, w + v,v) 
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Here $ = w — (z' w (x, v , w+v))~ 1 z(x, v, w+v) G T*. By ( 1621) . $ vanishes 
to third order at i> = w = 0. Clearly, v = = z at a critical point of 
</?. It follows that v = w = and r] = $ = !; define the critical points. 

Now apply the method of stationary phase to (!72|) and deduce that 
c G S"" 1 ^ is a symbol which, moreover, has an asymptotic expansion 



(73) c ~ ^L^ft,^) + (^,^» i ( e *(*.»>/fc co ) 



f;=u>=0,?7=i9=£ ' 



Using that $ does not depend on 77 and and vanishes to third order 
at t> = w = 0, we infer that the summands with j > 1 belong to 

5 m-2,fc-2 It follows that 

ah — ih{d v a, d v bM) — iha(d$, d w )bM, 

evaluated at the critical point, is the leading symbol of C. Here b = 
r [x^-™i "{ViO- We have d w b — at v — w — 0. This follows from 
C = which is a corollary of 2 = w at f = 0. The derivatives of M 
with respect to v and w vanish at v = w = 0. Using T^^p „i ^4 = 
Id^ +C(|v| 2 ) at w = 0, we derive 

&t v = w = 0. Summarizing the computations, (I7T1) follows. □ 

Remark. The proofs of propositions [27] and [28] follow those in |Sha05a[ 
ISha05bj closely with only minor modifications. Our derivation of the 
asymptotic expansions of the symbols of adjoints and products may 
be somewhat shorter, however. We differ in defining the adjoint with 
respect to the volume element rather than by using half-densities. No- 
tice that the symbol expansions ([TO]) and (173]) depend only on the given 
symbols and on the geometry. In the formulas (1681) and (ITT]) , we ex- 
tracted the leading symbols. 

For the purposes of the present paper it suffices to assume X com- 
pact. A symbol calculus on general (complete) Riemannian manifolds 
needs to take the injectivity radius into account and handle mapping 
properties more explicitly. 

It is well-known that a pseudo-differential operator acting on half- 
densities has an invariantly defined subprincipal symbol; see [SZ02, 
Appendix] for a proof in the semiclassical case. We relate the subprin- 
cipal symbol to the leading geometric symbol. Equip the half-density 
bundle f2 1//2 -> X with the inner product (u\v) = u ■ vj dVx, where 
the operations on the right are in the sense of densities. The connec- 
tion given by V fil/2 dV^ 2 = is metric with respect to the Hermitian 
structure of Vt 1 / 2 . 
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Corollary 29. Let A G ^ m ' k (X; f2 1/2 ). The leading symbol of A equals 
that of the corresponding scalar operator A G \l/ m ' fc (X) which is given 
by Au = dV^ 1 ^ 2 A(u dV^J 2 ). If the geometric symbol a of A is classical, 
a ~ Sj>o ^~ h aj, cij G S m ~i , then h~ k a$ is the principal symbol of A, 
and 

a su b = h l ~ k (ai + i v V a . h V a /2) 
is its subprincipal symbol. 

Proof. Consider the multiplication operator dV^J 2 G \I/ ' (X; C, Vt 1 / 2 ). 

The Hom(C, fi 1 / 2 )-valued symbol 7r* dV^J 2 is the leading symbol of this 
operator. Note that its horizontal and vertical derivatives vanish. The 
equality of the leading symbols of A and A now follows from Proposi- 
tion EH 

Let a u denote the local symbol of A in a geodesic coordinate chart 
U centered at a given point x. We use normal coordinates centered at 
x. Assume a classical, h k a = a + ha\ + 0(h 2 ). Then a u is classical, 
and h k a u = Oq + 0(h). Moreover, it follows from Lemma [26] that 
h k a u = a + hd\ + 0{h 2 ) at x. The subprincipal symbol equals, by 
definition, h 1 ~ k (ai + 2~ 1 i^- cPao/dXjd£j). The horizontal derivative 
in the j-th coordinate direction equals, at x, the partial derivative with 
respect to xj. The formula for the subprincipal symbol follows. □ 
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RAYLEIGH-TYPE SURFACE QUASIMODES 
IN GENERAL LINEAR ELASTICITY 



SONKE HANSEN 

Abstract. Rayleigh-type surface waves correspond to the charac- 
teristic variety, in the elliptic boundary region, of the displacement- 
to-traction map. In this paper, surface quasimodes are constructed 
for the reduced elastic wave equation, anisotropic in general, with 
traction-free boundary. Assuming a global variant of a condition 
of Barnett and Lothe, the construction is reduced to an eigenvalue 
problem for a selfadjoint scalar first order pseudo-differential oper- 
ator on the boundary The principal and the subprincipal symbol 
of this operator are computed. The formula for the subprincipal 
symbol seems to be new even in the isotropic case. 



1. Introduction 

Rayleigh [R] discovered the existence of surface waves which propa- 
gate along a traction-free flat boundary of an isotropic elastic body and 
which decay exponentially into the interior. The propagation speed of 
the surface wave is strictly less than that of body waves. Barnett and 
Lothe [LB 76] showed that Rayleigh-type surface waves can also exist 
at flat boundaries of anisotropic elastic media. 

The goal of this paper is to construct, for elastic media which are 
not necessarily isotropic, Rayleigh-type surface quasimodes which are 
asymptotic to eigenvalues or resonances. We use a geometric version 
of semiclassical microlocal analysis. 

The Rayleigh wave phemomenon of isotropic elastodynamics was 
explained by Taylor |Tay79| as propagation of singularities, over the 
elliptic boundary region, for the Neumann (displacement-to-traction) 
operator. Nakamura [Nak91j generalized this to anisotropic media, 
using the theory of Barnett and Lothe. Assuming isotropy of the elas- 
tic medium, Cardoso-Popov |CP92j and Stefanov [SteOOj constructed 
Rayleigh quasimodes. 
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Let (M, g) be an oriented Riemannian manifold with non-empty com- 
pact smooth boundary X. The (infinitesimal) displacement of an elas- 
tic medium occupying M is a vector field u on M. The Lie derivative of 
the metric tensor is a symmetric tensor field, Def u = C u g/2, called the 
deformation (strain) tensor caused by the displacement u. The elastic 
properties are defined by the elasticity (stiffness) tensor. This is a real 
fourth order tensor field C G C°°(M; End(T°' 2 M)), e \-> Ce, which 
maps into symmetric tensors and vanishes on antisymmetric tensors. 
We assume positive definiteness of C, i.e., (e\f)c = (Ce\f) defines an 
inner product on the space of symmetric tensors e and /. Here (-|-) 
denotes the inner product on tensors induced from g. This assumption 
is often called the strong convexity condition. If coordinates are given, 
then the components of C satisfy symmetries, C 1 ^ = & = C i , 
and eijCkt > if is a nonzero symmetric tensor. (We use the 
summation convention.) Denote the Riemannian volume elements on 
M and X by dV^ and dVx, respectively. The elasticity operator L 
and the traction T are defined, on compactly supported vector fields, 
by 

(1) / (Def u\ Defy) c dV M = [ (Lu\v)dV M + [ (Tu\v)dV x . 
Jm Jm Jx 

A positive mass density p G C°°(M) and the elasticity tensor C de- 
fine the material properties of the elastic medium. If the surface X is 
traction-free, then vibrations of the medium are solutions of the follow- 
ing eigenvalue problem: Lu = \ 2 pu in M, Tu = at X. See [MH83J 
for linear elasticity in the language of Riemannian geometry. 

The principal symbol of L, and of the /i-differential operator h 2 L, 
equals the acoustic tensor, c(£) = c(£, £) G End(CT x M), £ G T*M; 
see fl3Tl) . Here the associated acoustical tensor c(£, rj) G End(CT x M), 
£, rj G T*M, is defined as follows: 

(2) (c(£,t))v\w) — (y (8) r}\w <g) 

(Using g, we identify vectors with covectors.) The ik-th covariant com- 
ponent of c(£, 77) equals C^^rji. 

The existence of Rayleigh waves depends on the characteristic va- 
riety, S, of the surface impedance tensor, z. To define z, we first 
recall the definition of the elliptic boundary region, £ C T*X. Let v 
denote the unit exterior conormal field of the boundary X. Identify 
T*X = u ± C T* X M. By definition, £ G S iff c(f + sp) - p Id is positive 
definite for real s. From the factorization theory of selfadjoint matrix 
polynomials one gets g(£) G End(CT x M), £ G £ n T*X, such that 

(3) c(t + sv)-pld= (sld-q*(0)c(u)(sld-q(0), 
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s G C. Moreover, the spectrum of q(£) lies in the lower half-plane, 
specg(£) C C_, and these properties determine q(£) uniquely. The 
surface impedance tensor z is defined as follows: 

(4) z{Z) = ic(v)q(Z)+ic(v,Z), £e£. 

The significance of z results from the fact, proved in Lemma [181 that z 
is the principal symbol of a parametrix of the displacement-to-traction 
operator. In physics, the meaning of z is that it relates the ampli- 
tudes of displacements to the amplitudes of tractions (forces) needed 
to sustain these. 

The surface impedance tensor is Hermitian, and positive definite for 
large |£|, [BL85| Theorem 6]. If dimM = 3, then 

(U) z(£), £ e £■> nas a t m ost one non-positive eigenvalue. 

This property expresses the uniqueness of Rayleigh-type surface waves, 
[BL85j Theorem 8]. In case dimX ^ 3, we shall assume (jUJ) as a 
hypothesis. The characteristic variety of z, 

E = {£g£ ;detz(£) = 0}, 

is a smooth hypersurface, transversal to the radial directions of the 
fibers of T*X. Compare |BL85(, Theorem 7]. Rayleigh waves exist 
only if E is not empty. We shall make the stronger assumption that E 
intersects every radial line: 

(El) EnR + ^0 if £ G T*X \ 0. 

Compare [BL851 Theorem 12], [Nak91l Theorem 2.2], jKNOOl (ERW)]. 
Assuming flU} and flEll . there exists a unique p G C°°(T*X\0), p > 0, 
homogeneous of degree 1, such that 

(5) S=p- 1 (l). 

See Proposition [6j Furthermore, the kernel of z defines a line bundle, 
kerz -> S, over the compact base E. We shall require that its first 
Chern class vanishes: 

(E2) ker z ->■ E is a trivial line bundle. 

In particular, the bundle is assumed to possess a unit section. Property 
(IE2p is stable with respect to homotopies in the material properties; 
see Corollary [7J In the case of isotropic elasticity with positive Lame 
parameters, ([U|) . flElj) . and flE2j) hold. Moreover, 

E = {c r |£| = 1} c £ = {c s \£\ > 1}, p(£) = c r \C\. 

Here c r is the propagation speed of the Rayleigh surface wave which is 
strictly less than the speeds of the body waves, < c r < c s < c p . See 
Example [HJ 
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Next we state the central result of this paper: The traction-free 
surface eigenvalue problem can be intertwined with a selfadjoint eigen- 
value problem on the boundary. We employ a semiclassical pseudo- 
differential calculus, with distributions and operators depending on a 
small parameter, < h < 1. We write Ah = Bh iff the Schwartz kernel 
of Ah — Bh belongs to C°° with seminorms satisfying Oc°°{h°°). 

Theorem 1. Assume dimM = 3 ; or ([U]). Assume (1E1I) . (1E2I) . Given 
a unit section v o/ker z -> £, there exists a selfadjoint, elliptic operator 
P E \E fl (X; fi 1 / 2 ), independent of h, and operators, 

B h : L 2 (X; CT X M) -» L 2 (M; CTM), \\B h \\ = 0(h 1/2 ), 

J h ,J h E fi 1/2 , CT X M), J£j h e/Mzc at £, 

(/i 2 L -p)B h = 0, TS ft J ft = J fe (P - /T 1 ), 

and -B^lx = H in o neighbourhood of E. T/ie principal symbol of P 
equals p of ([S]). Furthermore, there is a formula, (1501) . /or t/ie subprin- 
cipal symbol p sn h of P. If v is changed to another unit section, e l(p v, 
then the subprincipal symbol changes to p su b + {p, <^} ; where {p, (p} de- 
notes the Poisson bracket. 

This result is known in the isotropic case, |CP92j . [SteOOj . except for 
the assertions about the subprincipal symbol. 

The operator Bh is a parametrix of the Dirichlet problem near S; 
see Proposition [T71 Its range consists of functions which are smooth 
in the interior of M, supported in a preassigned neighbourhood of the 
boundary, and which decay like e ~ 5dmt x/h j n ^ interior. 

Ignoring finitely many eigenvalues the spectrum of P consists of a 
sequence of positive eigenvalues pj t oo. Applying Theorem [1] to an 
associated orthonormal system of eigenvectors we obtain, in Proposi- 
tion [211 a sequence of quasimode states: Luj — p 2 puj = Oc°°{hf) with 
boundary tractions equal to zero. Moreover, the quasimode states are 
well-separated. The construction also works when starting with a se- 
quence of almost orthogonal quasimode states of P. 

The unbounded operator D L 2 (M; CTM; p dV M ), u H> p~ l Lu, 
with domain D = {u 6 ; Tu = 0}, is symmetric and nonnegative. 
The associated quadratic form is given by the left-hand side of (jTJ). 
Denote Lt the Friedrichs extension of this operator. For a selfadjoint 
operator A with spectrum consisting of a sequence of eigenvalues ac- 
cumulating at +oo, denote Na(X) the usual counting function for the 
eigenvalues of A. The following lower bound on Nl t (X) is an example 
application of our results. 
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Corollary 2. Assume M compact, dimM = 3, and flEll) . ( 1E2I) . Lei 

P be the selfadjoint operator given in Theorem^ For every m > 1, 
Nl t {^) — Np(X — A~ m ) is bounded from below. 

Rayleigh waves have been studied in several papers with the empha- 
sis of getting information about resonances in scattering theory, [S V95j . 
[S V96j . |SjV97| , |Ste00j . and, for anisotropic media, |KN00] . Stefanov 
[SteOOj uses Rayleigh quasimodes to derive lower bounds on the num- 
ber of resonances. See the remark at the end of section [9] about going 
from quasimodes to resonances. 

The subprincipal symbol p su b affects the eigenvalue asymptotics of 
P, |DG75j . and it enters quasimode constructions, |CP92j . The sub- 
principal symbol occurs in the final formulas via integrals, such as 
f s , x Psub and J 7 Psub) where 7 is a closed bicharacteristic. We point 
out that these integrals do not depend on the choice of the unit section 
v in Theorem [H although p su t, itself does. An important aim of the 
present work is to give explicit formulas for the subprincipal symbol 
of P. These seem to be new even in the isotropic case which is dealt 
with in more detail in Proposition [25J The main difficulty comes from 
the fact that an invariant notion of subprincipal symbol has only been 
available for scalar operators. To overcome this obstacle we adapt and 
systematically use the geometric pseudo-differential calculus of Shara- 
futdinov [Sh a05al ISha05bj which assumes given a differential geometric 
structure. The principal and subprincipal symbol levels are contained 
in the leading symbol of a (pseudo-) differential operator. 

The paper is organized as follows. In section [2] the surface impedance 
tensor is studied; in particular, a selfcontained treatment of Barnett- 
Lothe theory is given. The leading geometric symbols of some differen- 
tial operators are computed in section [31 In section H] we geometrically 
decompose the elasticity operator near the boundary into normal and 
tangential operators, keeping track of leading geometric symbols. Sec- 
tion [5] gives, microlocally at the elliptic region S, a factorization of 
h 2 L — p into a product of first order operators. Using the factorization, 
we construct in section O a parametrix for the Dirichlet problem mi- 
crolocally at £ . The displacement-to-traction operator Z is defined in 
section and its leading geometric symbol is determined. In section [S] 
we derive a diagonalization of Z, and we prove Theorem [TJ In section M 
we construct localized traction-free surface quasimodes, and we prove 
Corollary [2j In section [10] we calculate, for an isotropic elastic medium, 
the subprincipal symbol of P. The appendix [A] contains a detailed ex- 
position of Sharafutdinov's geometric pseudo-differential calculus in a 
semiclassical setting. 
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2. The Surface Impedance Tensor 

First we collect some well-known facts about spectral factorizations 
of selfadjoint matrix polynomials. Refer to [ GLR8 2, Chapter 11]. Let 
V be a finite-dimensional complex Hilbert space, and f(s) = as 2 + bs + 
c G End(V) a quadratic polynomial in the complex variable s. The 
spectrum of / is the set of s G C such that ker/(s) ^ 0. Assume 
that the leading coefficient of /, a, is nonsingular. Then the spectrum 
is finite. Assume that / is selfadjoint, f(s)* = f(s), and that, in 
addition f(s) is positive definite for real s. The spectrum of / is a 
disjoint union <j + U cr_, where a + and cr_ are contained in the upper 
and lower half-planes, respectively. There is a unique q G End(V) such 
that f(s) — (s — q*)a(s — q), and the spectrum of q equals <r_. If 7 is 
a closed Jordan curve which contains cr_ in its interior and a + in its 
exterior, then 

(6) ql f(s)- 1 ds= (fsfisy'ds. 

The integral on the left is nonsingular. Jordan-Keldysh chains are a 
means to compute q. In particular, one has qv = sv if f(s)v = and 
Ims < 0. Moreover, the solvency equation f(q) = holds. 

The following representation of the factor q by integrals is important. 
We shall also apply it later to establish symbol properties. Denote 
i = a/— 1 the imaginary unit. 

Lemma 3. Let f and q be as above. Then 

(7) aqfo = -7rzld+/i, 

where fo = /(s) -1 ds is selfadjoint and positive definite, and 

h= f saf(s)- 1 ds+ [ s-^a-ZOO)/^)" 1 ^. 

J\s\<l J\s\>l 

The integrals converge absolutely in End(V). 

Proof. Let 7^ denote the negatively oriented closed contour composed 
of the semicircle {\s\ = R,lms < 0} and the interval [— R, R]. The 
integral representation ([6]) holds with 7 = 7^ if i? is sufficiently large. 
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We have /(s)" 1 = s~ 2 a~ l + 0(|s|~ 3 ) as \s\ — > oo. It follows that 



Urn j> f(s)' 1 ds = J™ f{s)- x ds, 



and 



lim d> saf(s) 1 ds = — 7ri Id + lim / saf(s) 1 ds. 
R ~*°° Jir R ^°° J-R 

Using s 2 af(s)~ 1 — Id = (s 2 a — f(s))f(s)~ 1 we obtain 

/ saf(s)- 1 ds= [ s -i( s 2 a _ /(s) ) /(s) -i ds . 

Jl<\s\<R Jl<\s\<R 

This proves the formulas. The remaining assertions follow from these 
and the positive defmiteness of f(s). □ 

Let £ G T*X, and denote v G T*M the unit exterior normal. Set 
a = c(v), cti(0 = c(z/, £), and a 2 (£) = c(£). Note that ai(£)* = c^;^)- 
The polynomial 

(8) /(s) = c(£ + sz/) - p = as 2 + (ai + a*)s + a 2 - p, 

/(s) = /(s,£), has values in End(CT.j.M). It is selfadjoint with real 
coefficients. By definition, £ G £ iff /(s) is positive definite for sel. 

Lemma 4. TTie elliptic region £ is an open subset ofT*X with compact 
complement. Moreover, £ is symmetric and star shaped with respect to 
infinity, i.e., G £ whenever £ G £ and t real, \t\ > 1. 

Proof. By positive defmiteness of C, there exists 5 > such that 
g(v,c(rj)v) > 5\v®7] + ri®v\ 2 for (co-)vectors ^,77. The symmetrization 
of a non-zero real elementary tensor is non-zero. Therefore, with a new 
5 > 0, in the sense of selfadjoint maps, 

(9) c{ V ) > 5\r]\ 2 Id. 

Since |£ + sz/| 2 = |£| 2 + s 2 the first assertions follow. The symmetry 
and the star-shapedness follow from c(trj) = t 2 c(r]). □ 

If £ G £, then ([7j) holds with q = g(£), = /_? (£)• The spectral 
factor g solves Q; using current notation: 

(10) as 2 + (ai + a{)s + a 2 — p = (s — q*)a(s — q). 

The spectrum of q lies in the lower halfplane, and q is uniquely deter- 
mined by these properties. Notice that q is a smooth section of the 
bundle vr* End(CT x M) -> £, where ?r : £ C T*X _> X denotes the 
canonical projection. 
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The surface impedance tensor, denned in (BJ, equals z = i(aq + ai). 
Lemma [3] implies 

(11) ^o = 7rId+z(A + aa/ ). 

Since the fj's are real, this gives the decomposition of z into real and 
imaginary parts. Following [MF04J, we shall use the Ricatti-type equa- 
tion 

(12) (z + ial)a~ 1 (z — ia\) = a 2 — p 

to deduce properties of z. Equation (|12p follows upon insertion of 
q = —a~ l {iz + ai) into the solvency equation associated with ffTU]) . 

(13) aq 2 + (ai + a\)q + a 2 - p = 0. 
A consequence of f|T2|) is 

(14) (iq)*z' + z'(iq) = <g + q*a[ + (a 2 - p)' + q*a'q, 

where the prime denotes the derivative with respect to some cho- 
sen parameter. The spectra of q and q* are disjoint. Therefore, the 
Sylvester equation {iq)*x + x{iq) = i{xq — q*x) = y has a unique so- 
lution x for given y. The solution is, in fact, given by an integral, 
x = J° exp(irq)*y exp(irq) dr. It follows that x is positive definite if 
V is- 

Proposition 5. The impedance tensor z(£), £ G £ , has the following 
properties, 
(i) z(£) is selfadjoint. 

(ii) z(£) is positive definite if |£| is sufficiently large. 
(Hi) Rez(^) is positive definite. 

(iv) z(£) has at least two positive eigenvalues if dim M > 3. 

(v) (d I ' dt)\ t= \t~ x z(t£) = z — z is positive definite. 

(vi) The complex conjugate z(£) = z(— £). 

We call £(£) = (d/dt)\ t= iz(t^) the radial derivative of z at £. It 
follows from (jvj) that i is positive definite on the kernel of z, ker z. 

Proof. To prove @ we follow the arguments in |MF04[ Theorem 2.2]. 
First note that ( |T2|) remains true if z is replaced by z* . Subtracting the 
two equations we get the Sylvester equation (iq)*(z— z*) + (z — z*){iq) = 
0, implying z — z* = 0. 

It follows from (fTTj) that Re z = 7t/q . This proves ( Imj) . 

Suppose dimM > 3. Aiming at an indirect proof of (jiv)l . assume that 
z(£), £ e h as a ^ m ost one positive eigenvalue. Then there exists 

w G CT X M such that z(£) is negative semidefinite on the orthogonal 
complement w 1 . Choose a real vector n ^ which is orthogonal to 
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both Rew and Imw. Then v G w 1 , and (Re z(£)v\v) = (z(£)v\v) < 0, 
contradicting the positive definiteness of Re z. 

Next we prove (jvj) following the method of [MF04, Theorem 2.3]. 
Since Oj(£) is homogeneous of degree j in £, equation (fr2|) implies, 

(r 1 ^) + iol(0)a- 1 (r 1 z(tO - iai(O) = a 2 (0 - t~ 2 P . 

Taking the derivative with respect to t at t = 1, we get 

— 2) + (i — = 2p. 

By the remarks following f JT4|) . we see that z — z is positive definite. 

We now prove Note /(s, -£) = f(s,£), /,•(-£) = ( -l) J 7j(fl , 
and Ol (-0 = -%(£). Using ([II]) we derive z(-£)/o(£) = *(0/o(0- 
Since /o is real and nonsingular the formula follows. 

It remains to prove (Jn|). Let 77 G T^X, |r/| = 1. It suffices to show that 
Zoo = li m t|oo t^zitrf) exists and is positive definite. Set q t = t^qitrf), 
t > 1 large. From ffTO]) deduce 

as 2 + (01(77) + a\(q))s + a 2 {rf) - t~ 2 p = (s — q* t )a(s - q t ), s G R. 

Using ([7]) and dominated convergence in the integrals giving fj we 
infer that q^ = lim^ q t exists. In particular, t~ l z(tr]) converges to 
Zoo = i(aqoo + a i(v)) as t t 00 • Let y G T X M such that (zoo2/||/) < 0. 
We must show y = 0. Set w(r) = exp(irg 00 )?/, r < 0. The solvency 
equation ( ITBl holds with g replaced by goo, p = 0. Therefore, aD 2 w + 
{ai+a\)D r w+a2W = holds, where we use the abbreviation a 3 - = 07(77). 
Take the inner product in CT X M with to and integrate. A partial 
integration gives 




(aD r w + aiw\D r w) + (D r w\aiw) + (a,2w\w) dr 



= i(aD r w + aiitj|iw)|_ = (zoayly) < 0. 

Set W(r) = w(r) <g> 77 + D r w(r) <S> v G End(CT x M). Recall a = c(u), 
ax = c(v,r)), a 2 = c(rj), and (j2J). We have shown: 

/ (W(r)\W(r)) c dr < 0. 
J —00 

Recall that C is real, and that ( | )q is an inner product on symmetric 
2-tensors. It follows that the symmetrization of W(r) vanishes for all 
r < 0. In particular, 

(15) (w(r) (g) i] + D r w(r) <g> v \ ( <g> v + v <g> C) = 

for C G CT^M, r < 0. Recall (77^) = 0. Setting ( = v, we derive 
-D r (w(r)|z/) = (D r w(r)|z/) = 0. Since w(r) — > as r — >■ —00, we obtain 
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(w(r)\u) = 0. Now, ( 1T5]) simplifies to (D r w(r)\() = 0. Since ( is 
arbitrary, this implies, successively, D r w = 0, w = 0, y = 0. □ 

If dimM = 3 then (JU|) holds. This follows from flrv]) . 

Proposition 6. Assume (jH])- T/ien t/ie characteristic variety of z, 
£ = {detz(£) = 0}, is a smooth hypersurface in £. Each radial line 
M + £ C T*X intersects £ in at most one point, and the intersection is 
transversal. The kernel of z\x defines a line bundlekerz -» S. Assume, 
in addition, (lEip . There is a unique p G C°°(T*X\0) , homogeneous of 
degree one, such that E = Moreover, p > 0, and p(— £) = 

Proof. From the assumption and ([v]) of Proposition [5] it follows that 
(d/dt) det z > if t£ G S, t > 0. In particular, zero is a regular 
value of det z. Hence £ is a codimension one submanifold transversal 
to the radial field. Since M + £ D £ is connected, a given radial line 
intersects S in at most one point. Because of (|U]) and the selfadjoint- 
ness of z, zero is simple eigenvalue of z. It follows that kerz -> E is a 
line bundle. Now assume also (lEip . Then each radial line intersects £ 
in a unique point. Define p as follows. For ^ £ G T*X set = 1/t 
if t£ G E, t > 0. Smoothness of p follows from the implicit function 
theorem. The evenness of p is a consequence of fjvTj) . The other proper- 
ties of p are obvious. Clearly, the homogeneity and p|g = 1 determine 
p uniquely. □ 

Corollary 7. Let p t and C t be homotopies of the mass densities and 
the elasticity tensors, < t < 1. Assume that the associated surface 
impedance tensors z t and their characteristic varieties £t satisfy ((HI) 
and (lElj) /or ever?/ 1. 27ie line bundles ker zo So one? ker zi -> Si are 
zsomorpazc. 

Proof. The factorization ()3]) and the definition of the impedance ten- 
sor imply that z t depends continuously on the homotopy parameter 
t. It follows from Proposition E] that the characteristic varieties are 
canonically diffeomorphic to the sphere bundle SX. We deduce that 
the Chern classes of the bundles ker z t -> SX do not depend on t. The 
assertion follows from this. □ 

Example 8. We consider, as special case, an isotropic elastic medium. 
We shall verify (JH]), f lEip . and (1E2j) . The elasticity tensor reads, in 
component notation, 

(16) C*" M = Aa V + H(9 ik 9 ji + 9 li 9 jk ), 
where A, fi denote the Lame parameters. Equivalently, 

(17) c(£,n) = Ae®n + /in®e + /io(e,")W. 
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Positive definiteness of C is equivalent to p > 0, AdimM + 2p > 0. 
We make the stronger assumption A,p > 0. Let £ G r*X. We list 
the eigenvalues s G C and the eigenvectors v G CT X M of the quadratic 
polynomial c(£ + si/) — p: 

(a) (A + 2p)(|£| 2 + s 2 ) - p = and u = £ + si/, 

(b) p(|£| 2 + s 2 ) - p = and v = s£ - 

(c) p(|£| 2 + s 2 ) — p = and v is orthogonal to £ and za 

Introduce c p = \/(X + 2p)/p and c s = \fpjp, the speeds of pressure 
and of shear waves, respectively. Assume that £ G £. This is equivalent 
to c s |£| > 1. The above eigenvalues and eigenvectors diagonalize q, 
q(£,)v = sv if Ims < 0. Denote V the subbundle of 7r*(CTxM) -> S 
spanned by v and £, and l^ -1 its orthogonal bundle. Fix the orthonormal 
frame u, £ = £/|£| of V, and choose an orthonormal frame of l/ -1 . In 
block decompositions of matrices we let the indices 1 and 2 correspond 
to V and l^- 1 , respectively. We denote (e)ij the block ij of the matrix 
which represents the endomorphism e. Observe that q leaves V and 
V 1 - invariant, (q)n — = (9)21- A simple computation gives 

lei 



(*?) 
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Uty/l — t —%{b — ut) 
i(b-t) ty/l - ut 

Here t = (c s |£|)~ 2 , u = (c s /c p ) 2 = p/(A + 2/i), 6 = 1- VI - ut^/T^t. 
Moreover, {iq)22 equals — £ times the unit matrix. The maps 

a = c(z/) and ai = c(i/, £) also leave V and l^ -1 invariant. We compute 



(19) (z)n 



ty/T^t -i (2b - 1) 
i(2b - t) ty/1 - ut 



and (2)22 = p (2(7)22- The determinant of z equals (p|£|\/l — t) dimV 
times 

(20) det(z) n = p 2 )^ 1 (V(l - t)(l - ut) - (2 - t) 2 ). 

Given k g]0, l/2[, the unique zero t g]0, 1[ is found as the solution of 
Rayleigh's cubic equation, (0 (24)], = ((t-2) 4 - 16(1 -t)(l - ut))/t. 
Define the Rayleigh wave speed c r = c s Vt G l7°°(X). Set p(f) = c r |£|. 
The characteristic variety £ equals {p(£) = 1}. Thus (jUJ) and (lElj) 
hold. Obviously, i(26- + - G kerz(£), £ G E. Observe that 

(21) 2(26 - t) = t(2 - t) on S. 
Thus 

(22) z(2-t)z/ + 2 v / l-t£Gker^, i = (c r /c s ) 2 , 

is a nowhere vanishing section of the kernel bundle. Hence also (1E2[) 
holds. This example is of course well-known. 
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Remark. The identity ( ITTT) goes back to Barnett and Lothe; compare 
[LB76[ (3.18)]. It is key to proving, in dimension three, the unique- 
ness of subsonic traction-free surface waves |BL85l Theorem 8]. The 
second assumption in Proposition [6] is needed to prove the existence 
of Rayleigh surface waves. Compare with [BL85, Theorem 12], where 
existence criteria are given in terms of the so-called limiting velocity 
which corresponds to the boundary of the elliptic region. See |Nak91t 
Theorem 2.2] for the Barnett-Lothe condition in a microlocal setting, 
and the real principal type property of the Lopatinski matrix it entails. 
See [Tan07] for a recent exposition of Barnett-Lothe theory, and for a 
treatment of isotropic and transversely isotropic media. 



3. Connections and Geometric Symbols 

The elasticity operator is defined in terms of the Levi-Civita con- 
nection and of the elasticity tensor. We use the geometric pseudo- 
differential calculus of Appendix [A] to define and compute the leading 
symbol of the elasticity operator. The leading symbol includes the prin- 
cipal and the subprincipal level. The calculus depends on the choice of 
connections. 

Equip M with the Levi-Civita connection of g. Let exp denote its 
exponential map. If x, y G M, then denote by [y <- x] the shortest 
geodesic segment from x to y, assuming its interior does not intersect 
the boundary, and that it is unique. 

Let E -y M be a (complex) vector bundle with connection V E . 
Denote G End(£ , a; , E y ) the parallel transport map along a given 
curve 7 in M from x to y, e.g., T ^ x y The connection can be recovered 
from its parallel transport maps: 

(23) Vfs(x) = -^\ t=Q T^ exPxtv] s(exp x tv). 

Denote n*E -> T*M the pullback of E to the cotangent bundle 
7r : T*M — > M. Let a be a smooth section of tt*E -> T*M. Following 
|Sha05a| [Sha05b] . we introduce the vertical and the horizontal covariant 



derivative of a. The vertical derivative v \/a(x, £) G E x ®T X M, at 
f G T*M, is the derivative of the map T*M ->■ E x , £ ^ a(x,£). The 
definition of the vertical derivative depends only on the linear structure 
of the fibers of T*M. The horizontal derivative h Va(x,£) G E X ®T*M 
is the derivative at v = of a map T X M — > E x , 

d 

(24) h Va(x, = — I _t£ ex „ v] a(exp x v, tLJ?^). 



dv 
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The horizontal derivative depends on the Riemannian structure and on 
the connection V B . In the scalar case, E — C, in local coordinates, 

h Va(x,0 = {d Xj a(x,0 +T%(x)Z k ds t a(x,Z)) dx* , 

where T^- denote the Christoffel symbols of the Levi-Civita connection. 
Writing a local section of n*E as a sum of products ai(x, £)a,2(x) where 
a\ is scalar and a 2 a section of E one readily derives local formulas 
for the horizontal derivative in terms of connection coefficients. The 
vertical and the horizontal derivative extend to first order differential 
operators, "V and ^V, which map sections of 7r*(E(g)T r ' s M) to sections 
of n*(E®T r+1 ' s M) and of 7i*(E®T r ' s+1 M), respectively. The operators 
"V and commute. It suffices to prove this when E is the trivial line 
bundle, E — C In this case the assertion is easily checked in normal 
coordinates. 

Let F -> M be another vector bundle. Let A : C°°(M; E) -> 
C°°(M; F) be a differential operator of order m. We introduce a small 
parameter, < h < 1, and we replace A by the /i-differential operator 
h m A. Then A G \I> m ' (M; E, F) as a semiclassical (pseudo-) differential 
operator. Refer to Appendix |A] for an exposition of Sharafutdinov's 
geometric pseudo-differential calculus in a semiclassical setting. The 
formula (167]) for the geometric symbol, <Jh{A) G S m '°, simplifies to 

(25) a h (A)(x,Os = M^' CXP ' ly)/ SU s )\ y = x > 

where £ e T*M, s G E x , and z = y/—l. The geometric symbol extends 
by continuity to the boundary of M. In symbol computations we track 
the leading geometric symbol, defined before Proposition [571 In the 
following, the symbol of an operator is always its geometric symbol. 

For the Laplace-Beltrami operator one has cr/ l (— h 2 A)(x, £) = |£| 2 . 
This is readily checked using normal coordinates. 

From ([25]) and flUD deduce 

(26) a h {-ihV E ){i)e = e®ieE x ®T*M. 

As before, to ease notation, we usually do not write the base point x 
into the arguments of tensors and symbols. 

If E -> M is a Hermitian vector bundle then we define, using the 
volume element <1Vm, the Hilbert space L 2 (E). Assume E and F are 
Hermitian vector bundles having metric connections. The leading sym- 
bol of the formal adjoint A* of A is given by 

(27) a h (A*) =a h (A)* -ihtT( v V h V(T h (Ay). 
See Proposition 1271 
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Equip the bundle E £g> T*M with the induced Hermitian structure 
and the induced connection. The connection is metric. Observe that 
the horizontal derivative of ah{—ihV E )* vanishes. Therefore, (1271) and 
flU imply 

(28) a h ((-ihV E y)(0(e®v)=9(Z,v)e, t,V G T^M, e e E x . 

By Proposition [28] the leading symbol of a composition is given as 
follows: 

(29) a h (AB)=a h (A)a h (B)-ihtr( v Va h (A). h Va h (B)). 

The trace is the contraction of the TM®T*M factor which is produced 
by a pair of vertical and horizontal derivatives. The dot terminates a 
differentiated expression, serving as a closing bracket. 

Let C e C°°{M;End(E ® T*M)). View C as an operator which 
acts by multiplication on sections of the bundle E ®T*M ->• M. Let 
V denote the connection on the bundle End(E £g> T*M) -> M induced 
from the Levi-Civita connection and from V s . Define sections c, dive 
of 7T* End(E) -> T*M as follows: 

c(Oe=(e,7r*C(e<8»0). 

(divc)(Oe = ^(^', (7r*V v .C)(e ® 0>, 
j 

where the angular brackets denote contractions on covectors, using g. 
Furthermore, (vj) and (rf) are any dual frames of TM and T*M. 

Lemma 9. o h (-h 2 V E * o C oV E ) = c - ihdiv c + 0(h 2 ). 

Proof. Observe that a h (C) = n*C, and h \7n*C = tt*VC. The symbol 
( 1281) is linear in £. Its vertical derivative is obvious. Using ( J29|) . the 
symbol a of —ihSJ E o C is found to be 

a(0(e®?7) = (e,7r*C(e®77))-i/iJ](77 J ',7r*V, ) .C(e®7 ? )). 

i 

Here (uy) and (r]- 7 ) are as in the definition of dive. The horizontal 
derivative of the symbol of —ihV E vanishes. Therefore, 

a h (-h 2 V E * o C o V E )(Oe = a(0a h (-zhV E )(0e = a(£)(e <g) 0, 

where we used ( 1291) . □ 

Now assume = CTM and C the elasticity tensor. Identify 

End(CT°' 2 M) = End(CTM ® CT*M). 
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Let L the elasticity operator defined in (TI]). Recall from Riemannian 
geometry the following relation between the Levi-Civita connection and 
the Lie derivative: 

(30) (£ u g)(v, w) = g(V v u, w) + g(v, V w u), 

for (real) vector fields u,v,w. Using the symmetries of the elasticity 
tensor we get 

L = Def oC o Def = (-iV)* o C o (-iV). 
We obtain the following corollary to Lemma [9j 

(31) a h (h 2 L~ p) = c - pld-ih&w c + 0(h 2 ). 

If respresent C with respect to some given local coordinates, then 
QUID reads 

a h (h 2 L - p)(0 lk = Cv%& - p5 lk - V^lh&^tt + 0{h 2 ). 

The vertical bar followed by j means covariant differentiation with 
respect to the j-th coordinate. If the elastic medium is isotropic then 
the leading symbol becomes 

*h(h 2 L - p)(0 = p{c 2 p \i\ 2 - 1)P(0 + p{c 2 s \i\ 2 - l)(Id ~P(0) 
-z/i(VA®£+(V/i®£)* + (£,V//)Id), 

where P(£) = | <8> | denotes the orthogonal projection to the propaga- 
tion direction £ = £/|£|. 

4. The Elasticity Operator in a Boundary Collar 

In a boundary collar, ] — e, 0] x X C M, we write the elasticity op- 
erator L in terms differential operators on X having coefficients which 
depend on r e /, the negative distance to X. 

Let N(x) G T X M denote the unit exterior normal at x G X. There 
exists e > such that, if we set J =] — e, 0], the exponential map of the 
Levi-Civita connection defines a diffeomorphism onto a neighbourhood 
of X in M: 

J x X — > M, (r, x) t-> y = exp(riV(x)). 

Essentially without losing generality, we assume that this map is onto 
M. The inverse map is y ^ (r,x), where — r = d(y,X) is the distance 
from y to X, and x = p(y) is the unique point in X closest to y. 
The distance function r satisfies the (eikonal) equation |Vr| = 1 in M. 
Extend X to M by X = Vr. Also introduce the unit conormal field 
v = dr. The level hypersurfaces 

M r = {yeM;r + d(y,X) = 0} 
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are diffeomorphic to X = Mo. The shape operator S = ViV is a 
field of symmetric endomorphisms of TM, g(Su,v) = g(u,Sv). The 
second fundamental forms of the level hypersurfaces M r assign (u, v) h-> 
—g(Su, v) (Weingarten equation). The dependency of the metric tensor 
on r is given by the formula {C^g){v,w) = 2g(Sv,w). This formula 
follows from (|30|) . Introduce J e C°°(I r x X), the solution of d r log J = 
trS, J\ r =o = 1. Then we have the following formula for the volume 
form of M: 

(33) / f{y)dV M (y)= [ [ f(exp(rN(x)))J(r,x)dV x (x)dr, 

J M J I J X 

f E Cp°(M). See |Pet98t Ch. 2] for the geometry of hypersurfaces using 
distance functions. 

Let E -> M be a vector bundle with connection V s . Denote E r -> M r 
the bundles induced by the inclusions M r C M, reJ. Set i?x = -So- 
Let u G C°°(M] E) be a section of i£. Using parallel transport in E 
along the geodesies which intersect the boundary orthogonally, define 
u:I^C°°(X;E x ), 

u( r )(x) = u(r,x) = T { ^ y] u(y), if y = exp(riV(x)). 

The map 

(34) C°°(M; E) — )■ C oc (J, C°°(X] E x )), u^u, 

is an isomorphism of Frechet spaces. The isomorphism commutes with 
bundle operations such as tensor products and contractions. 

The covariant derivative in normal direction is transformed into d r 
under the above isomorphism: 

(35) V|w(r) = d r u{r), r el. 

To see this, consider the geodesic / —> M, r \-t y(r) = exp(rN(x)). 
The tangent vectors are y(r) = N(y(r)). Using f l23|) . it follows that 

(V£ (v(r)) u)(y(r)) = ^\ s=r r§ ir ^y( s)] u{y{s)) 
j£ d I , ., > 

This implies (ESJ). We have V N N = SN = 0. It follows that d r N = 0, 
and d r v = 0. Abusing notation, we write d r to denote Vfr. 

Define e u e = e g) v and i v {e ® r\) = (rj, v)e. Notice that e v and i v 
commute with d r . 

Let F -> M be a another vector bundle with a connection. Let 
5 : C°°(M;E) -> C°°(M;F) be a differential operator. Assume that 
i? is tangential. This means, by definition, that i? commutes with the 
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distance function r, [B,r] = 0. Then, for every r £ I, B restricts 
to an operator B r : C°°(M r ;E r ) -»■ L7°°(M r ;F r ), B r U = (Bu)\ Mr , 
where u is a section of E -> M which extends a given section [/ of 
£7 r -> M r . The assumption [-B,r] = implies that B r is well-defined. 
Parallel transport along the geodesies orthogonal to X defines bundle 
isomorphisms E r = Ex and F r = Fx- Via these isomorphisms the B r 's 
induce differential operators B(r) : C°°(X;E X ) -> C°°(X;Fx-), called 
associated with B, such that Bu(r) = B(r)u(r), r G /. Each B(r) is 
a differential operator having coefficients which are C°° with respect 
to r. Conversely, an operator B is tangential if it is given in this way 
by a family of differential operators {B(r) ; r G /} with coefficients 
depending smoothly on r. 

Lemma 10. Let E M be a real vector bundle with connection V E . 
Then 

(36) V E = e u d r + B, 
where B is tangential. Moreover, B(0) = V Ex . 

Here E cg> T*M carries the induced connection. The lemma extends, 
by decomposition into real and imaginary parts, to complexifications of 
real bundles with connections. In particular, it holds for complexified 
tensor bundles with the Levi-Civita connection. 

Proof. Let P ± ,pW E C°°(M; End(TM)) denote the orthogonal pro- 
jectors onto the span of iV and onto its orthogonal complement, iV -1 , 
respectively. Identify E®T*M with Hom(TM, E). Let u e C°°(M; E). 
We have the following decomposition in C°°(M; Hom(TM, E))\ 

\/ E u = (V E u)P L + (V E u)P 11 = (V%u) ®v + Bu. 

This defines B, and implies fl36l) . Note that B is tangential. We have 

S(0)u(0) = Bu\ x = {{V E u)pW)\ x = (V Ex (u\ x )){P ll \x). 

This proves the asserted formula for -B(O). □ 

Assume E ->■ M a Hermitian bundle with a metric connection. Us- 
ing fl33l) . and the fact that parallel transport preserves inner products, 
we have 

(37) / ( u \v) E dV M = I I {u\v)E x JdV x dr, 

J M J I JX 

if u, v G C^°(M; E). Formal adjoints of differential operators on M 
are taken with respect to these inner products. The inner product 
of sections u and v of E x -> X is J x (u\v)e x dVx- Formal adjoints 
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of operators A(r) associated with a tangential operator A are defined 
with respect to this inner product. 

Next we prove a formula which expresses the elasticity operator L as 
a quadratic polynomial in D r = —id r with tangential coefficients. Now 
assume E = CTM, and let B as in (I36I) . Define tangential operators, 

A = L u Ce U} A 1 = -ii u CB, A 2 = B*CB. 

The order of Aj is j. Moreover, A\ = iB*Ce u . 

Proposition 11. The elasticity and traction operators defined in (CQ) 
are as follows: 

L — (D r — i tr S) (A D r + A x ) + A{D r + A 2 , 
-iT = A (0)D r + Ai(0). 
Furthermore, A*(0) = A^O)* . 

Proof. Let u,v £ C 6 °°(M; CTM). It follows from ([30]) and the symme- 
try properties of the elasticity tensor that 

f ( Def u | Def v) dV M = [ (CVu \ Vv) dV M - 

J M J M 

Inserting ( 1361) and using the definition of Aj, the right-hand side equals 

/ {i v CVu\d r v) dV M + [ {B*CVu\v)dV M 
Jm Jm 

(A d r u + iAiu\d r v) dVu + / (— iA\d r u + A 2 u\v ) dVu 
M Jm 



Partial integration with respect to r is done, using (1371) . as follows: 

/ (w\d r v)dV M = [ [ {w\d T v)J dV x dr 
Jm J i Jx 

= [ (w{0)\v(0)) dV x - [ [ ((d r \ogJ)w + d r w\v)JdV x dr. 
Jx J i Jx 

Summing up we have 

/ (Def« | Defv) c dV M 
Jm 

= / ((D r - itiS)(A D r + A x )u + A\D r u + A 2 u \ v) dV M 
Jm 

+ / (A o (0)(d r u)(0)+iA 1 (0)u(0)\v(0))dV x . 



' x 

Comparing with the formulas for L and T follow. The last assertion 
follows because J = 1 at X. □ 
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Next we compute the leading symbols of the operators (associated 
with) Ay The symbols are r-dependent sections of 7r*End(CTxM) -> 
T*X. Dropping tildes, the symbol of A equals 

a h (A ) = a = c{u) G C°°{I, C°°{T*X; tt* End(CT x M)). 

Introduce the divergence of the acoustic tensor restricted to X as fol- 
lows: 

(div x c)(C)v = (V a , {^ Va C){v ® C)>, 

if C ^ T*M, v G T X M, x G X. Here (v a ) and (j] a ) are any dual frames 
of TX and T*X. If local coordinates are chosen such that r is one 
coordinate and the other coordinates are constant along the geodesies 
orthogonal to X, then (divx c)(C)* fc = C l \ a Ct- Here the summation 
convention is used with Latin indices refering to all coordinates, and 
Greek refering to all coordinates except r. We also need the contraction 
(C,S) G C°°(M; End (TM)) of the elasticity tensor with the shape 
operator, in coordinates, 

//"> G\*^ riiikl o c 

\y,o) = u bj£, bj£ = Uj\£. 

(Because of VatS 1 = one can also write Greek coordinates instead of 
j and £.) 

Lemma 12. Let a% and ci2 denote the principal symbols of the h- 
differential operators hA\ and h 2 ' A2, respectively. At r = 0: ai(£) = 
c(z/, £), and 02(C) = c (0- On the leading symbol level, o~h{hAi) = ax, 
o~h{hA\) = a* — ihax-, and ah(h 2 A 2 ) = a 2 — iha 2 - + 0{h 2 ), where, at 
r = 0, 

ai - = (dw x c)(u) + 7r*(C,S), a 2 _(£) = (div* c)(£). 

Proof. By Lemma [TUl we have 

hAi(0) = l u C o (-iTiV), h 2 A 2 (0) = (-ihV)* 0C0 (-ihV), 

where V = V TX is the Levi-Civita connection of the boundary. We 
compute the leading symbol of hAi(0) using the composition formula 
Recall The vertical derivative of the symbol of i v C vanishes, 
Hence 

a h (hAx)(0)(0 = o- h (Mi(0))(0 = c(v,£), £ G T* X M. 

The formula for ah{h 2 A2(0)) follows from Lemma [91 In view of ff2Tj) . 
a\- = tr "V^VaJ. Since a\{£) = c(£,v) = (£, n*{Ce v )) is linear in £, 
its vertical derivative is immediate. Hence 
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Now, V v (Ce v ) equals {S7 v C)e v plus a contraction of C with V^z/ = Sv, 
proving the formula for ai_. □ 

If the elastic medium is isotropic, ( Tl6|) . then a straightforward com- 
putation shows that, at r = 0, 

(div x c)(C) = (VA <g> C) + (V/i ® 0* + (C, V/i) Id, 
(C, 5) = (A + //)£' + (/itrS) Id. 

Here VA, V/i G TX C TxM are the gradients of the Lame parameters 
restricted to X. 



5. Microlocal Factorization 

We factorize, microlocally in the elliptic region, the /i-differential 
operator h 2 L — p into a product with right factor hD r — Q, where Q is 
a tangential /i-pseudo-differential operator such that the spectrum of 
its principal symbol is contained in the lower halfplane, C_. 

As in the previous section we identify M with a boundary col- 
lar / x X, and sections of CTM -> M with r-dependent sections of 
CTxM -¥ X. Operators are polynomials in D r /h = hD r with tan- 
gential /i-(pseudo-)differential operators as coefficients. The latter are 
quantizations Bh = Op h (bh) G ^tang, of tangential symbols, 

h e SZn g = C°°{I, S m ' k {T*X; vr* End(CT x M))). 

By Proposition [Til the principal symbol f(s,£) = c(£ + sv) — p of 
h 2 L — p at £ + is a second order polynomial in s. View s as the 
symbol of D r /h- The coefficients are /i-independent tangential symbols. 
By (jnj), there exists a constant 5 > such that 

(38) /( S ,0>^(l + |s| 2 + |£| 2 )ld, s6K, 

holds if £ is sufficiently large. If F C £ is closed and R > 0, then 
F\ {|£| > -R} is compact. Hence there exist < e',5 such that (1551) 
holds uniformly for (r, £) G [— e', 0] x F. We say that a property holds 
at the elliptic region £ if it is true in every open subset of I x £ where 
fl38|) holds uniformly. 

Recall from section [2] that we have a unique spectral factorization 
([TO]) at £. 

Lemma 13. Lei q = £ G £ ; i/ie unique solution of the spectral 
factorization /(s, £) = (s — <7(£)*)a(s — <z(0)> s P eC( K0 c C_. T/ien 

9 £ ^tang Ot£. 
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Proof. By Lemma [3] we have aq = —nif^ 1 + fif$ l with integrals fj = 
/j(0 defined there. Using (!38l . we can estimate /o(£) = / /(s,£) _1 ds 
as follows: 

/oo 
5- 1 (l + \s\ 2 + \er 1 ds = 7 r/5(0, 
-oo 

(£) = (1 + l^l 2 ) 1 / 2 . The integrand /(s,0 _1 remains integrable after ap- 
plying d r , h \7, and v V finitely many times. Therefore these derivatives 
can be interchanged with the integral. In view of the symbol properties 
of /, we deduce, using estimates as above, fo G S^ng a ^ Using an 
upper bound /(s,£) < (Hfls) 2 + (£> 2 ) Id , w e derive /„(£) > 5(f) -1 Id, 
again in the sense of selfadjoint maps. Therefore /o is an elliptic sym- 
bol, and /q -1 G 5^ at £. 

Write f\ = f w + f n , where / 10 (£) = J] s |< 1 sa/(s, £) _1 ds, 

J|s|>l 

Recall s 2 a — f(s) = —s(ai + a*) — (a 2 — p). Reasoning as in the proof 
of f G S^g, we see that the integrand of /n and its derivatives are 
integrable. Moreover, we deduce fn G S^ang- It is easy to see that 
fio G S^g. Therefore, at £, /i G S^g. The lemma follows. □ 

For a /i-tempered family (w^) G h~°°C~ 0O (X) the semiclassical wave- 
front set WF h (u h ) C T*X U S*X is defined, [Ger88j . [SZ02] . Below we 
deal with operators associated to symbols which are not defined on all 
of T*X but only at £. These operators are defined microlocally in £ 
by letting them operate on the subspace of distributions (uh) which 
satisfy WF h (uh) C S, modulo the space h^C 00 . 

Lemma 14. Let q be as in Lemma [XI Microlocally at £, 

(39) h 2 L - p = (D r/h - Q")4)(A. /fc - Q), 

where Q,Q^ G ^tan g ? suc h that Q — Op h (q),Q^ — Op h (q*) G ^tang- #ere 
A is as in Proposition [ill 

Proof. Initially we set Q = Op h (q) and = Op h (q*). At £, 

(40) h 2 L-p= (D r/h - Q*)A (D r/h -Q)+R 1 + R D r/h , 

where Rj G ^tang- Here we used the formula for L given in Proposi- 
tion [TTJ Observe that, if A G ^tan g , then the commutator [D r / h , A] 
belongs to ^tang -1 - Aiming at an inductive construction, we assume 
that (|4"01) holds for some positive integer k with Rj G ^(ang ~ k - The 
spectra of q and q* are disjoint. It follows that the equation sq — q*s = r 
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has, at £, for every symbol r £ S' m a unique solution s £ S*" 1 Ap- 
plying this construction to the principal symbols of the Rj's, we find 
operators Sj £ ^tang"* sucn that SjQ — Q^Sj — i?j £ ^tang Set 

Qx = Q - Ao^SoQ + 50, Q} = Q s + (Q»S + Si)V- 

Then 

{D r/h -Q{)A Q {D r/h -Q 1 ) 

= (D r/h - Q s )A (D r/h - Q) 
+ (S Q - Q i S )D r/h + (S 1 Q - Q i S 1 ) 
+ [Dr/h, SoQ + Si] - (Q*S + S 1 )Ao 1 (S Q Q + Si). 

Replace Q and Q$ by Q\ and Q\, respectively. Then, by the symbol 
calculus, (flOl) holds with smaller errors, Rj £ ^tara' - *" 1 - The proof is 



V J c * tang 

completed using asymptotic summation. □ 

It follows from the foregoing construction that the symbol of Q is 
classical. 

6. A Dirichlet Parametrix 

Microlocally at £, we solve, constructing a parametrix, Bf = u, the 
Dirichlet problem h 2 Lu — pu = 0, u\x = f ■ We adapt the method of 
|Tay96[ 7.12] to our setting. 



Denote S™ is C C°°([— 1, 0], C°°(T*X; tt* End(CTxM))) the space of 
symbols b(s,T]), — 1 < s < 0, r\ £ T*X, which satisfy the estimates 



for all nonnegative integers r, j, and I. Let S^ is denote the correspond 



pois 



ing space of /i-dependent symbols b h . Observe that g(s(r])) £ S® ojs if 



g(t) = |tpe e *, £ > 0, j a nonnegative integer. 

We continue to work in a collar I x X G M. Choose a cutoff Xo 



as in (164T) . Let 6 > 0. Given 6^ £ S 1 ™', introduce the operator B 



pois 



h 



Op 5h (b h (r/h)) as follows: 

B h f(r,y) = (2nh)- n / / r —/-<<-,,, ; :,/< ^ r , 

(41) 




■ 6^(r//i, y, rj)T^l x J£ vV] f(ex.p y v) dv drj, 

r & I, n — dimX. We call Bh a Poisson operator with symbol bh and 
(exponential) decay 5. The arguments in |Tay96 Ch. 7 Prop. 12.4] ap- 



ply to give B h : L 2 (X) -> H~ m+1/2 (IxX) with norm 0(h- k+l / 2 ). (The 
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Sobolev spaces are defined using hD instead of D.) If < 8' < S, 
j E N, then r j B h E Op a> S^ ,k ~ s . Moreover, B h f E C°° in r < 0, 
and Bhf(r) decays together with its derivatives as e s ' r ^ h , uniformly if 
/ ranges in a bounded subset of L 2 (X). We call /i-dependent oper- 
ators negligible if they have Schwartz kernels which are smooth and 
O c °°(MxX)(h°°). We write A = B iff A - B is negligible. Note that 
Bh in f l4Tj) is negligible if there exists e > such that bh(s,rj) = if 
-e < s < 0. 

We need to handle the composition of a Poisson operator with a 
tangential operator. The following lemma deals with this when the 
symbols are classical, i.e., they possess asymptotic expansions in powers 
ofh. 

Lemma 15. Let < 5' < 5. Let A = Op h a(r) and B = Op Sh b(r/h), 
where a = a(r, r]) E and b(s, rf) E S™ ois are h-independent symbols. 
Then AB = Op s , h c(r/h), where c = Ch E S™£ s '° has an asymptotic 
expansion c ~ Ylj>o & c j> c j £ S™oL ~^ ■ The principal term equals 

c (s lV ) =a(0, 77)6(5, r,)e^- 5 '^. 

Proof. Using Taylor expansions, a(r,rj) = J2j<n r ^ a j(v) + r7Va Ar( r ; v)i 
and the properties of r^B noted above, we may assume without loss 
of generality that a does not depend on r. Arguing as in the proof 
of Proposition [281 we can write, at least formally, AB = Op Qh c(r/h), 
where 

c(s, x, = (2nh)-^ J e^ h a(x, ^r^^h^ y, C)e^> 

T x xT*xT x xT* 

■ M(x, w + v, v) d(v, r), w, 

ip as in (1721) . We use the standard arguments in handling compositions 
of symbols: dyadic decompositions and the method of (non-)stationary 
phase. We infer that there exist e > and dj E S™^ ~ J , d (s,r]) = 
a(r])b(s,r)), such that for every N, 

eM = ( h j dj{s,rj))e SB U> + d Nh (s, V )e es ^ , 

j<3N 

where d^h £ ^^a~ ~ ■ Observe that (0 / (v) is uniformly bounded 
from below if £ and rj range in the same dyadic shell. Above we have 
chosen e less than 5 times this bound. Define Ch(s, rf) as the product of 
an asymptotic sum 

T,j>o hjd A s ,v) with the symbol e**"*'^ E S° ois . 
It follows that AB — Op s , h c{r/h) belongs to Op e h S™J~^~ N ~ f° r every 
N. Thus AB = Op s , )h c(r/h). " □ 
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Let q and Q as in Lemma [TU If r\ G £ ranges in a set having a 
positive distance to the complement of the elliptic region, then there 
exist positive constants 8$ and M such that 

(42) | e »*9(o,»7)| < Me sSo{r >\ s<0. 

This follows from the fact that the spectrum of q(0, rj)/ (rj) is contained 
in a compact subset of the lower halfplane then. We shall solve (D r m — 
Q)B = 0, B\ r=0 = Id, microlocally at £. On the symbol level we have 
to solve linear ordinary differential equations with constant coefficient 
matrices. The following assertions are true microlocally in £ where 
(g2D holds. 

Lemma 16. Let < 5 < 5 . Let r G and v G S m . Let b(s,rj) be 

the solution of the initial value problem 

(43) d s b( Sl ri) = (iq(0,r]) - 5(r]))b(s,T]) +r(s,T]), -1< s < 0, 
and 6(0, rf) = viji). Then b(s, rf) G S 



m 

pois ' 



Proof. Note that the coefficient matrix of ff43l) does not depend on s. 
Representing b by Duhamel's formula and using fT42l) we derive the 
estimate 

r° 

\b(s,r))\ < M\v(r))\ + M / e (<5 °-' 5)s<r ' ) |r(s,?7)|ds 



< M\v(r])\ + (M/(5 - 5)) sup \r{s,r])\/(r)). 

s<0 

Moreover, we can estimate d s b(s, rj) by estimating the right-hand side 
of (j4"3l) . Differentiating fj4*3l we derive linear ordinary differential equa- 
tions for 9j( 1 'V)' 7 ( ?i V) £ &('5, r]). These equations are of the same struc- 
ture as ( 14"3"|) with the same coefficient matrix. The asserted symbol 
estimates are obtained recursively. □ 

Proposition 17. Let < 5 < 5q, and e > 0. There exists B G 
Op<5,h Spois with Schwartz kernel supported in — e < s < 0, such that, 
microlocally at £, (D r /h — Q)B = and B\ r= o = Id. Moreover, (h 2 L — 
p)B = 0. 

Proof. It follows from Lemma [T5l that, for a classical symbol b G 5*™'^, 
< 5' < 5, modulo negligible operators, the composition (D r / h — 
Q)Op Sh b(r/h) equals Op s , h c(r /h), c G S^ s 1,k . Moreover, c is classi- 
cal, and, modulo S™^ 1 , 

c(s, 77) = ( - id s b(s, 77) - i5(r})b(s, rf) - q(0, r])b(s, r]))e {S - s ' )s{r,) . 
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Fix a sequence (Sj), 5 < Sj + i < 5j. Using Lemmas [15] and [16] we 
recursively find /i-independent symbols bj G S &i| r =o = Id, bj\ r=0 = 
if j > 1, such that -Bj = ft/ -1 Op Sjh bj(r/h) satisfy 

(D r/h -Q)( Bl + ... + Bj) G Op,, h Sl;t J , <5'< 5 r 

Now B is contructed using asymptotic summation. The last assertion 
follows from the factorization ([391) . □ 

7. The Displacement-to-Traction Operator 

In this section we deal with operators on the boundary X. Therefore, 
in the following, operators and symbols £1X6 ; cLS cL rule, evaluated at 
r = 0. 

Let B denote the Dirichlet parametrix given in Proposition [T7] and 
T the traction defined in ([T]). The operator Z = hTB is called the 
semiclassical displacement-to-traction operator, or Neumann operator, 
at £. By Propositions [TT1 and [TTl we have, if WF/ l (/) C 8, 

Zf = (iA D r/h Bf + ihA x Bf)\ x = iA (0)Q(0)f + ihA^f- 

Therefore, Z = iA Q + ihAi, and Z is, microlocally in £, a pseudo- 
differential operator of class ty 1,0 . The symbol of Z is classical since 
the symbols of Aj and Q are. 

Lemma 18. The displacement-to-traction operator Z is, in S, up to 
a negligible operator, formally selfadjoint. The principal symbol of Z 
equals the surface impedance tensor 

(44) z = i(aq + a x ) G S\£\ vr* End(CT x M)). 
The leading symbol of Z is z + hz~, where z_ G S* , 

(45) z_q - q*z_ = itv(S)z + id r z - a 2 _ - a x _q + tr ( v \7 q* -a h V q) ■ 

Proof. Let f h f 2 G L 2 (X;CT X M), WF h (fj) C £, and set Uj = Bf). 
By®, 

/ {Zf\\f 2 )&V x - [ (h\Zf 2 )dV x 
J x J x 

= h^ 1 / (ui\h 2 Lu 2 — pu 2 ) — (h 2 Lui — pui\u 2 ) dV^-- 

J M 

It follows from Proposition [T7] that the right-hand side is 0(h°°), uni- 
formly if the /j's range in a bounded set and have /i-wavefronts con- 
tained in a common closed subset of £. Thus Z* = Z in £. Recalling 
Z = iAoQ + ihAi, we infer from the symbol calculus that z = i{aq + a{) 
is the principal symbol. 
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It remains to prove the formula for z_. Write the leading symbols 
of Q and as q + hq^ and q* + hq'L, respectively. It is easy to see 
that z_ = iaq-. Recall the formula for L in Proposition [TTJ The 
factorization ( 1391) is equivalent to 

(D r/h - ihti(S))hA 1 + {hA\ - ihtr(S)A )D r/h + h 2 A 2 - p 

= -D r/h A Q - Q*A D r/h + Q*A Q. 

This in turn is equivalent to the following two equations of tangential 
operators: 

hA x + hA\ - ihtr{S)A + A Q + Q*A = 0, 

[D r/h , A Q + hA^ - ihtT(S)hA! + h 2 A 2 -p- Q^AqQ = 0. 

On the principal symbol level these equations become a± + a\ + aq + 
q*a = and a 2 — p — q*aq = 0. These equations agree with (JTDj) . On 
the leading symbol level the equations become, after division by h, 

-idx- -itx(S)a + aq^ + qta - ztr ( v \7q*- h Va) = 0, 
— d r z — itr(S)ai — ia 2 - — q*aq^ — qi_aq + ztr [ v \/ q* . h \/ aq) = 0. 
Elimination of gl from these equations gives 

{aqJ)q — q*(aqj) = ia^q + tr(S)z + d r z + ia 2 _ —iti ( v \7 q* .a hj \7 q). 
Formula (H5j) for Z- = iaq_ follows. □ 

In principle z_ is found as the unique solution of the linear equation 
(|45|) . The right-hand side of the equation consists of known quantities 
and their first order derivatives. Refer to section [10] for an algorithm 
computing Z- if the elastic medium is isotropic. 

8. DlAGONALIZATION OF Z 

Assume ((HI) and (1E1I) . By Proposition [6] the kernel ker z defines a 
line bundle over the characteristic variety E = of the surface 

impedance tensor z. Since zero is a simple eigenvalue of z at E, there 
exist e > and an open neighbourhood K C £ of E such that 2(£), £ G 
has exactly one eigenvalue A (£) of modulus < e. (In the following, 
K is to be replaced by a smaller neighbourhood when necessary.) The 
line bundle E = ker(z — A ) -> if is a subbundle of n*CT x M = 
Hom(C, tt*CTxM). The orthoprojector onto this bundle is given by a 
contour integral, u = (27ri)~" 1 ^j A | =e (A — z)^ 1 dA. Denote u\ = Id — u 
the orthoprojector onto the orthogonal bundle, E\. 

Assume also (1E2j) . Choose a unit section v of kerz E, |u| = 1. 
Using uq, extend v to a unit section of E -> K. Call this section also 
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v. Clearly, Uq = v <E> v*. If R G \1/ ' denotes the inverse of a square 
root of the scalar operator Op h (v)* Op h (v), then V = Op h (v)R satisfies 
V*V = Id, i.e., V is an isometry. 

Lemma 19. Choose V G ^f°'°(K; C, CTxM), with principal symbol v, 
such that V*V = Id. Set U = VV* , U x = ld-U . There exist B G 
^" 1 '- 2 (X; CT X M), B* = B, and R G tf" 1 >~\X\ CT X M), R* + R = 0, 
such that, microlocally in K, 

(46) (Id - R*)Z (Id -R) = U (Z + B)U + U X (Z + B)U X . 
In particular, 

(47) (Id -R*)Z (Id -R)V = VV*(Z + B)V. 

The leading symbol of the scalar operator V*(Z + B)V G iff 1 ' equals 

(48) A + h(z-v\v) - thtr (v* v Vz. h Vv + v Vv* . h V X -v) . 

Here, as in Lemma lT^ z + hz- denotes the leading symbol of Z . 

Proof. To prove (j4"6l) we adopt ideas of jSteOOj . The operators Uq and 
XJ\ are orthogonal projectors, U* = Uj = Uj, and U\Uq = 0. Write Z = 
UqZUq + U\ZUi + B, where B = U ZUi + U\ZU§. Since UjZ = zuj and 
u x Uo = we have B G \I/ '~ 1 . Let hb, b = b* G S°, denote the principal 
symbol of B. Define the section Zj = z\e 4 of End(Ej). The spectra of z 
and Z\ are disjoint. Therefore the Sylvester equation sz — z x s = Uibu 
has a unique solution s which is a section of Hom(£'o, E x ). We extend 
s to a section of End(7r*CTx-M) by s = u\SUq. Then sz — zs = u x buo, 
and s G S' 1 . Define S = Op h (hs) and R = UqS*^ - U^Uq. Then, 
R* = -R and B = U Q BU X + U x BUq = R*Z + ZR modulo if- 1 '- 2 . 
Therefore, with a different B G <ff~ N+1 ~ N , N = 2, and Z = Zi = Z, 
we have 

(49) (Id -R*)Z(Id -R) = U Z U + U 1 Z 1 U 1 + B. 

If iV > 2 then, using the same construction as before, we find -Ri G 
q-N,-N^ R * = _ Ri guch that UqBUi + c/jSC/q = + modulo 
1 j r -A',-Jv-i_ H ence we g e ^ (|49p with R and Zj replaced by R + -Ri and 
^ + B, respectively. The new error B belongs to \fr- N - N - 1 , Iterating 
this construction and using asymptotic summation (|46|) follows. Since 
U V = V, gSD implies (H7jl . 

Observe that the leading symbols of V*(Z + .B)V and V*ZV are 
equal. The principal symbol equals (v \zv) = Ao because \v\ = 1. We 
write the leading symbol of V as (l+Zj/y^+Ztw;, where v*w = (w\v) = 0. 
Note (v\zw) = 0. A straightforward symbol computation, using (1681 
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and ( I7TT) . gives 

a h (V*ZV) = X + h(z^v\v) + /i( 7 + 7)A 

-ihtx (v* v Vz. h Vv+ v V h Vv*.zv+ v Vv*. h Vzv) 

modulo 0(h 2 ). From V*V = 1 it follows that the leading symbol of 
V*V equals unity. Since \v | 2 = 1 is the principal symbol, this implies 

Kl + 7) - ihtT ("V ft Vw*.w + w V«*. h Vu) = 0. 
Therefore the expression for the symbol of V*ZV simplifies to 
a h (V*ZV) = \ + h(z_v\v) 

+ ihtr (\ v Vv*. h Vv - v* v Vz. h Vv - v Vv*. h Vzv) 
modulo 0(h 2 ). Using h \J zv = \ h Vv+ ft VA .v we deduce ( 148ft . □ 

Denote ^phg the class of /i-independent pseudo-differential operators 
A with polyhomogeneous symbols, a ~ Sj<m a i> a j homogeneous of 
degree j. When regarded as an /i-dependent operator, A G \l/ m ' m has 
the classical symbol Ylj<m^ l ~^ a j- m ^ ne nex t lemma, following |PV99j 
and [SteOO] . we use this relation to conjugate the scalar operator con- 
structed in Lemma [T9l into hP — 1, where P is h- independent. 

Recall that f2 1//2 -> X denotes the bundle of half-densities. 

Lemma 20. There is a selfadjoint operator P G \l/ ph (X; fi 1 / 2 ) wift 
principal symbol p, and an operator A G \l/ ' /rom half-density sections 
to scalar functions, elliptic nearH, such that A*V*(Z+B)VA = hP—1 
in a neighbourhood o/E. T7ie subprincipal symbol of P equals, on E ; 

Psub = (if|i>) _1 ( Re(,2_t;|t;) + Imtrff^V^-^Vv)) 
(50) \ 

+ Imtr h Vp-"Vf*-f- 

Here z denotes the radial derivative of z. If instead of v another unit 
section v = e tLp v of ker z -> E is used to define V , and thus P, then 
the principal symbol of P remains unchanged, whereas the subprincipal 
changes to p su b = p su b + {p, V 9 } onT,. Here {p, <p} denotes the Poisson 
bracket. 

Obviously, P is elliptic and bounded from below. 

Proof. The radial derivatives of p and of A = (zv\v) are, at E, equal 
to 1 and (zv\v) > 0, respectively. Therefore, near E, Oq\ = p — 
1 for some a G C°°, a > 0. Set Z = A*V*{Z + B)VA , A = 
Op h (ao). Choose Pi G \E f1,0 (formally) selfadjoint with leading symbol 
p — z/itr^V ft Vp)/2. The selfadjoint operators Z and Pi — 1 have the 
same principal symbol, p — 1. Therefore, the imaginary parts of their 
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leading symbols are equal. It follows that the principal symbol go of 
Q = Z — (Pi — 1) G equals, on E, = (zv\v)~ l times the real 
part of the coefficient of h in (HB1) . 

Define p G C°°(T*X \ 0), homogeneous of degree 0, and r_i G S^ 1 
such that go — Po + 2(p — l)r_i holds in a neighbourhood of E. Then 

Op ft (r_i)*)£(l - /i P/l (r_!)) = A + hP - 1 + 

where Po is selfadjoint with principal symbol po- Proceeding induc- 
tively, we obtain selfadjoint operators Pj G \I/ JI ' with classical symbols 
such that, for N < 1, 

(1 - hR* N )Z(l - hR N ) = h h ~ j Pj ~ 1 + h ~ N QN, 

N<j<l 

where Qn G ty N '°, R N g \l/ _1,0 . Therefore, there is an /^-independent 
operator P G ^ hg such that (1 - hR*)Z(l - hR) = hP - 1 near E. 
Moreover, P = Pi + hP modulo ^/~ 1, ~ 2 . The symbol of P equals 
p — itr^V h Vp)/2 +po modulo It follows from Corollary l29l or 
rather its analogue for /i-independent operators, that p is the principal 
symbol of P and p su b = Po its subprincipal symbol. By construction 
Po — go on E. Formula ( )50l) follows from the formula for go mentioned 
earlier. 

Note {p, if} = tr ( v Vp- h V<f- h Vp- v Vf)- The last assertion of the 
lemma follows from fl50|) . using v* v \7 z.v = ^V^o- □ 

Proof of TheoremUi The following assertions hold microlocally in a 
neighbourhood of E. It follows from lemmas [19] and [20] that, if A~* 
denotes a parametrix of A*, (Id -R*)Z(ld -R)VA = VA~*(hP - 1). 
Define J h = (ld-R)VA and J h = {Id-R^VA^*. We have J h , J h G 
* 0,0 , J h - Jh e Moreover, JJJJ h is elliptic. By definition of Z, 

TBhJh — Jh{P — h^ 1 ), where R>h is the Dirichlet parametrix given in 
Proposition \T7\ Combining the results in section [6] with lemmas [19] and 
1201 the theorem follows. □ 

9. Construction of Quasimodes 

Given P of Theorem [1] we associate to the sequence of positive eigen- 
values of P a sequence of quasimodes of Lt- We follow |Ste00l sect. 4], 
differing in some details, however. 

Let P, Bh, and Jh as in Theorem [1] Assume given a sequence of 
quasimodes, (fij), with almost orthogonal quasimodes states: 

(51) Pfj - fijfj = O c ~{hf), (fj\f k ) - S jk = 0((hj + h k )°°), 
fj G L7°°(X; n 1 / 2 ), 0<Hj< ->• oo, h, = fij 1 . 
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We define quasimode states for the traction-free boundary problem. 
By Theorem [T] the traction tj = TB h .J hj fj = Oc°°{h°°)- Choose u'j = 
Oc^iti**) satisfying A (0)d r u'j\x + tj = and u'j[x = 0. Define Uj G 
C?(M;CTM), 

(52) u j = hJ 1/2 (B h .J h .f j +u' j ). 
By Theorem [TJ 

(53) L Uj - \i)puj = O a oo{hf), Tuj = 0, 

and ||mj||l 2 = C(l)- We can assume that the Uj are supported in 
a given neighbourhood of X. Using the ellipticity of L, we deduce 

\\ Uj \\ H 2=0(hj 2 ). 

To go from quasimodes to eigenvalues or, in scattering theory, to 
resonances, it is desirable to be able to decompose the quasimodes 
into well-separated clusters. In addition, the quasimode states of each 
cluster should be linearly independent, and remain so after applying 
small perturbations. 

Proposition 21. Let the assumptions of Theorem U\ hold. Assume 
given quasimodes [ij = hj 1 > of P as in flBTl) . and define Uj as in 
(|5"2"|) . Then (|5"3"|) holds. Let m > dimX. There exist 5 > and a 
covering of {fij} by a sequence of intervals [ak,bk] C M +; such that 

b k + 25b k - m ~ dhaX < a k+1 , b k -a k < b k m . 

Let Wj E H 2 (M; CTM) such that, for some N > 0, 

IK-Uh* = 0(hJ 2 - N ), Wj - Uj = L2 (h 2dimX+N ). 

Then, for large k, {ujj}a k <^<b k is linearly independent. 

Proof. Property fl53|) is clear by the arguments already given. 

It is well-known that a quasimode sequence (151]) is asymptotic to a 
subsequence of the sequence of eigenvalues of P. The latter satisfies 
the Weyl asymptotics. Hence we have a Weyl estimate j < C/j, dimX . 
It follows that every interval [a, b], 1 < b, of length > L has a subin- 
terval of length > Lb~ dimX jC which contains no quasimode fij. The 
existence of intervals [a k , b k ] having the stated properties follows from 
this observation. Compare |Ste99l Proof of Theorem 2] . Define the set 
of indices of the fc-th cluster: I k = {j ; \ij G [ak, b k ]}. 

Choose a left inverse K h G *°'°(X; CT X M, fi 1 / 2 ) of J h , K h J h = Id 
at S. Since J^Jh is elliptic at E, Kh is readily found. 

1 /2 

Denote 7 : v >->■ v\x the trace map. By (T521) . h-' ^Uj = Jhjj + ju'j- 
From (J3Q) it follows that WF hj fj C E. Therefore, 

h) l2 K hjl Uj = fj + Oc~{hY). 
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By the remark after Lemma [26] we can assume that there exists a 
constant C such that for all j, £ G l k , k G N, 

\\Kh e - -K/iJl 2 ^l 2 < Cb k \ht — hj\. 

Using bk\li£ — hj\ < b k a k 2 \fii — < a k 2 b k m+l , it follows that 

h] /2 \\{K he - K hj ) 1Uj \\ L i = 0(b k m ), Jje h, 

if k is sufficiently large. The assumptions on Wj imply \\vjj — Uj\\n^ = 
0(h 1 j +dimX ). Here we use the estimate \\v\\ 2 m < C\\v\\ L 2\\v\\ H 2. Apply- 
ing the trace theorem, W'jWj — 'jUjW^ = 0(h^ +dimX ). Summarizing the 
estimates, we have shown that, for some e > 0, 

\\h] /2 K hm - ftWv = Om +dhaX ), j, £ g I k . 

Because of almost orthogonality of the fj and the Weyl estimate, we 
can apply |Ste99t Lemma 4]. We obtain, for every I G I k , the linear 
independence of {K he -fWj}j € i k when k is sufficiently large. Since K^pj 
is linear, also {wj}j £ j k is linearly independent. □ 

Proof of Corollary^ We apply Proposition [211 with fij f oo the se- 
quence of positive eigenvalues of P, counted with multiplicities, and 
{fj} a corresponding orthonormal system of eigenvectors. Fix m > 
dim X. Let [a k , b k ] be the intervals, clustering {fij}, given in the propo- 
sition. The quasimode states defined in (IB^|) belong to the domain of 
the selfadjoint operator L T . Let n k denote the spectral projector for L T 
of the interval [a' k , b' k ], where a! k = a k -5b k m - dimX , b' k = b k +5b k m ~ dimX . 
The intervals [a' k , b' k ] are pairwise disjoint. Set Wj = n k Uj if fij G [a k , b k ). 
A well-known argument, using the spectral theorem, gives 

5V w - 2dimX IK - «j\\h < W(Lt - ^)u 3 \\h = 0{b k °°) 

if fij G [a k ,b k ]. Since L T is elliptic, we have ||iUj||#2 = 0(fi 2 ). Now 
Proposition (2T], with = 0, implies that, for k sufficiently large, the 
rank of n k equals z{j ; /ij G [a k ,b k ]}. Hence an increase by n of Np 
over [a k , b k ] leads to an increase > n of Nl t over [a^.,6y. Taking into 
account the widths of the intervals, the corollary follows. □ 

Remark. The foregoing arguments also apply to give lower bounds for 
the counting function of resonances. In this case, ir k is the projector 
onto the space of resonant states which correspond to resonances in 
rectangles [ak, b k ]+i[0, s k ]. To satisfy the assumptions in Propositionl2Tl 
for Wj = 7i k Uj, one establishes resolvent estimates. See [SV96j . [TZ98J, 
[Ste99j . and |Ste0 0j. for ways from quasimodes to resonances. The 
clustering method was developed in this context, |Ste99j . to handle 
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multiplicities appropriately. Resolvent estimates for anisotropic elastic 
systems are given in |KN00j . 



10. The Isotropic Subprincipal Symbol 



In this section we assume that the elastic medium is isotropic. We 
evaluate the subprincipal symbol of P, p su b, starting from the general 
formula (150]) . 

We continue with Example El referring to the notation introduced 
there. The kernel bundle ker z is a line subbundle of V, the subbundle 
of CT X M spanned by u, £ = £/|£|. Abbreviate ([19]) and (1151) as follows: 



Ci 

<2 



<2 



(iq) 



ii 



IK 2 \ 



-ZAC12 
^22 



It will be convenient to use the slownesses relative to the Rayleigh 
wave speed, a s = c r /c s and <j p = c r /c p . Then t = a 2 s , ut = ai on 

We first show how to evaluate (z-v\v), v G ker z, z_ as in (1451) . 



^ = { c r|£| = 1}- Moreover, we set r s = (1 



Lemma 22. Set K = (iq) n . Define Yj by (JHSD, and (EH]). Let 

X = (xjk) the selfadjoint 2x2 matrix which is the unique solution of 



(54) XK + K*X = -2Y 1 

Let v — V\V + G ker 2. T/ien 



^ - y* + r 3 



I 3 ■ 



2 Re(z_f |f ) = £ii|t>i| 2 + X22\v2\ 2 + 2 Rexi2fif2- 

Proof. Set x = z_+z*_ . Then 2Re(^_t'|t') = (xv\v). By (j45j) . x satisfies 
the uniquely solvable Sylvester equation x(iq) + {iq)*x = iy + (iy)*, 
where y equals the right-hand side of (]45]) . Since q leaves V and V 1 - 
invariant, X = (x)n = (xjk) is the unique solution of (1541) provided 
the right-hand side of the equation equals (iy + (iy)*) . The latter 
holds if 

Y 1 = (ti{S)z + d r z) u , Y 2 = (oi_ig) lx , F 3 = {itr v V(iq)*.a h Viq) n . 

Observe that the term of P5j) drops out because of the skewness 
of (z<22_)ii. In the following we derive formulas for Yj. 

The basis vectors v and £ do not depend on r. Therefore, (d r z)\\ = 
d r (z)n. We obtain 



(55) 



Y x = tr(S) 



Ci ~i& 

<2 C3 



+ 



d r Ci -id r ( 2 
id r ( 2 d r ( 3 
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Using Lemma [12] and the remark following it we obtain a formula 
for (ai-)n. Clearly, (a 1 _iq) u = (ai_)n(ig)n. We derive 



(56) 



Y, 



(t VA) n tr S + ( A + n) (1 SO \ [ i«2i ^22 
It remains to determine I3. Fix an orthonormal frame (77^-) of T%M, 
r/i — u, r]2 — i- To compute the contraction we use the frame (T)j)j>2 
of T*X, and the dual frame. We compute derivatives of 

iq = \£\Vl - t(Id -1/ ® 1/ - 1 <g> 

Set Sjk = (Sr}j,r}k). A calculation using h \7 v = S and ^V^ = gives 
( h V j iq) u = h Vj(iq)u + s^b^M, j > 2, where 

(ut-b)y/T=t 
i{ut — t) 



M 



21 



{ut-b)y/T=t 

Regard the coefficients Kj k as functions of c s ,c p , |£|. Then h Vj(iq) 11 = 
(rjj, Vc s )K s + (r/j, Vcp)K p , where K s and denote the partial deriva- 
tives of (ig)n with respect to c s and c p , respectively. In particular, 

Define »i = [(^ — 6)^/1 - £, — i(b— ut)). The row fc > 2 in ( h Vj 
equals Sj^b^l^Wi. 

The vertical derivative of a function k which, when restricted to a 
fiber depends only on |£|, is given by its radial derivative: 

(57) v V n K=\£\-%ri)k. 

A calculation using v \7 u = and "VC = — £ <8> gives 

("Viigjn = "Vi(ig)n = \^%jK, j > 2, 
where we have set 

Ku —iK-13 

Define w 2 = [i(b - *), y/l-ut - VI - t\. The row k > 2 in ("Viig) 21 
equals Sji t b~ 1 W2- 

Denote A = (a) n = diag[A + 2jj,,fi\. Note that (0)22 equals \x times 
the unit matrix. Summing over j > 2 we derive 

F 3 = i K*A(\i\- 1 (i,Vc s )K s + Ifr 1 ^, Vc p )if p + 8226"^) 

+ i>6 -2 |f Ktr(S') - S22M <8> iui, 

evaluated at S. □ 
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Denote v the unique unit section of ker(^ — A ) satisfying (£\v) > 0, 

w = 7" 1 «2^+(Ci-A )e), 

7 > such that \v\ = 1. We compute the w-dependent terms in the 
right-hand side of ( l50i) . 

Lemma 23. On E, Imtr h \7p. v \7v*.v = 0, and 

I67 2 Imtrty^V^V^) 
(59) = m 3 /U - 1 c r V s 6 (4 - a s 2 )(2 - a 2 )(2r s C 3 - (2 - a 2 )^)^ 

+ 2m 3 c r a s 6 (2 - o?)(5^ - 4 - a s 4 ) tr'(S), 

where tr'(S) = tr(5) — S22, S22 = (£, anc ^ m — 

Proof. Set 71 = (2/7 and 72 = (Ci — A0V7. We continue to use the 
frame (77^- ) . For j > 2 we have 

Note that ^Vj v*.v is real. Hence Imtr h \7 p. v \7 v* .v = 0. We need the 
vertical derivative of z. To compute it we proceed in the same way as we 
did when computing the derivatives of iq. Recall that z equals Id 
on V ± , where C" 1 = A*|f |V1 = *• We obtain { v V j z) ll = v Vj(z) n . 
Moreover, the column k > 2 in ("Vj2)i2 equals <5jfc|£| _1 times the 
transpose of the row vector [—1(2, (3 — C" 1 ]- We get 

Im v* v Vj z. h Vj v = 71 Re v* v Vj z.Srjj 

= 71(72 B VjCs- 7i "ViC2)s 2 j 

+ 7l|er 1 (72(C3 - C X ) - 7lC 2 )^(l - Sy). 

Summing over j > 2 we obtain 

7 2 Imtr(^V^V^) = CiCa'Vrf Cs - C^V^G 

+ c r C 2 (Ci(C3-C ± )-C2 2 ) tr'(S). 
The first term on the right equals 

m 2 \i\- l s 22 {2b - ^(tVT^Cs - (2b - *)C a ). 
Moreover, using the definition of b, we calculate 

C 3 - C X = m(v / T^ - v 7 ! 3 *)- 
Using (|2TJ), 46 = t(4 - i), we derive ([59]). □ 

The restriction to X of the radial derivative of the eigenvalue Ao = 
(zv\v) = clq 2 (p — 1) equals A = (zv\v) = Oq 2 because p = 1 on E. 
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Lemma 24. On S 7 

7 2 A = m 3 a!?(4 - (j 2 )t s (t p /t s + c s t s /c p t p + a 2 - 2), 
where m = fi\^\/b. 

Proof. The section w = i&v + Cil; equals jv on E. Therefore, to second 
order on E, 7 2 A = (zw\w) = Ci det(z)n. Inserting (|20|) . 

7 2 A = m6tv / T 3 t(4 v / (l - t)(l - ut) - (2-t) 2 ). 

Recall p = c r |£| = a s t~ 1/2 , £ = {t = o" 2 }. The rule of de l'Hospital 
gives 

, 4:J(l-t)(l-Ut) - (2-t) 2 t n, . , 9 \ 

I™, — ,_ 1/2 \ = (Tp Ts + C S T S /C P T P + a] - 2 . 

Summarizing, the formula for 7 2 Ao = 7 2 Ao/(p — 1) follows. □ 

Inserting the formulas of the lemmas of this section into the general 
formula fl50|) for the subprincipal symbol of P we obtain a formula for 
the subprincipal symbol in the isotropic case. 

Proposition 25. Denote X = (xjk) the 2x2 matrix solving (|54|) . Se^ 
N = t s (t p / t s + c s T s /c p T p + a 2 — 2) . Lei P fre £/ie operator of LemmaWJ\ 

determined by the unit section v of ker z having positive £ component. 
The subprincipal symbol of P is given as follows. 

IQNpsnb = (cr/2fj)(x u (2 - a 2 ) 2 + 4x 22 (l - a 2 ) 2 + 4Imx 12 (2 - ct 2 )t s ) 

+ / i" 1 c 2 (2 - a 2 ) (2t s ( 3 - (2 - a 2 )( 2 ) (S& |) 

+ 2c r (4 - a 2 y\2 - a 2 )(5a 2 - 4 - o*)(ix(S) - (SH)) 

Proof. OnE,ti) = 7D = (mi/2) (z(2-t)z/+2 v / l - *Q • Using lemmas[22] 
and [241 we calculate 16iVRe(z_if |u>)/ 7 2 Ao. The result is the first term 
on the right-hand side of the claimed formula. Similarly, we obtain the 
other terms combining the lemmas [23] and [2H □ 

The constituents of the above formula for p su b are curvature and 
velocities (Lame parameters), assumed known. It seems difficult to 
analyze the formula further unless it is specialized to particular cases. 
However, it should be noted that the formula allows explicit numerical 
evaluation of p su b- Therefore it can be used when solving transport 
equations for Rayleigh wave amplitudes numerically with a (seismic) 
ray tracing program, say. Formulas for the amplitudes of Rayleigh 
waves were given by Babich and Kirpichnikova in [BK04] . 
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Appendix A. Geometric Pseudo-differential Calculus 

Pseudo-differential operators on manifolds are usually introduced by 
reducing to the euclidean case via partitions of unity, |Hor85j 18.1], 
[EZ] . The principal symbol of a pseudo-differential operator is invari- 
antly defined. If the operator acts on sections of the line bundle of half- 
densities then there also is an invariantly defined subprincipal symbol, 
[Hor85l Theorem 18.1.33], [SZ021 Appendix]. 

In the body of the paper we explicitly track, down to the subprin- 
cipal level, symbols of operators acting between vector bundles. To 
achieve this we use Sharafutdinov's geometric pseudo-differential cal- 
culus |Sha05a| ISha05bj . The purpose of this appendix is to recall this 
calculus, presenting a semiclassical variant. Since we have to refer, in 
the main part of the present paper, to proofs of the calculus, we give 
a rather detailed presentation. The calculus depends on a symmetric 
connection of the manifold and on metric connections of the (Hermit- 
ian) bundles. We make the stronger assumption that the manifold is 
Riemannian and that the symmetric connection is the Levi-Civita con- 
nection. The important features of the calculus are a symbol isomor- 
phism modulo order minus infinity, and complete symbol expansions 
for products and adjoints given solely in terms of geometric data. Us- 
ing connections to develop a pseudo-differential calculus and to prove 
the existence of a complete symbol isomorphism was done earlier by 
Widom, [Wid80j. This was further developed by Pffaum who gave a 
convenient quantization map from symbols to operators, [Pfl98j . Shara- 
futdinov gave symbol expansions in terms of geometric data. 

Let X a compact Riemannian manifold without boundary, dim X = 
n. The exponential map, exp, of the Levi-Civita connection defines 
a diffeomorphism, (x,v) i— > (x,y) = (x,exp x v), between a neighbour- 
hood of the zero-section of the tangent bundle T = TX and a neigh- 
bourhood of the diagonal in X 2 . In the proofs of the propositions 
below we need the following properties of exp. In local coordinates the 
exponential map satisfies 

(60) (exp x . vY = x i + v i - r jk (x)v j v k /2 + 0(\v\ 3 ), 

where T l - k denote the Christoffel symbols. Normal coordinates centered 
at x satisfy (exp x t>) 1 = v l . There exist < r < R < inj(x), the 
injectivity radius of X, such that the equation 

(61) exp exp ^ z = exp x w 

defines, for every v G T x = T X X, \v\ < R, a diffeomorphism w i— > 
z = z(x,v,w) from an open neighbourhood of the origin, contained 
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in {|iu| < R} C T x , onto the ball {\z\ < r} C T y , y = exp x v. 
This map is used below to change variables of integration. Obviously, 
z(x,0,w) = w. A computation in normal coordinates centered at x 
shows that 

(62) (z' w )- 1 z = w-v + 0{{\v\ + \w\) i ) asv,w^0. 

Recall, from section El the notation for segments and for parallel trans- 
port maps. In local coordinates, 

(63) {T^^wy = w*- V) k {x)wh k + 0(\v\ 2 ). 

Let E -y X and F -* X be Hermitian vector bundles with metric 
connections. Recall from section [3] the definition ( T2"4"l) of horizontal 
derivatives and the definition of vertical derivatives. A C°° section 
a of the bundle n* Hom(E, F) -> T*X is called a Hom(£', F)-valued 
symbol of order m G R, a G S m = S m (T*X; vr* Hom(£, F)), iff for all 
nonnegative integers j and £, 

sn P (l + \^\y~ m \( v Vy( h vYa(x,0\<oc. 

These are the usual type 1, symbol estimates. The symbol space 
S m is a Frechet space. The space S m ' k = S m ' k (T*X; n* Hom(E, F)) 
of /i-dependent symbols of order m and degree k is the Frechet space 
of families G S m such that {h k ah ; < h < 1} is bounded in 
S m . We call a G S' ni ' fc classical if there exists an asymptotic expansion 
a ~ ^ . h?~ k dj with /^-independent symbols a 3 - G S m ~K 

In the following lemma we define, in a semiclassical setting, the quan- 
tization of symbols according to Sharafutdinov's geometric pseudo- 
differential calculus. We relate this definition of /i-pseudo-differential 
operators to the definition in the euclidean situation. For semiclassical 
analysis, in particular, for the class \l/ m ' fc = Op ft S m > of /i-pseudo- 
differential operators, including mapping properties, and for frequency 
sets (/i-wavefront sets), refer to |Ger88j . [Tv^8] . [DS99] . [SZ02] . [EZ] . 
The class of negligible operators, vl/ -00 - -00 , consists of /i-dependent op- 
erators whose Schwartz kernels are C°° with 0(h°°) seminorms. 

Fix xo £ C°°(TX), real-valued, |u| < r on the support of xo(%,v), 
such that Xo = 1 m a neighbourhood of the zero-section in TX. 

Lemma 26. Let ah G S m ' k be a Hom(E, F) -valued symbol. Then 



(64) 




A h u h {x) = (2irh)- n / / e-^'^xo^v) 



T* TV 



a h (x, r])T [x ^ exPxV] u h (exp x v) dv dr], 
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defines an h-pseudo- differential operator Ah G \& m,fc (X; E, F). Given 
a point x there exists a geodesic ball U centered at x, and a symbol 
a\ G S m,k such that, for Uh compactly supported in U , 



where y = exp^f and y' = exp x w. Moreover, at x, = modulo 
gm-2,k-2^ Every h-pseudo-differential operator is, modulo negligible 
operators, of the form f )64|) . 

The measures in (1641) are the normalized Lebesgue measures of the 
euclidean spaces T x and T*. 

Proof. We shall drop the subscript h from the notation. Fix x £ X. Let 
U denote a geodesic ball with center x and radius < R. In the following 
we assume that the support of u is a compact subset of U. In (IMj) we 
replace the variables x,v,rj by y,z,(- Next we change variables in the 
integral Au(y) such that the domain of integration does not depend on 
y. Set y = exp x v. Define z = z(x, v, w) by flBTl) . Using the symplectic 
map (w,$) i — y (z, C), C — t (. z ' w )~ 1 $i we S e t 



ip = (C,z) = (z'^^z). Since z = if and only if v = w, we have 
ip{v,w,-&) = (ip(v,w)$,w — v). Here if) = ld+0(\v\ 2 + |w| 2 ) by ( 162]) . 
Decreasing the radius of U and making the linear change of variables 

9 = i[)(v, w)$, we get 



Ji(v,w) = 1 + 0(M 2 + \w\ 2 ). It follows that A restricted to U is a 
/i-pseudo-differential operator of class fy m > h . As it stands the symbol 
depends on v, 9, w. Using the standard symbol reduction procedure we 
obtain a u (exp x v, 9). Moreover, the asymptotic expansion implies that, 
at v = 0, a u - a G S m ~ 2 ' k - 2 . 

Note that Au(y) = if the distance between y and suppw is > r. 
Using a partition of unity, we infer that the class of operators given by 




T* T a 




where the kernel K is given by 
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( |64|) equals the class of /i-pseudo-differential operators with Schwartz 
kernels supported in small neighbourhoods of the diagonal. □ 

Standard arguments show that up to a negligible operator Ah = 
Op h (dh) does not depend on the choice of the cutoff xo- The space 

V m > k {X; E, F) = Op h S m ' k + $- 00 '- 00 . 

is the space /i-pseudo-differential order m and degree k. We denote the 
geometric symbol ah(A h ) = a h . 

Remark. Let Ah = Op h ( a h) £ ^ /0 '°- Then Ah is L 2 bounded, uni- 
formly in h. Assume, in addition, that ah depends differentiably on h 
with dh<ih £ S 0,0 . Changing variables in (1611) from 77 to £ = 7]/h, we 
obtain A hl - A ho = h' 1 Op h (b h ) dh, where b h G S* ' , b h (x,r}) = 
hdhah(x,rj) + v \7 ri ah(x,ri). This implies the following useful Lipschitz 
estimate: 

\\A hl - A h J L 2^ L 2 < Ch^hx -h \ if h < h, 

where || Op h (°h) \\l 2 ->-l 2 < C < 00. The assumption holds if a h is 
classical and given as a Borel sum. 

In the following, we often suppress from writing the /i-dependence 
of symbols, operators and distributions. Moreover, when dealing with 
integrals like fl64l) . we move, without explicitly writing this, the x- 
dependency from the domain of integration into the integrand using 
arguments as in the proof of the lemma. 

Lebesgue measure dv on T X X and Riemannian volume are related 
by f f(y)dVx(y) = f f(exp x v)J (x,v)dv, y = exp x t>, with Jacobian 
J = 1 + C(M 2 ) at v = 0. Let A = A h as in (JM}- The Schwartz kernel 
K A of A, 

Au{x) = / K A (x,y)u(y)dV x (y), K A (x,y) e Eom(E y ,F x ). 
Jx 

equals in a neighbourhood of the diagonal a partial Fourier transform 
of the symbol, 

(65) K A (x,y) = (27rh)- n J e*to a »- 1 M h a(x,ri) d V tP(x,y)r^ y] . 

Here ip(x, y) = Xo( x i v )/ Jo(x, v), y = exp x v. The symbol a is recovered 
via the inverse Fourier transform: 

(66) a(x,0 = J e^ h ( Xo Jo)(x,v)K A (x,exp x v)T { l Px ^ x] dv 
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modulo S~ oc, ~ co . The correspondence between an operator A = Op h (a) 
and its full symbol a, named the geometric symbol of A, defines the 
complete symbol isomorphism 

^ m ' fc (X; E, F) Z^- 00 -- 00 S m ' k (T*X- Rom(n*E, Tr*^))/^ 00 ^ 00 . 

The geometric symbol can also be computed by applying the operator 
to suitable testing functions as follows. 

(67) a(x,Z)s = ^(e'^'^xo^exp^y)^^)!^. 

Here A y means that A acts on functions of the variable y. In particular, 
in case E — C, the geometric symbol is obtained at the center of normal 
coordinates when A is applied to eJ&i* 3 l h and evaluated at x J =0. 

We derive symbol properties and expansions using the method of 
stationary phase: 

(det(H/2nih)) 1/2 J e i( ^ x))/h a(x) dx 

= exp (2~ 1 ih(H- 1 d, d)) (e lp(x)/h a(x)) \ x=Q 

= E { ^(H-%dy(e^y h a( X ))\^ + O(h N ), 

j<3N 

if <p G C°°, real-valued, <p'(x) — iff x = 0, H — ip"(Q) non-singular, 
<p(0) = 0. The remainder p(x) = (p(x) — (Hx,x)/2 vanishes to third 
order at x = 0. The expansion has the advantage, when compared to 
that obtained using the Morse lemma, of giving an efficient algorithm 
for computing the asymptotic series. See |Hor83l Theorem 7.7.5] where 
the expansion is arranged in powers of u^ 1 = h. 

We are mainly interested in the leading symbols of operators. We 
call the residue of a in S m ' / S m ~ 2 ' k ~ 2 the leading symbol of an operator 
OP/i( a ) £ \fr m < k . The principal symbol is, of course, the residue in 

Proposition 27. Let A = Op h (a) as in with geometric symbol 
a e S m ' k . The formal adjoint A* 6 ^ m < k (X; F, E) has the geometric 
symbol 

(68) 6 = o*-i^tr ,, v' l Va* mod T" 2 ^ 2 . 

If a is classical then so is b. 

Notice that 11 V h V a* is a section of vr*(Hom(F, E) <8> T ® T*). The 
trace is taken of the T Cg> T* part. 
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Proof. The formal adjoint of A is denned by 

/ ( Ul (x)\Au 2 (x)) F dV x (x)= [ (A* Ul {y)\u 2 (y)) E dV x (y). 
Jx Jx 

The Schwartz kernel satisfies Ka*(x,v) = Ka(v,x)*. Recall that par- 
allel transport preserves inner products. It follows from (l65l) that 

K A .(x,y) = (2nh)- n [ e^^ 1 '""i^afa vT d V ^(y,x), 



T* 
V 



and K A *(x,y) = if the distance between x and y is > r. Set y = 
exp x v. Define z G T y by exp y z = x. After a linear change variables 
from rj G T* to ( = t (exp' x (v))r] G T* we have 

K A .{x,y) = (2nh)- n [ e^V^o^, rj)* d( ^p(y,x)/J 1 (x,v), 



with Jacobian Ji(x, v) = 1 + 0(|v| 2 ). Define 

b(x,0= [ e^/ h ( X oJo)(x,v)K A *(x,y)T^ x] dv. 



Inserting K A * we have 

(69) 6(x,0 = {2ixh)- n J Je iv/h UdCdi 



T x T* 



where 



<p = (Z,v) + (r),z) = -((-£,v) + ((,$), 

5 = T [^ y ]a(y,vT^x] = T fr™} F ' E)a (y,vY, 

J = xo(x, v)J (x, v)ip(y, x)/ J x (x, v) = 1 + 0(\v\ 2 ), 

and $ = <$>(x,v) = exp^f)" 1 ^ + v. A computation in normal coordi- 
nates centered at x shows that $ = 0(\v | 3 ) as v — > 0. If ip'^ = then 
z = 0, hence v = 0. It follows that the critical points of ip are defined 

by v = o, c = e- 

Apply the method of stationary phase to (16"9"j) and deduce b G if> m ' fe . 
Moreover, the following asymptotic expansion holds: 

( 70) i-E^-*.*^**"*)!^- 

Differentiation of the exponential factor brings out a non-zero factor 
only if it comsumes at least three derivatives with respect to v and 
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at most one derivative with respect to (. It follows that the sum is 
asymptotic. Moreover, b is determined modulo S m ~ 2 ' k ~ 2 by the terms 
in the asymptotic sum with j < 2, b = a* — ih{d^, d v )a. Observe that 

r [Lcx P ^] ° exp!» = Wt. +0(\v\ 2 ) as v ->• 0. 

It follows that 9„a|^_ = k V«*(s, £). Hence b = a* — ihtr v \7 h \7a*. 
The Schwartz kernels of Op fe (6) and A* are equal in a neighbourhood 
of the diagonal. Therefore A* - B G □ 

Proposition 28. Let A G ^ mA ' kA (X; F,G) and B G ^ mB ' ks (X; E, F) 
with geometric symbols a and b, respectively. Set k = k^ + ks, rn = 
ttia + mB- Then AB G ty m,k (X; E, G) with geometric symbol 

(71) c = ab-ihtr( v Va h Vb) 

modulo S m ~ 2 ' k ~ 2 . If a and b are classical then so is c. 

Again the trace is taken of the T Cg> T* part, and the dot terminates 
differentiated expression. 

Proof. Setting y = exp x v , C = AB is given by 



Cu{x) = {2uh)~ 2n //// e- l ^ v)+ ^ z))/h a{x^) 




T x xT*xT y xT; 



| (%, O^cxp, z] u ( ex P y z)) dz dC dv dry. 

Here and in the following we do not write the cutoff factors. Let 
z = z(x, v, w) be the solution of exp z = exp x w. The symplectic 
change of variables (w, i?) h-> (z, C), C — *(< 2 4) _1 '#, preserves the volume 
form. We get Cu(x) = f T K c (x,exp x w)u(exp x w)J (x,w)dw, with 
Schwartz kernel 

Kc(x, exp x w)Jo(x, w) 

= {^h)- 2n J e-^+^^co d(0, v, 77) T[l exp 



_ M" 1 

T* xT x xT* 



c = a(x,r j )r^ E ' F) b(y,C)M(x,w,v). Here M{x,w,v) G GL(E X ) 
denotes the parallel transport in E along the geodesic triangle x — > 
exp x w — > exp x v — > x. It follows that the symbol of C equals 

(72) c(x, = (2nh)- 2n J e*< h Co d(v, 77, w, 0), 
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<f = (£, w) — (77, v) — (C, z). We introduce w — v as a new variable, u>. 
Then (J72D holds with 

<p = -(ri-Z,v)-{&-Z,w) + {&,$), 

c = a(x, ri)r^°^ E ' F) b(y, C)M(x, w + v,v), 

Here $ = to— (^(a;, u, w+v t>, w + w) G T*. By (|02"]) . $ vanishes 

to third order at v — w — 0. Clearly, t> = = z at a critical point of 
ip. It follows that f = w = and r\ = $ = £ define the critical points. 

Now apply the method of stationary phase to ATS]) and deduce that 
c G S"" 1 ^ is a symbol which, moreover, has an asymptotic expansion 



( 7 3) c~ a„> + <3„ s„» 3 (««♦«/»«*) |_,„. fc{ . 

i 

Using that $ does not depend on rj and and vanishes to third order 
at v = w = 0, we infer that the summands with j > 1 belong to 
5 m-2,fc- 2 _ It follows that 

ab — ih{d v a, d v bM) — iha(d$, d w )bM, 

evaluated at the critical point, is the leading symbol of C. Here b = 
r [x^-™] "{ViO- We have d w b = at v = w = 0. This follows from 
(' w = which is a corollary of z = w at v = 0. The derivatives of M 
with respect to v and u> vanish at v — w = 0. Using tT „ , o = 

Id^ +C(|v| 2 ) at if = 0, we derive 

^ = ^[^SKexp^/K)- 1 ^) = 'V&M), 
at v — w — 0. Summarizing the computations, (|71[) follows. □ 



Remark. The proofs of propositions [27j and [28] follow those in |Sha05a 



ISha05b] closely with only minor modifications. Our derivation of the 
asymptotic expansions of the symbols of adjoints and products may 
be somewhat shorter, however. We differ in defining the adjoint with 
respect to the volume element rather than by using half-densities. No- 
tice that the symbol expansions (170]) and (173]) depend only on the given 
symbols and on the geometry. In the formulas (168]) and (ITT]) , we ex- 
tracted the leading symbols. 

For the purposes of the present paper it suffices to assume X com- 
pact. A symbol calculus on general (complete) Riemannian manifolds 
needs to take the injectivity radius into account and handle mapping 
properties more explicitly. 

It is well-known that a pseudo-differential operator acting on half- 
densities has an invariantly defined subprincipal symbol; see [SZ02, 
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Appendix] for a proof in the semiclassical case. We relate the subprin- 
cipal symbol to the leading geometric symbol. Equip the half-density 
bundle fi 1/2 -> X with the inner product (u\v) = u ■ vj dVx, where 
the operations on the right are in the sense of densities. The connec- 
tion given by V Ql/2 dV x 2 = is metric with respect to the Hermitian 
structure of Vt 1 ^ 2 . 

Corollary 29. Let A G ^ m > k (X; VL l/2 ). The leading symbol of A equals 
that of the corresponding scalar operator A G \I/ m ' fc (X) which is given 

~ — 1 /2 1/2 

by Au = dV x A(udV x )■ If the geometric symbol a of A is classical, 
a ~ 2~2j>o h? dj, dj G S ,m_J ; then h~ k ao is the principal symbol of A, 
and 

a su b = h l ~ k (ai +i v Va . h Va /2) 
is its subprincipal symbol. 

Proof. Consider the multiplication operator dV x 2 G \I/ '°(X; C, Vt 1 / 2 ). 

The Hom(C, fi 1 / 2 )-valued symbol tc* dV x 2 is the leading symbol of this 
operator. Note that its horizontal and vertical derivatives vanish. The 
equality of the leading symbols of A and A now follows from Proposi- 
tion E3 

Let a u denote the local symbol of A in a geodesic coordinate chart 
U centered at a given point x. We use normal coordinates centered at 
x. Assume a classical, h k a = a + ha\ + 0(h 2 ). Then a u is classical, 
and h k a u = ao + 0(h). Moreover, it follows from Lemma [26] that 
h k a u = a + ha\ + 0{h 2 ) at x. The subprincipal symbol equals, by 
definition, h 1 ~ k (ai + 2~ l i d 2 ao/ dxjd^j). The horizontal derivative 
in the j-th coordinate direction equals, at x, the partial derivative with 
respect to Xj. The formula for the subprincipal symbol follows. □ 
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